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1 .  INTRODUCTION 


1.1  Summary 

Primary  activities  during  the  third  year  consisted  of  (i) 
expanding  the  laboratory  capabilities  for  processing  and  testing  of 
composites,  ( i i )  conducting  research  in  accordance  with  the  Statement 
of  Work  given  in  Section  1.3,  (iii)  preparing  thirteen  technical  papers 
and  M.S.  theses,  and  (iv)  various  interactions  of  the  faculty  with  the 
technical  community  through  presentation  of  papers,  participation  as 
members  of  technical  committees,  etc. 

Sections  2-4  summarize  the  research  activities.  The  professional 
personnel  associated  with  the  project  and  the  outside  activities  of  the 
faculty  related  to  composites  are  given  in  Section  5.  Papers  completed 
and/or  published  during  the  year  are  reproduced  in  Appendix  A.  A 
brochure  describing  the  student  activities  is  in  Appendix  B. 

We  provide  only  a  brief  discussion  of  most  activities,  as  they 
are  detailed  in  existing  reports  (Appendix  A)  or  theses,  or  are  in 
early  stages  of  development.  However,  the  processing  and  testing 
facility  has  been  greatly  expanded  in  1980,  but  its  increased 
capability  is  not  described  elsewhere;  thus  a  relatively  detailed 
description  is  given  (Section  2.1). 

1.2  Discussion 

The  general  objective  of  the  research  program  is  to  develop 
improved  understanding  of  fibrous  composite  material  as  a  basis  for 
improving  its  structural  performance.  Integral  to  achievement  of  this 
objective  is  the  determination  of  constitutive  equations  and  time  (or 
cycles)  to  fracture  for  complex  loading  patterns  and  hostile  environments 
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such  as  high  relative  humidity.  This  in  turn  requires  a  better  under¬ 
standing  of  the  viscoelastic  deformation  and  mechanisms  of  damage 
growth  and  fracture  in  composites  including  delamination,  cumulative 
damage  via  microcrack  initiation  and  growth,  and  global  fracture  criteria. 

The  importance  of  micromechanisms  of  viscoelastic  deformation  and 
fracture  of  resin  matrix  composites  has  only  recently  become  widely 
appreciated;  for  example,  graphite/epoxy  composites  are  quite  sensitive 
to  processing  and  to  service  environments,  especially  moisture  [l]t 
Furthermore,  because  -...'ing  is  seldom  taken  to  completion  in  the  epoxies 
used  in  graphite/epoxy  composites  [2],  the  degree  of  cure  is  a  processing 
variable.  Physical  aging  (densification)  of  the  epoxy  matrix  may  occur 
in  service,  giving  significant  changes  in  mechanical  properties  [3], 

The  post-cure  cool-down  path  affects  the  magnitude  of  residual  stresses 
in  a  composite  [4].  Moisture  absorption  can  change  the  glass  transition 
temperature,  giving  a  significant  change  in  mechanical  properties  [1]. 

The  rate  of  moisture  absorption  appears  to  depend  on  the  degree  of 
curing  and  physical  aging,  and  moisture  content  may  itself  affect  the 
rate  of  physical  aging.  The  lack  of  consistency  in  experimental 
results  and  micromechanistic  modelling  appearing  in  the  composite 
materials  literature  today  is  believed  a  consequence  of  the  fact  that 
much  of  the  earlier  experimental  research  was  conducted  without  fully 
controlling  all  of  the  important  variables  which  affect  resin  matrix 
composites. 

It  is  anticipated  that  the  continuing  research  program  at  Texas 
A&M  University  described  herein  will  make  a  major  contribution  toward 
more  unified  methods  of  characterizing  behavior  of  composite  materials  by 

♦Numbers  in  brackets  indicate  References  on  p.  48. 
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identifying  and  accounting  for  the  important  variables.  In  view  of  the 
large  number  of  variables  to  be  considered  experimentally,  it  is  essential 
that  the  mechanical  testing  be  guided  by  mechanistic  modelling  to  keep 
the  total  amount  of  experimental  work  at  a  manageable  level.  A  consider¬ 
able  part  of  the  effort  at  Texas  ASM  University  is  being  so  directed. 
Additionally,  some  of  the  theoretical  research  is  leading  to  analytical 
methods  of  predicting  deformation  and  fracture  behavior  of  laminates  for 
use  in  design  analysis  of  those  structures  that  must  withstand  severe 
environments  over  long  periods  of  time. 

1.3  Statement  of  Work 

"a.  Investigate  the  effects  of  cure  cycle  parameters  on  the 
mechanical  characteristics  of  resins,  composites  and  composite 
structural  specimens: 

(1)  Study  the  curing  process. 

(2)  Investigate  the  effects  of  cure  cycle  parameters  on 
physical  aging. 

(3)  Investigate  generation,  relief  and  effects  of  residual 
stresses  in  laminates. 

b.  Investigate  deformation,  damage  growth  and  fracture  behavior 
for  resins,  composites  and  composite  structural  specimens: 

(1)  Develop  and  verify  constitutive  equations. 

(2)  Develop  and  verify  ply  damage  and  delamination  models. 

(3)  Develop  an  automated  structural  material  property 
characterization  system  and  data  base." 
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2.  Processing  and  Testing  of  Resins  and  Composites 
2.1  New  Equipment  and  Software* 

During  the  year  1980,  several  pieces  of  new  equipment  were  acquired 
for  expansion  of  the  experimental  research  laboratory  to  broaden  the 
capability  in  processing  and  testing  of  composite  materials.  The  majority 
of  this  equipment  was  purchased  with  funds  from  the  Texas  Engineering 
Experiment  Station;  supplemental  funds  were  provided  by  AFOSR.  The  follow¬ 
ing  equipment  was  purchased,  delivered,  and  installed  in  the  composite 
materials  experimental  research  laboratory  (McNew  Building)  during  the 
year: 

curing  oven 
vacuum  oven 
process  controller 

cold  stage  for  differential  scanning  calorimeter 

diamond  tools 

humidity  monitor 

digital  thermometer 

laboratory  oscilloscope  and  modules 

closed  loop  servohydraulic  mechanical  test  system 

MINC-11  laboratory  data  acquisition  and  control  computer 

HP-45  desktop  computer 

digital  oscilloscope 

ultrasonic  C-scan  system 

charge  amplifier 

HP  digital  plotter 

HP-85  desktop  computers  (2) 

data  acquisition  and  control  units  (2) 

synthesizer/function  generator 


*Prepared  by  Dr.  K.  L.  Jerina 


5 


The  laboratory  oscilloscope  and  modules  are  for  general  purpose 
laboratory  use  and  the  other  equipment  has  been  incorporated  into  specific 
functional  systems  for  processing  and  testing  of  composite  materials. 

The  capacity  for  processing  and  machining  of  test  samples  has  been 
enchanced  by  the  addition  of  a  vacuum  oven  and  convection  oven.  These 
ovens  are  used  to  cure  and  post  cure  resins,  adhesives  and  composites 
fabricated  from  prepreg  systems.  An  assortment  of  new  diamond  tools  for 
the  Micromech  cutting  machine  has  increased  the  capacity  and  quality  of 
machining  specimens  with  hard  fibers  such  as  glass  and  graphite. 

The  processing  press  has  been  used  extensively  in  glass/epoxy  and 
graphite/epoxy  specimen  fabrication  for  student  and  faculty  research 
projects.  During  the  reporting  period,  a  microprocessor-based  process 
controller  was  installed  on  the  press  and  control  parameters  optimized. 

The  process  controller  makes  it  possible  to  generate  a  wide  variety  of 
temperature,  vacuum  and  pressure  profiles  for  cure  cycles.  New  cure 
cycles  are  easily  programmed  in  an  automatic  fashion.  Heating  or  cool¬ 
ing  rates  of  3°F/min.  are  possible  while  controlling  within  an  absolute 
temperature  reading  of  +3°F.  The  best  accuracy  and  control  for  an 
arbitrary  temperature  profile  is  obtained  at  a  rate  of  3°F/Min.,  although 
the  press  is  capable  of  higher  heating  and  cooling  rates.  The  control 
scheme.  Fig.  1,  consists  of  a  primary  control  loop  on  laminate  temper¬ 
ature  and  two  secondary  control  loops  on  the  upper  and  lower  platen 
temperatures.  The  actual  laminate  temperature  is  blended  with  the 
command  temperature  profile  in  a  cascade  control  scheme  to  command  the 
upper  and  lower  platen  heaters  and  cooling  water  valves.  Optimization 
of  the  controller  was  accomplished  through  an  empirical  process  of  ad- 


Block  Diagram  of  Press  Controller 
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justing  the  individual  loop  parameters  for  gain,  rate  and  reset.  The 
controller  is  interlocked  so  that  an  orderly  shut  down  occurs  on  an 
error  such  as  over  temperature  or  thermocouple  burn  out.  The  press 
controller  has  proven  successful  in  increasing  production  of  laminate 
samples  and  in  improving  the  accuracy  and  repeatability  of  cure  cycles. 
Additionally,  the  ability  to  study  the  effects  of  cure  cycle  variation 
is  greatly  enchanced  by  the  press  controller. 

The  Perkin  Elmer  Differential  Scanning  Calorimeter  II  has  an  extended 
temperature  range  made  possible  through  the  addition  of  a  liquid  nitrogen 
cold  stage.  Also  a  MINC-11  laboratory  data  acquisition  and  control  computer 
has  been  interfaced  to  the  DSC,  Fig.  2.  The  temperature  and  specific  heat 
signals  from  the  DSC  are  now  acquired  by  the  computer  through  a  digital 
input  module  and  analog-to-digital  converter.  Under  control  of  a  real 
time  computer  program  the  specific  heat  of  a  resin  sample  can  be  monitored 
as  a  function  of  temperature.  Once  acquired  by  the  computer,  the  data 
can  be  easily  manipulated  to  subtract  baseline  calibrations  and  then 
permanently  stored  as  part  of  a  data  base  on  floppy  disk  for  later  review 
and  analysis.  The  data  can  be  plotted  on  a  video-graphics  terminal  or 
digital  plotter.  Automation  of  the  DSC  has  increased  the  accuracy  and 
repeatability  of  experiments  involving  glass  transition  temperature  studies 
of  resins  used  in  composite  materials. 

New  instrumentation  has  aided  the  project  in  development  of  the 
piezoelectric  "duomorph"  complex  modulus  gage.  A  charge  amplifier, 
digital  oscilloscope,  function  synthesizer  and  desktop  computer  have  been 
incorporated  into  the  instrumentation.  Fig.  3.  The  task  of  data  acquisition 
and  reduction  for  a  typical  experiment  has  been  reduced  from  several  weeks 
to  less  than  a  day  through  automation.  The  electrical  excitation  and 


Figure  2  Block  diagram  of  differential  scanning  calorimeter  and  laboratory  instrument  computer 


Figure  3  Block  diagram  of  duomorph  instrumentation 
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response  of  the  duomorph  are  monitored  by  the  digital  oscilloscope. 

Under  control  of  a  real  time  program  on  the  HP-85  desktop  computer, 
excitation  frequency  is  automatically  programmed  on  the  function 
synthesizer  and  sampled  values  of  the  gage  excitation  and  response 
are  transferred  to  the  computer  for  amplitude  and  phase  angle  analysis, 

Fig.  4.  Automation  has  reduced  the  data  acquisition  and  analysis  time 
for  the  duomorph  gage, making  it  possible  to  consider  the  duomorph  as  a 
real-time  cure  process  monitor  if  feasability  studies  are  successful. 

Management  of  a  mechanical  property  data  base  and  the  analysis  of 
experimental  data  has  been  enchanced  by  an  HP-45  desktop  computer  and 
HP  9872B  digital  plotter.  Fig.  5.  The  system  is  designed  to  allow  analysis 
of  experimental  data  stored  in  a  data  base  with  a  powerful  computational 
capability,  video  graphics,  digital  plotter,  hard  copy  graphics,  and 
communications  to  other  laboratory  instrument  computers.  This  system 
has  been  installed  during  the  reporting  period  and  software  is  currently 
being  developed  to  bring  the  system  up  to  its  full  usefulness  and  potential 
as  an  interactive  graphics  system.  Fig.  6  shows  the  graphics  capability  of 
the  system  for  a  [+45°] 2S  graphite/epoxy  creep  experiment.  The  greatest 
utility  of  this  data  analysis  system  will  be  to  allow  interactive  access 
and  analysis  of  data  base  information  with  high  quality  graphics  capability. 

Data  acquisition  for  the  five  channel  temperature/humidity  creep/recovery 
system  has  been  accomplished  with  an  HP-85  desktop  computer  and  an  HP-3497A 
data  acquisition  and  control  unit,  Fig.  7.  Also,  a  multi -point  digital 
thermometer  and  humidity  gage  monitor  the  environmental  conditions  in  the 
test  chambers.  Creep  and  recovery  data  are  acquired  by  the  data  acquisition 
unit  under  program  control  from  the  computer.  The  program  samples  two 


Figure  4  Data  analysis  of  duomorph  excitation  and 
for  amplitude  and  phase  angle  (theta) 


Figure  6  Creep  of  [+45] s  graphite/epoxy 
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Figure  7  Block  diagram  of  five-station  automated  creep/recovery 
tempera tu re/ humidi ty  apparatus 
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channels  of  strain  data  from  each  creep  station  at  logarithmic  time  intervals 
The  data  are  tabulated,  Fig.  8,  and  stored  as  a  file  as  part  of  the  lab¬ 
oratory  data  base.  Implementation  of  this  hardware  configuration  has  allowed 
more  cost  effective  and  efficient  use  of  laboratory  instrumentation. 

Capability  for  non-destructive  inspection  of  composites  has  been  added 
to  the  laboratory  during  the  reporting  period.  A  Testtech  scanning  bridge, 
tank,  and  plotter,  Panametrics  transducers  and  tranducer  conditioner,  and 
an  HP  delayed  sweep  oscilloscope  give  the  laboratory  an  NDT  capability 
for  composite  laminates.  Fig.  9.  A  selection  of  transducers  in  the  range 
of  2.5  to  25  MHz  allow  inspection  of  graphite/epoxy  and  glass/epoxy  laminates 
Initial  tests  of  the  system  show  that  it  is  possible  to  detect  small  density 
changes,  such  as  an  extra  roving,  in  cured  laminates.  The  system  is  being 
used  for  both  instructional  purposes  in  the  composite  materials  curriculum 
and  inspection  of  research  material. 

Capability  for  static  and  fatigue  testing  of  composites  has  been 
added  during  the  reporting  period.  An  MTS  Systems  Corporation  20,000 
pound  capacity  servohydraulic  testing  machine  and  a  MINC-11  laboratory 
instrument  computer  have  been  intergrated.  Fig.  10.  Data  acquisition  and 
control  programs  have  been  developed  so  far  for  static  tests  and  single 
cycle  (load  and  unload)  tests.  Structuring  of  a  data  base  has  been 
initiated  with  test  data  from  fracture  tests  on  a  short  glass  fiber/ 
polyester  composite  at  different  temperatures  and  strain  rates.  This  initial 
data  base  will  be  used  as  a  guide  for  structuring  and  developing  data  base 
concepts.  Fortran  Software  for  general  purpose  data  analysis,  including 
curve  fitting  and  graphic  plotting  (both  video  and  hard  copy),  has  been 
developed  for  the  system.  The  program  reads  data  files  from  the  data  base 
(these  files  are  written  by  real-time  data  acquisition  and  control  programs), 
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Figure  8  Tabulation  of  creep  and  recovery  data  for  an  epoxy  resin. 

Column  1  is  time  (sec).  Column  2  longitudinal  strain,  and 
Column  3  transverse  strain. 
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Figure  9  Block  diagram  of  ultrasonic  C-scan  nondestructive  inspection  system 
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provides  interactive  curve  fitting  and  data  editing  and  provides  graphics 
output  (video  or  hard  copy)  of  actual  data  and  curve  fit  parameters. 

An  algorithm  and  Fortran  software  were  developed  to  fit  creep 
data  to  a  power  law  for  creep  compliance,  D(t)  =  Dq  +  D^tn,  using  a 
minimum  squared  error  criterion.  Current  effort  involves  development 
of  data  acquisition  and  control  software  for  fatigue  experiments  and 
related  expansion  of  the  laboratory  experimental  data  base  and  analysis 
capabil ity. 

2.2  Delamination  Fracture* 

Recently  completed  work  on  delamination  fracture  of  a  glass/epoxy 
composite  (Scotchply)  and  a  graphite/epoxy  composite  (AS/3502)  is 
described  in  the  following  paper  (cf.  Appendix  A)  and  thesis  (cf.  Section  4 
for  Abstract),  respectively: 

(i)  Devitt,  D.  F.,  Schapery,  R.  A.,  and  Bradley,  W.  L.,  "A  Method 
for  Determining  the  Mode  I  Delamination  Fracture  Toughness  of 
Elastic  and  Viscoelastic  Composite  Materials",  Journal  of 
Composite  Materials,  Vol .  14,  Oct.  1980,  pp.  270-285. 

(ii)  Hulsey,  R.  C.,  "Delamination  Fracture  Toughness  of  a  Uni¬ 
directional  Graphite/Epoxy  Composite",  M.S.  Thesis,  Texas 
A&M  University,  Dec.  1980. 

Delamination  in  the  opening  mode  of  fracture  was  investigated  using 
a  split  laminate  loaded  as  a  double  cantilever  beam,  which  is  shown  in 
the  inset  in  Fig.  11;  the  fibers  are  parallel  to  the  beam  axis.  Here 
we  shall  discuss  primarily  the  graphite/epoxy  study,  as  the  investigation 
on  Scotchply  is  detailed  in  Appendix  A  (Ref.  (i)). 

Initial  experiments  on  the  graphite/epoxy  composite  were  conducted 
at  ambient  temperature  and  humidity,  and  consisted  of  measuring  load,  P, 
versus  displacement,  2a,  for  various  crack  lengths,  a.  The  beam  analysis, 
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which  includes  geometric  nonlinearity  due  to  large  rotations,  was 

verified  using  two  of  the  measured  quantities  to  predict  the  third, 

which  was  also  measured.  Subsequent  experiments  were  then  run  at 

higher  temperatures  in  a  controlled  atmosphere  where  measurement  of 

instantaneous  crack  length  was  not  convenient;  the  analysis  was  used 

with  P  -  a  data  to  predict  crack  length  as  a  function  of  time.  The 

results  of  these  experiments  indicated  a  critical  energy  release  rate 

2 

of  approximately  1  in-lb/in  except  at  the  highest  temperature  and 
humidity  (200°F  and  95%RH)  (cf.  Fiq.  12)  where  a  20"  increase  in  the  critical 
energy  release  rate  was  observed.  This  increased  toughness  at  higher 
temperature  and  relative  humidity  could  be  interpreted  as  resulting 
from  a  softening  of  the  matrix  as  the  glass-transition  temperature, 

Tg,  approached  the  test  temperature,  T;  thus,  with  much  more  deformation 
taking  place  in  the  vicinity  of  the  crack  tip,  more  energy  would  be 
dissipated  as  the  crack  propagated.  A  significant  variation  in  energy 
release  rate,  G,  with  crack  growth  rate,  3,  was  not  observed.  However, 
a  somewhat  serrated  load-time  record  as  well  as  fractographic  analysis 
of  the  broken  specimens  in  the  scanning  electron  microscope  suggested 
strongly  the  possibility  of  discontinuous  crack  extension  which  might 
obscure  the  true  G  versus  3  relationship. 

Although  more  energy  may  be  dissipated  in  the  material  surrounding 
the  crack  tip  with  an  increase  in  temperature  and  moisture  level,  as 
noted  above,  the  more  severe  environment  may  in  some  cases  reduce  the 
total  amount  of  external  work  required  for  the  local  resin  fracture 
process;  e.g.  crack  propagation  in  a  thermorhelogically  simple  elastomer 
is  facilitated  by  an  increase  in  temperature  [5].  Consequently,  there 
could  be  competing  processes,  and  without  a  good  understanding  of  all 
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Figure  12.  Energy  release  rate  versus  delamination  crack  speed  for  graphite/ 
epoxy  (AS/3502)  in  various  environments  (saturated  state). 
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significant  phenomena  related  to  crack  growth,  one  cannot  predict 
behavior  with  any  real  confidence.  In  order  to  address  this 
question,  a  preliminary  theoretical  viscoelastic  analysis  of  delamination 
crack  growth  was  made  using  the  approximate  method  in  [6],  and  the 
results  are  shown  in  Fig.  13.  (The  model  is  that  of  a  thin  cracked 
resin  layer  bonded  between  orthotropic  beams.)  The  energy  release  rate  G 
is  plotted  against  the  resulting  crack  speed  3  using  dimensionless  ratios. 
The  intrinsic  fracture  energy  required  for  material  separation,  r, 
was  assumed  constant,  and  the  resin  was  assumed  to  be  thermorheolocal ly 
simple;  therefore,  all  temperature-dependence  of  the  value  of  G  required 
to  produce  any  given  speed  3  is  due  to  the  effect  of  the  familiar 
time-temperature  shift  factor  on  the  resin’s  creep  compliance  [5]. 
Recognizing  that  a  major  effect  of  moisture  in  epoxy  (and  many  other 
polymers)  on  mechanical  response  is  analogous  to  that  of  temperature, 
we  have  indicated  in  Fig.  13  the  predicted  effect  of  an  increase  in 
both  temperature  and  moisture. 

The  existence  of  a  maximum  in  the  G  -  3  curve  and  the  subsequent 
negative  slope  shown  in  Fig.  13  for  graphite/epoxy  is  due  entirely  to 
the  mechanical  interaction  of  resin  and  fibers.  In  order  to  illustrate 
the  effect  of  fiber  modulus,  we  predicted  the  G  -  3  curve  for  glass/ 
epoxy,  in  which  the  same  resin  and  fracture  properties  were  employed; 
the  G  -  3  slope  turned  out  to  be  nonnegative  (cf.  Fig.  13).  However, 
other  reasonable  choices  for  resin  properties  indicate  that  the  slope 
for  glass/epoxy  can  be  negative  in  some  speed  ranges. 

The  analysis  was  based  on  continuous  crack  growth.  However, 
when  the  G  -  3  curve  has  a  negative  slope,  one  can  show  the  actual 
behavior  is  unstable.  As  the  actual  overall  specimen  geometry  and 


Theoretical  prediction  of  normalized  energy  release  rate  versus  normalized  delaminat 
crack  speed  based  on  representative  fiber/matrix  properties;  the  same  matrix  propert 
used  for  both  composites. 
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loading  condition  (constant  crosshead  rate)  employed  in  our  tests  is 
such  that  average  crack  growth  is  stable  (cf.  Ref.  (i)),  one  can 
expect  crack  growth  to  occur  in  very  small  steps  of  high  speed  prop¬ 
agation  and  arrest.  Further,  at  any  given  speed,  the  energy  release 
rate  is  predicted  theoretically  to  increase  with  an  increase  in 
moisture  and  temperature  whenever  the  intrinsic  G  -  a  curve  has  a 
negative  slope,  as  indicated  in  Fig.  13. 

These  preliminary  theoretical  findings  are  fully  consistent 
with  the  experimental  results  for  graphite/epoxy  discussed  previously. 
Figure  11  compares  graphite/epoxy  and  glass/  epoxy  results;  the  data 
shown  for  the  former  composite  are  from  one  test  sample,  but  are  typical . 

The  noticeable  decrease  in  G  with  crack  speed  for  graphite/epoxy  occurs 
over  a  very  small  speed  range  compared  to  that  in  Fig.  13.  This  latter 
difference  may  be  due  to  the  actual  step-growth  phenomenon,  in  which 
each  small  step  consists  of  very  slow  to  very  high  speed  propagation, 
whereas  the  theory  assumed  continuous  growth. 

Finally,  it  should  be  added  that  epoxy  resins  (without  reinforcement) 
in  some  cases  exhibit  unstable,  step-growth  [7],  Such  behavior  may  be 
predicted  from  a  viscoelastic  fracture  model  in  which  the  material  in 
the  neighborhood  of  the  crack  tip  is  very  soft  compared  to  the  surrounding 
continuum.  We  should  emphasize  that  the  theoretical  models  discussed  here 
are  only  tentative  representations  of  the  actual  material,  and  much  more 
theoretical  and  experimental  work  is  needed  in  order  to  achieve  an 
acceptable  level  of  understanding  of  time-dependent  fracture  in  resins 
and  composites. 

Further  work  to  clarify  the  delamination  fracture  process  in  the 
graphite/epoxy  composite  is  currently  underway;  a  large  portion  of 


this  research  activity  is  leading  to  Master's  theses  for  some  of  the 
current  M.S.  students.  Included  in  this  effort  is  a  study  of  delamination 
under  axial  compression,  as  a  follow-up  to  work  reported  in  the  thesis 
by  C.  D.  Poniktera  (cf.  Section  4).  Also,  mixed  mode  fracture  toughness 
testing  has  been  initiated;  the  split-beam  geometry  in  Fig.  11  is  used, 
but  the  left  end  is  held  in  place  and  different  loads  are  applied  to 
the  top  and  bottom  beams.  More  extensive  fractographical  analysis  of 
specimens  broken  at  TAMU  and  at  General  Dynamics  is  planned  for  this 
year  also.  Furthermore,  a  special  stage  has  been  ordered  for  the 
scanning  electron  microscope,  which  will  allow  composite  specimens 
to  be  broken  in  the  microscope  while  under  observation;  this  should 
greatly  assist  the  modelling  efforts  currently  underway. 

2.3  Effect  of  Physical  Aging  on  Creep  and  Recovery  of  Resin* 

The  recent  work  on  physical  aging  is  discussed  in  the  M.S.  thesis 
by  D.  H.  Metz  (cf.  Section  4);  it  is  a  continuation  of  an  investigation 
started  in  1979  [8].  The  changes  in  mechanical  properties  of  3502 
epoxy  resin+,  as  measured  by  creep/recovery  tests,  were  determined  as 
a  function  of  time  after  quench  from  above  the  short-time  glass  transition 
temperature.  Quench  rate  and  aging  temperature  were  also  varied.  Aging 
occured  most  rapidly  for  a  high  quench  rate,  a  higher  aging  temperature 
and  short  times  after  quenching.  The  physical  aging  was  observed  to 
cause  a  10  -  20%  change  in  the  compliance  of  the  epoxy  (with  the  material 
getting  stiffer  with  aging  time)  at  the  fastest  cooling  rate  (100°F/min) 
and  the  highest  aging  temperature  (300°F).  At  cooling  rates  more  typical 

♦Prepared  by  Dr.  W.  L.  Bradley 
tProvided  by  Hercules,  Inc.,  Magna,  Utah. 


of  some  commercial  processing  cycles  (5°F/min)  and  at  ambient  aging  temperature 
such  as  one  might  also  expect  commercially,  a  neglibly  small  amount  of 
physical  aging  was  observed.  The  epoxy  material  was  thought  to  physically 
age  much  less  than  thermoplastic  materials  and  many  other  thermosetting 
materials  because  the  very  high  degree  of  chemical  crosslinking  constrains 
the  material  on  a  molecular  level  such  that  less  free  volume  expansion 
is  observed  above  T^;  thus,  less  excess  free  volume  is  trapped  on  sub¬ 
sequent  cooling.  Since  physical  aging  is  associated  with  time-dependent 
free  volume  collapse,  this  effectively  reduces  the  effect. 

2.4  Molecular  Structure  -  Property  Studies  of  Resin* 

The  resin  in  glass  and  graphite  fiber  composites  is  the  source  of 
most  of  the  time-dependent  and  ambient  sensitive  properties  and  so 
studies  of  the  resin  properties  are  very  basic  to  an  understanding  of 
composite  properties.  An  important  interface  between  the  studies  on 
molecular  structure  and  those  on  the  mechanical  properties  of  glassy 
resins  is  in  the  study  of  the  glass  transition  temperature  and  related 
effects.  Here  one  hopes  to  be  able  to  relate  the  phenomena  simultaneously 
to  both  the  molecular  structure  and  to  the  mechanical  properties.  We 
have  been  studying  these  effects  on  3502  epoxy  resin  with  the  differential 
scanning  calorimeter  (DSC)  in  an  effort  to  improve  our  understanding  of 
the  effects  of  curing  conditions,  aging,  moisture  and  molecular  structure 
upon  the  resin  properties. 

Most  of  our  effort  has  been  spent  in  an  attempt  to  improve  the 
reproducibility  and  quality  of  the  DSC  runs  since  the  intended  studies 
seem  to  push  the  instrument  to  the  limit  of  its  capabilities. 


♦Prepared  by  Dr.  J.  S.  Ham 


28 


It  is  often  observed  that  the  first  trace  on  a  DSC  sample  is 
ragged  and  so  cannot  be  used,  while  subsequent  traces  are  of  high  quality. 
For  many  purposes,  the  later  scans  are  adequate,  but  to  study  aging 
phenomena  the  first  trace  must  be  used.  This  raggedness  seems  to  be  due 
to  moisture  evaporating  from  the  sample  so  we  are  now  trying  to  keep 
our  samples  extremely  dry  with  some  success. 

One  of  the  problems  considered  is  how  to  measure  the  glass  transition 
temperature  of  the  proprietary  resins  (e.g.,  3502  and  others  used  in  ad¬ 
vanced  composites)  which  is  above  the  decomposition  temperature.  One 
concept  is  to  lower  by  dissolving  moisture  into  the  sample  and  then 
measure  the  lowered  T  .  Unfortunately,  escaping  moisture  obscures  the 
trace,  so  that  we  cannot  maintain  the  desired  level  of  moisture  while 
measuring  T  .  If  the  moisture  level  is  uncertain,  then  the  value  of 
Tg  is  uncertain  so  that  the  changes  we  seek  are  lost. 

The  experimental  procedure  that  offers  the  most  practical  method  to 
determine  Tg  in  these  samples  is  to  scan  rapidly  to  a  high  temperature 
and  then  immediately  cool  back  down  before  much  decomposition  has  taken 
place.  If  one  can  do  this  with  different  scan  rates,  both  up  and  down, 
one  should  be  able  to  correct  for  a  small  rate  of  thermal  decomposition. 
This  approach  requires  some  data  analysis  which  accounts  for  the  time 
constant  of  the  response  of  the  equipment. 

In  principle,  there  are  a  variety  of  response  time  constants  since 
the  heater,  thermometer  and  cell  are  geometrically  separated  and  each 
has  its  own  heat  capacities  and  thermal  conductivies  to  other  parts  and 
the  surroundings.  Since  all  of  the  decay  times  are  exponential,  we  are 
concerned  primarily  with  the  longest  one.  We  have  measured  the  response 
rate  constants  under  various  conditions  in  an  effort  to  develop  a  method 
to  unfold  the  data  so  that  much  of  the  data  is  saved  while  the  equipment 
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is  approaching  a  steady  state.  This  work  is  continuing. 

In  related  work,  I  spent  the  summer  (1980)  at  the  Materials 
Laboratory,  Air  Force  Wright  Aeronautical  Laboratories, investigating 
problems  associated  with  composites.  The  Materials  Laboratory  has 
underway  a  program  to  develop  acetylene  terminated  polymers  for  use  as 
high  temperature  adhesives  and  composite  matrices.  The  thermal  stabil¬ 
ity  of  these  materials  is  excellent  and  the  chemistry  is  very  flexible 
so  that  many  different  polymers  may  be  made.  Unfortunately,  there  are 
few  guidelines  as  to  which  polymer  will  best  satisfy  the  eventual  end 
uses.  Since  the  cost  is  very  high  to  develop  the  relevant  chemistry, 
there  is  a  danger  that  a  premature  choice  will  be  made, resulting  in  a 
material  far  from  optimum.  Therefore,  they  have  developed  a  program  to 
collect  a  variety  of  data  on  each  polymer  or  version  of  the  polymer, 
utilizing  extremely  small  amounts  of  material.  One  of  the  many  items 
measured  is  the  fracture  toughness.  While  this  information  cannot  yet 
be  tied  with  precision  to  the  properties  of  the  adhesive  or  composite 
matrix,  it  is  clear  that  there  is  a  relationship  and  that  one  seeks  a 
material  with  high  fracture  toughness. 

My  task  was  to  seek  methods  to  enhance  the  fracture  toughness  of 
resins  so  that  appropriate  monomers  will  be  produced  to  be  studied  in 
detail.  My  conclusions  are  that  while  a  heterogeneous  material  may  be 
the  eventual  solu*'on,  the  development  of  such  a  material  is  difficult 
and  lacks  clear  guideposts  along  the  way.  On  the  other  hand,  it  should 
be  possible  to  loosen  the  crossl inked  network  of  the  polymer  by  intro¬ 
ducing  some  dangling  ends  so  that  the  fracture  toughness  is  improved. 
This  will  reduce  certain  properties,  such  as  the  glass  transition 
temperature,  but  this  temperature  is  already  so  high  that  some  tradeoff 
probably  can  be  made.  (The  length  of  these  loose  ends  was  estimated  by 
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the  entanglement  length  determined  from  viscous  flow  in  the  corresponding 
thermoplastic  material).  The  materials  being  developed  have  an  extremely 
short  entanglement  length  so  that  only  short  dangling  chains  should 
suffice.  To  enhance  fracture  toughness,  this  work  will  be  pursued  at  AFWAL 
this  summer,  where  it  is  hoped  that  my  concepts  can  be  tested.  Work 
conducted  during  the  summer  1980  will  soon  be  published  as  a  Materials 
Laboratory  report. 
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3.  ANALYSIS  OF  COMPOSITES 


3.1  Effects  of  Environment  on  Response  of  Composites* 

The  investigations  include  studies  of  the  deformation  of  non-symmetric 
laminated  plates  due  to  thermal  cool-down,  the  diffusion  of  moisture  into 
laminates  under  fluctuating  ambient  envi ronments ,  and  optimization  of 
cool-down  time-temperature  paths  in  cross-ply  laminates.  The  latter  study 
accounts  for  the  temperature  dependence  of  the  coefficients  of  thermal 
expansion. 

The  following  papers  and  theses  on  this  work  were  prepared: 

(i)  Douglass,  D.A.  and  Weitsman,  Y.,  "Stresses  Due  to  Environmental 
Conditioning  of  Cross-Ply  Graphite/Epoxy  Laminates."  In  Advances 
in  Composite  Materials  (Proceedings  of  3rd  International  Con¬ 
ference  on  Composite  Materials,  Paris,  France).  A.  R.  Bunsell 

et  al . ,  Editors,  Vol .  1,  pp.  529-542,  Pergamon  Press  (1980). 

(ii)  Harper,  B.D.  and  Weitsman,  Y.,  "Residual  Thermal  Stresses  in 
an  Unsymmetrical  Graphite/Epoxy  Laminate".  To  appear  in  the 
proceedings  of  the  22nd  SSD  &  Materials  Conference,  Atlanta, 

GA,  April  6-8,  1981. 

(iii)  Weitsman,  Y.,  "A  Rapidly  Convergent  Scheme  to  Compute  Moisture 
Profiles  in  Composite  Materials  Under  Fluctuating  Ambient 
Conditions",  Texas  A&M  University  Report  No.  MM3724-81-6. 

(iv)  Harper,  B.D.,  "Residual  Thermal  Stresses  in  an  Unsymmetrical 
Cross-Ply  Graphite/Epoxy  Laminate",  M.S.  Thesis,  Aug.  1980. 

(v)  Lott,  R.S.,  "Moisture  and  Temperature  Effects  on  Curvature 

of  Anti -Symmetric  Cross-Ply  Graphite/Epoxy  Laminates",  M.S.  Thesis, 
Dec.  1980. 


An  additional  investigation  on  "Optimal  Cooling  of  Cross-Ply  Laminates 
and  Adhesive  Joints"  is  presently  in  preparation  by  Y.  Weitsman  and  B.D. 
Harper. 
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The  essential  finding  of  Ref.  (i)  are  illustrated  in  Figs.  3-5 
therein  (Appendix  A).  If  we  consider  in  particular  the  state  of  the 
residual  stresses  upon  termination  of  moisture  and  temperature  condi¬ 
tioning,  we  note  that  the  stress  profile  depends  most  significantly  on 
the  conditioning  history.  When  we  observe  the  solid  and  dashed  curves 
at  the  times  when  our  laminate  was  "Returned  to  Room  Temperature",  we 
note  that  the  stresses  are  of  entirely  reversed  sign,  in  spite  of  the 
fact  that  the  current  levels  of  moisture  and  temperature  are  indentical. 
The  significance  of  the  time-dependent  material  response  is  illustrated 
in  Fig.  5,  which  demonstrates  that  elastic  predictions  may  be  different 
in  both  sign  and  magnitude. 

In  Ref.  (ii),  we  managed  to  demonstrate  that  the  time-dependent 
behavior  is  amenable  to  experimental  detection.  In  addition,  we  also 
observe  that  the  residual  stresses  can  reach  critical  values,  suffi¬ 
cient  to  produce  large  cracks  in  composite  laminates.  Additional 
analytical  and  experimental  work  is  currently  in  progress,  extending 
the  findings  of  Ref.  (ii). 

Reference  (iii)  exhibits  the  coupling  between  the  moisture 
diffusion  process  and  temperature  in  composite  materials.  Figure  1 
therein  (Appendix  A)  demonstrates  that  when  exposed  to  identical, 
fluctuating,  ambient  relative  humidities,  entirely  different  internal 
moisture  profiles  can  develop  within  laminates  under  different  temper¬ 
ature  histories. 


3.2  Analysis  by  the  Finite  Element  Method:  Constitutive  Models  and 

Fracture  Mechanics* 

The  nonlinear  viscoelastic  constitutive  model  developed  during 
previous  years  hasbeen  incorporated  into  the  finite  element  program, 
AGGIE.  The  model  has  been  implemented  and  tested  for  the  plane  stress, 

plain  strain,  and  axi-symmetric  cases.  The  3/D  model  has  been  written 
but  has  not  been  tested  yet.  A  number  of  isothermal  and  nonisothermal 
creep  and  creep-plasticity  problems  have  been  investigated  and  results 
compared  to  available  experiments.  A  report  on  this  work  should  be 
completed  by  June,  1981. 

Related  experimental  work  on  the  nonlinear  viscoelastic  character¬ 
ization  of  composites  was  initiated;  it  is  described  in  the  thesis 
(Section  4) : 

Kerstetter,  M.S.,  "Nonlinear  Viscoelastic  Characterization  of 

AS/3502  Graphite/Epoxy  Composite  Material,"  M.S.  thesis,  Texas 

A&M  University,  Dec.  1980. 

The  finite  element  program  has  been  modified  to  calculate  nodal 
forces  in  the  vicinity  of  crack  tips  and  the  resulting  intensity  factors 
and/or  J-integral  values.  Three-dimensional  calculations  of  Mode  I,  II, 
and  III  cracking  are  thus  possible  with  the  program.  Although  not  yet 
tested,  elastic-plastic  calculations  for  the  J  integral  may  be  possible 
with  the  program. 

3.3  Damage  Growth  and  Fracture  of  Composites* 

Recently  completed  work  on  theoretical  models  for  damage  growth 
and  fracture  is  described  in  the  following  two  papers  (cf.  Appendix  A): 
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(i)  Schapery,  R.A.,  "On  Constitutive  Equations  for  Viscoelastic 
Composite  Materials  with  Damage,"  Proc.  NSF  Workshop  on 
Damage,  Cincinnati,  April,  1980. 

(ii)  Schapery,  R.A.,  "Nonlinear  Fracture  Analysis  of  Viscoelastic 
Composite  Materials  based  on  a  Generalized  J  Integral  Theory", 
Proc.  1st  Japan--U.S.  Conference  on  Composite  Materials, 
Tokyo,  Jan. ,  1981 . 

In  Reference  (i),  the  material  is  assumed  to  be  linearly  visco¬ 
elastic  for  a  fixed  state  of  damage.  A  rigorous  derivation  is  then 
given  to  obtain  the  effect  of  microcracking  and  rejoining  of  crack 
surfaces  (with  or  without  healing)  on  the  macroscopic  constitutive 
equation.  Residual  stress  effects  due  to  temperature  and/or  moisture 
are  included.  The  results  have  served  to  provide  guidelines  for  the 
development  of  more  general  nonlinear  viscoelastic  constitutive  equations 
with  damage  and  a  new  theory  for  the  prediction  of  initiation  of  crack 
growth  and  crack  speed  in  such  composites  [9].  Some  results  and  appli¬ 
cations  of  this  fracture  theory  are  discussed  in  Ref.  (ii). 

Experimental  work  on  microcrack  growth  in  a  viscoelastic  composite 
(Scotchply)  is  described  in  the  thesis  (cf.  Section  4), 

Lehman,  M.W.,  "An  Investigation  of  Intra-Ply  Microcrack  Density 

Development  in  a  Cross-Ply  Laminate",  M.S.  Thesis,  Texas  A&M 

University,  Dec.  1980. 

The  initial  microcrack  state  and  subsequent  growth  due  to  load 
application  are  described;  a  florescent  dye  penetrant  is  employed  to 
help  identify  microcrack  geometry.  A  comparison  is  made  between  theo¬ 
retically  and  experimentally  determined  stiffness  changes  due  to 
microcracking.  This  activity  represents  a  portion  of  the  work  underway 

to  develoo  analytical  models  of  damage  growth  in  viscoelastic  composites 
through  a  combination  of  experimental  and  theoretical  approaches. 


4.  GRADUATE  RESEARCH  ASSISTANT  ACTIVITIES 


4.1  Suirmary 

The  second  group  of  graduate  engineering  students  to  participate 
in  the  AFOSR  research  project  entered  the  program  in  September,  1979 
and  graduated  with  a  Master  of  Science  degree  by  December,  1980.  Results 
of  their  research  are  reported  in  the  following  theses: 

1.  Harper,  8.D.,  "Residual  Thermal  Stresses  in  an  Unsymmetrical 
Cross-Ply  Graph: te/Epoxy  Laminate." 

2.  Hulsey,  R.C.,  "Delamination  Fracture  Toughness  of  A  Uni¬ 
directional  Graphite/Epoxy  Composite." 

3.  Kerstetter,  M.C.,  “Nonlinear  Viscoelastic  Characterization 
of  AS-3502  Graphite/Epoxy  Composite  Material." 

4.  Lehman,  M.W.,  "An  Investigation  of  Intra-Ply  Microcrack 
Density  Development  in  a  Crossply  Laminate." 

5.  Lott,  R.S.,  "Moisture  and  Temperature  Effects  on  Curvature 
of  Anti -Symmetric  Cross-Ply  Graphite/Epoxy  Laminates." 

6.  Metz,  D.H.,  "Experimental  Investigation  of  Free  Volume 
Concepts  in  Relationship  to  Mechanical  Behavior  of  an 
Epoxy  System  Subjected  to  Various  Aging  Histories." 

7.  Poniktera,  C.D.,  "Application  of  Energy  Release  Rate 
Principles  to  Compression  Debonding." 


Abstracts  are  given  in  Section  4.2.  Copies  of  the  theses  will  be 
provided  upon  written  request  to  the  Principal  Investigator  (R.A.  Schapery). 

The  current  group  (starting  September,  1980)  consists  of  seven  (7) 

M.S.  students  and  one  (1)  Ph.D.  student.  The  topics  are  listed  below; 
most  studies  involve  both  experimental  and  theoretical  work. 
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M.S.  Theses-- 

1.  Shear  deformation  effects  in  highly  anisotropic  laminates 
(Coulter/Wei tsman) . 

2.  Micromechanisms  of  delamination  fracture  (Williams/Bradley). 

3.  Mixed-mode  delamination  fracture  (VanderKley/Bradley) . 

4.  Compression-induced  delamination  (Earley/Jerina ) . 

5.  Delamination  under  complex  loading  histories  (Cull en/Jerina ) . 

6.  Delamination  fracture  analysis  including  effect  of  matrix 
damage  (Arenburg./Schapery) . 

7.  The  effect  of  elliptical  hole  shape  on  the  design  of  pin 
loaded  filament  wound  fiberglass  tension  lugs  (Braswell/ 
Alexander) . 

Ph.D.  Dissertation-- 

8.  Envi ronmental  effects  in  unbalanced  laminates  (Harper/Wei tsman ) . 

The  next  quarterly  report  will  contain  thesis  proposals  describing 
the  research  planned  under  each  of  these  topics.  Besides  conducting  this 
research,  the  students  are  involved  in  academic  courses,  as  described 
in  the  brochure  in  Appendix  8. 

In  addition  to  graduate  students,  a  few  select  undergraduate 
students  assist  in  the  various  research  activities.  The  participation 
of  these  undergraduate  students  aids  in  the  research  and  helps  to  acquaint 
them  with  composite  materials  prior  to  enrolling  in  the  graduate  program. 

4 . 2  Abstracts  of  M.S.  Theses  Comp 1 e ted  i n 1 980 
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ABSTRACT 

Residual  Thermal  Stresses  in  an  Unsymmetrical  Cross-Ply 
Graphite/Epoxy  Laminate.  (August  1980) 

Brian  Douglas  Harper,  B.S.,  Ag.  En.,  Texas  A&M  University 
Chairman  of  Advisory  Committee:  Dr.  Y.  Weitsman 

This  thesis  presents  a  method  for  determining  the  residual 
thermal  stresses  in  AS-3502  graphite/epoxy  laminates  due  to  cool-down 
from  their  cure  temperature.  Also  included  is  a  method  for  de¬ 
termining  the  optimal  time-temperature  path  that  will  minimize  these 
residual  stresses. 

The  analysis  considers  the  time-dependent  behavior  of  the 
material  and  all  calculations  employ  recent  data  on  the  thermo¬ 
viscoelastic  response  of  the  AS-3502  graphite/epoxy  system. 

The  viscoelastic  analysis  is  verified  through  curvature  measure¬ 
ments  of  unsymmetric  cross-ply  plates  fabricated  from  the  AS-3502 
graphite/epoxy  material. 
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ABSTRACT 

Delamination  Fracture  Toughness  of  A  Unidirectional 
Graphite/Epoxy  Composite.  (December  1980) 

Roy  Charles  Hulsey,  B.S.,  Texas  A&M  University 
Chairman  of  Advisory  Committee:  Dr.  Walter  L.  Bradley 

The  opening  mode  delamination  fracture  toughness  of  a  graphite/ 

epoxy  composite  under  varied  temperature,  humidity  and  crack  growth 

rates  is  investigated  experimentally.  Energy  release  rate  for  a 

stably  growing  crack  (Gv)  is  determined  using  a  double  cantilever 

beam  specimen  and  a  linear  elastic  fracture  mechanics  analysis  coupled 

with  nonlinear  beam  theory.  The  temperature  range  of  75F  to  200F  and 

0%  RH  to  95%  RH  has  been  studied.  A  significant  effect  of  temperature 

and  humidity  on  is  observed  only  at  200F  and  95%  RH  for  the  system 

2 

studied.  Gy  values  of  0.955  in-lb/in  (standard  deviation  =  0.08)  are 

2 

determined  in  general  with  a  value  of  1.14  in-lb/ in  determined  at 
the  200F  -  95%  RH  condition.  The  energy  release  rate  is  not  found  to 
be  significantly  affected  by  crack  growth  rates  in  the  range  0.01  to 
10.0  in/min. 
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ABSTRACT 

Nonlinear  Viscoelastic  Characterization  of  AS-3502 
Graphite/Epoxy  Composite  Material.  (December  1980) 

Michael  Scott  Kerstetter,  B.S.,  Texas  A&M  University 
Chairman  of  Advisory  Committee:  Dr.  K.L.  Jerina 

The  objective  of  this  paper  is  to  study  the  creep  and  recovery 
response  of  a  composite  subjected  to  several  high  stress  levels  in 
a  high  humidity  environment.  Strain  versus  time  data  obtained  from 
uniaxial  creep  and  recovery  tests  were  used  to  characterize  visco¬ 
elastic  deformation  in  a  [i45]2^  laminate.  A  description  and  examin¬ 
ation  of  the  effectiveness  of  two  data  collection  schemes  are  pre¬ 
sented  along  with  a  discussion  of  some  important  experimental  aspects 

such  as  generation  of  the  test  environment  and  mechanical  and  humidity 
conditioning.  Data  were  obtained  from  creep  and  recovery  for  several 

stress  levels  which  were  reduced  using  graphical  shifting  procedures 

from  which  nonlinear  material  parameters  were  determined.  Conclusions 
and  recommendations  for  future  research  are  presented. 
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ABSTRACT 

An  Investigation  of  Intra-Ply  Microcrack  Density 
Development  in  a  Crossply  Laminate.  (December  1980) 

Michael  William  Lehman,  B.S.,  C.E.,  Texas  A&M  University 
Chairman  of  Advisory  Committee:  Dr.  R.  A.  Schapery 

An  investigative  technique  is  qualified  as  an  experimental  tool 
to  aid  in  the  quest  involving  identification  of  parameters  controlling 
damage  accumulation  in  laminated  composites.  Specifically,  the  flaw 
state,  encompassing  geometry,  size,  and  density  as  a  function  of  a 
monotonic  load  history,  is  characterized  for  a  [DO^/O^  glass  epoxy 
test  laminate.  Founded  upon  direct  and  indirect  experimental  observa¬ 
tions,  a  subjective  description  for  flaw  or  crack  development  on  both 
a  local  and  global  scale  is  also  presented. 

Theoretical  bounds  for  strain  energy  change  are  established  for 
the  investigated  test  laminate.  These  bounds  are  used  to  compare  and 
qualify  the  experimental  data  in  the  range  for  which  the  relationship 
between  residual  modulus  and  transverse  crack  density  is  linear. 
Finally,  possible  candidates  controlling  the  global  rate  of  strain 
energy  change  for  accumulating  damage  are  identified. 


ABSTRACT 


Moisture  and  Temperature  Effects  on  Curvature  of 
Anti -Symmetric  Cross-Ply  Graphite/Epoxy  Laminates.  (December  1980) 
Randall  Stephen  Lott,  BAE,  Georgia  Institute  of  Technology 
Chairman  of  Advisory  Committee:  Dr.  Y.  Weitsman 

This  thesis  presents  a  method  for  analytically  determining  the 
mid-plane  strains  and  curvatures  of  anti -symmetric  cross-ply  graphite/ 
epoxy  laminated  plates  which  are  exposed  to  high  humidities  at 
elevated  temperatures. 

The  analysis  considers  temperature  dependent  moisture  diffusion 
and  time/temperature/moisture  dependent  stress  relaxation.  Recent 
data  on  the  hygrothermal -viscoelastic  behavior  fo  the  AS/3502 
graphite/epoxy  system  is  employed  in  the  calculations. 

Results  of  both  elastic  and  viscoelastic  analyses  are  presented 
and  compared  to  measured  curvatures  of  anti -symmetric  cross-ply 
plates  fabricated  from  the  AS/3502  system  and  exposed  to  high  temp- 
erature/humidi ty  envi ronments . 
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ABSTRACT 

Experimental  Investigation  of  Free  Volume  Concepts  in  Relationship 
to  Mechanical  Behavior  of  an  Epoxy  System  Subjected 
to  Various  Aging  Histories.  (December  1980) 

Daniel  Hugh  Metz,  B.  S.  University  of  Illinois 
Chairman  of  Advisory  Committee:  Dr.  W.  L.  Bradley 

An  epoxy  resin  commonly  used  in  advanced  composite  materials  for 
aerospace  application  was  tested  for  changes  in  viscoelastic  behavior 
after  being  quenched  from  above  T^  to  300°F,  200°F  and  75°F  and  then 
isothermally  aged.  The  qualitative  correlation  between  the  changes 
in  viscoelastic  response  and  free  volume  is  discussed.  In  general 
more  total  physical  aging  and  a  more  rapid  physical  was  observed  at 
higher  temperatures.  Data  obtained  from  these  experiments  is  useful 
in  contributing  to  an  overall  understanding  of  factors  important  to 
the  optimization  of  processing  parameters  in  the  manufacture  of 
composite  material  components. 


A 


1 
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ABSTRACT 

Application  of  Energy  Release  Rate 
Principles  to  Compression  Debonding.  (August  1980) 

Christopher  Dale  Poniktera,  B.A.  Engineering  Science  with  a 
Specialization  in  Applied  Mechanics,  University  of  California 
Chairman  of  Advisory  Committee:  Dr.  K.  L.  Jerina 

The  objective  of  this  research  effort  was  to  provide  a  critical 
assessment  of  the  state-of-the-art  of  the  analysis  of  compressive 
debonding  of  fiber-reinforced  composite  materials  and  to  investigate 
the  applicability  of  using  an  energy  release  rate  approach  as  a  means 
of  analyzing  this  phenomenon. 

The  method  developed  utilizes  linear  beam-column  theory  and 
assumes  prior  knowledge  of  the  critical  strain  energy  release  rate 
of  the  system  analyzed. 

An  experimental  program  employing  a  polymethylmethacrylate  (PMMA) 
model  system  was  conducted  to  verify  the  proposed  analytical  technique. 
The  model  system  consisted  of  a  symmetrical  beam-column  arrangement 
subjected  simultaneously  to  axial  and  lateral  loads. 

Results  of  the  study  indicated  good  agreement  between  experi¬ 
mental  and  theoretical  model  displacement  predictions,  but  predic¬ 
tion  of  debond  propagation,  based  on  know  Mode  I  critical  strain 
energy  release  rate  values  (GjC),  was  not  obtained.  The  author 
offers  several  possible  reasons  for  this  apparent  discrepancy. 
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5.  PROFESSIONAL  PERSONNEL  INFORMATION 


5.1  Faculty  Research  Assignments 

Each  participating  faculty  member  is  responsible  for  the  research 
conducted  in  at  least  one  specific  area  of  investigation,  as  shown  below. 
In  addition,  most  serve  as  chairmen  of  one  or  more  of  the  graduate 
advisory  committees  for  M.S.  students  and,  as  such,  direct  their 
students'  research  project.  The  faculty  also  contribute  to  other  research 
activities  on  the  project  by  serving  on  student  advisory  committees, 
through  technical  meetings,  informal  discussions,  and,  in  some  cases, 
through  specific  research  work. 

The  Principal  Investigator  (R.  A.  Schapery)  has  responsibility  for 
overall  technical  direction  and  coordination  and  for  project  management. 

In  addition  he  has  direct  responsibility  for  certain  research  work,  as 
noted  below. 


Faculty  Member/Departmental  Affiliation 

Dr.  Walter  Bradley/Mechanical 
Engineering 

Dr.  Walter  Haisler/Aerospace 
Engineering 

Dr.  Joe  Ham/Physics 

Mr.  Bob  Harbert/Civil  Engineering 

Dr.  Ken  Jerina/Civil  Engineering 

Dr.  Richard  Schapery/Aerospace  and 
Civil  Engineering 


Primary  Research  Responsibility 

Physical  Aging  Behavior, 
Delamination  Fracture  Properties. 

Development  of  Finite  Element 
Models. 

Curing  and  Aging  Studies. 

Duomorph  Gage 

Experimental  Data  Base, 

Mechanical  and  Failure  Property 
Characterization. 

Principal  Investigator  and 
Theoretical  Models  for  Physical 
Aging,  Damage  Growth,  and 
Fracture. 


Dr.  Jack  Weitsman/Civil  Engineering 


Constitutive  Relations, 
Environmental  Effects. 


5.2  Additional  Professional  Staff 


Mechanics  and  Materials  Center 

Mr.  Carl  Fredericksen  -  Electronics  Technician 
Mr.  William  Eue  -  Computer  Programmer 


5.3  Spoken  Papers  and  Lectures  at  Conferences  and  Other  Professional 
Activities  of  the  Faculty  Related  to  Composite  Materials 
(1  January  1980  -  31  De  ember  1980): 

W.  L.  Bradley 

Invited  Lectures  and  Conference  Presentation: 

"J-Integral  Fracture  Toughness  Studies  of  Cast  Iron,"  American 
Founderymen  Society,  Annual  Meeting,  St.  Louis,  April  1980. 

"Fracture  Toughness  of  Nodular  Cast  Iron,"  Sandia  Laboratories, 

Sept.  1980. 

"Size-Effect  on  Toughness  Measurements  of  Nodular  Cast  Iron," 

American  Institute  of  Metallurgical  Engineers  (AIME),  Annual 
Meeting,  Las  Vegas,  Feb.  1980. 

Technical  Committee  Membership: 

ASTM  E-24  Committee  on  Fracture,  Metal  Properties  Council,  sponsored 
jointly  by  ASME,  ASM,  and  ASTM. 

Awards : 

Best  Materials  paper  of  the  year  by  the  American  Nuclear  Society: 
"Corrosion  and  Mechanical  Behavior  of  Iron  in  Liquid  Lithium," 

Nuclear  Technology,  Vol.  39,  1978,  pp.  75-83. 

Named  Haliburton  Professor  of  Mechanical  Engineering. 

K.  L.  Jerina 

Invited  Lecture  and  Conference  Presentation: 

"Effective  Moduli  of  Three  Dimensionally  Reinforced  Fibrous 
Materials,"  Gordon  Conference  lecture  on  Composite  Materials, 

Santa  Barbara,  Jan.  1980. 

"Viscoelastic  Characterization  of  a  Random  Fiber  Composite  Material 
Employing  Micromechanics,"  Short  Fiber  Reinforced  Composite  Materials 
ASTM,  Minneapolis,  May  1980. 
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Technical  Committee  Membership: 

ASTM  E-9  on  Fatigue,  Corresponding  Secretary 

ASTM  E-24  on  Fracture 

ASTM  D-30  on  Composite  Materials 

SAE  Fatigue  Design  and  Evaluation  Committee 
Chairman  of  Task  Force  on  Composite  Materials 

R.  A.  Schapery 

Invited  Lectures  and  Conference  Presentations: 

"Fracture  and  Fatigue  of  Viscoelastic  Materials,"  University  of 
Wyoming,  Jan.  1980. 

"On  Constitutive  Equations  for  Viscoelastic  Composite  Materials 
with  Damage,"  NSF  Workshop  on  Damage,  Cincinatti,  May  1980. 

"A  Complex  Modulus  Gage-  The  Duomorph",  Annual  Meeting  of  the 
Acoustical  Society  of  America,  Atlanta,  April  1980. 

"Nonlinear  Approximate  Analysis  of  Solid  Propellant  Grains", 

Structures  and  Mechanical  Behavior  Meeting  of  Rocket  Propulsion 
Group,  Sept.  1980. 

"Composite  Materials  for  Structural  Design",  Sixth  Annual 
Mechanics  of  Composites  Review  Meeting,  Oct.  1980. 

"Application  of  a  Generalized  J  Integral  to  Fracture  of  Linear 
and  Nonlinear  Viscoelastic  Composite  Materials",  17th  Annual 
Meeting,  Society  of  Engineering  Science,  Atlanta,  Dec.  1980. 

Technical  Coimittee  Membership: 

National  Materials  Advisory  Board,  National  Academy  of  Sciences, 
Committee  on  "High  Temperature  Metal  and  Ceramic  Matrix  Composites". 

Awards: 

Named  Alumni  Professor  and  Distinguished  Professor  of  Aerospace 
and  Civil  Engineering. 

Y.  Weitsman 

Conference  Presentation: 

"Stresses  Due  to  Environmental  Conditioning  of  Cross-Ply  Graphite/ 
Epoxy  Laminates",  3rd  International  Conference  on  Composite  Materials, 
Paris,  Aug.  1980.  Also  participated  in  panel  discussion. 
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Technical  Committee  Membership: 

AIAA  Subcommittee  on  Design  Allowables  for  Composite  Materials. 

In  addition  to  the  above  activities,  the  faculty  attended  several 
conferences  on  composites,  published  papers  on  other  projects,  and  worked 
as  consultants  to  industry. 
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APPENDIX  A 


Recently  Completed  Reports  and  Publications 


1.  "Stresses  Due  To  Environmental  Conditioning  of  Cross-Ply 
Graphite/Epoxy  Laminates"  by  D.A.  Douglass  and  V.  Weitsman. 

2.  "Residual  Thermal  Stresses  in  an  Unsymmetrical  Cross-Ply 
Graphite/Epoxy  Laminate"  by  B.D.  Harper  and  Y.  Weitsman. 

3.  "A  Rapidly  Convergent  Scheme  to  Compute  Moisture  Profiles 

in  Composite  Materials  Under  Fluctuating  Ambient  Conaitions" 
by  Y.  Weitsman. 

4.  "On  Constitutive  Equations  for  Viscoelastic  Composite  Materials 
with  Damage"  by  R.A.  Schapery. 

5.  "A  Method  for  Determining  the  Mode  I  Delamination  Fracture  Tough¬ 
ness  of  Elastic  and  Viscoelastic  Composite  Materials"  by  D.F.  Devitt, 
R.A.  Schapery,  and  W.L.  Bradley. 

6.  "Nonlinear  Fracture  Analysis  of  Viscoelastic  Composite  Materials 
Based  on  a  Generalized  J  Integral  Theory"  by  R.A.  Schapery. 


STRESSES  DUE  TO  ENVIRONMENTAL 
CONDITIONING  OF  CROSS-PLY 
GRAPHITE/EPOXY  LAMINATES 


by  David  A.  DOUGLASS 

Sr.  Structures  Engineer,  Bell  HtHcopter  Textron 
Fort  Worth,  Texas 


and  Yechiel  WEITSMAN 

Professor 

Mechanics  and  Materials  Center 
Texas  AAM  Unhrerxity 
College  Station,  Texas 


Thi*  paper  M’mctnH  thi-  Internal  residual  stresses  due  Co  the  environmental 
.ul  It  Ion  Inn  ° f  ct»*Hs-jily  graph  it  e/epoxy  laminates,  where  Moisture  le  Induced  into 
:ns  material  hv  exp  .sum*  t »»  h’gh  relative  humidity  at  elsvatad  temperatures . 

fh-  stress,  fields  rfbuiii-.g  l  too  . .ndir  toning  at  two  temperature  levels  are 
*  vs  1>.  .  ffu  (or  a  bo  1  m.  *d,  swrot-crK  cross-ply  graph  It  f  /  «poxy  laminate  by  Beans  of 
h«>th  elastic  and  v l soocl.ibt  1.  analyses.  Th«*  formulation  considers  temperature- 
dfpenJent  moisture  diffusion  and  time,  temperature  snd  Bolsturo  dependent  stress 
•  fixation.  The  computation*  are  bused  upon  recent  data  snd  employ  realistic  values 
l  juaterinl  param*  (  «m  *  . 

Tin.  analysis  shows  chat  the  vl«r  _,elsst  n  stresses  ere  much  smaller  then  those 
predicted  bv  elastic  analysis,  and  that  conditioning  at  150*F  results  in  residual 
tresses  whwb  are  Kuuller  than  thos*  reached  during  conditioning  at  130*F. 
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IKTROQUCTION 


Increasing  cocicena  with  the  effect*  of  temp*r*turs  and  molature  on  the  par- 
foneenct  of  composites  dwend»  the  develupmant  of  accslsratad  conditioning  tech¬ 
niques.  Theee  conditioning  scheme*  attempt  to  simulate  «xtr«M  expoeuree  to 
molscure  and  temperature  which  nay  be  encountered  In  service  life. 

Since  **>1  elute  diffuses  extremely  alow  rhe  conditioning  process  is  ec celereted 
by  employing  elevated  toaperatureb  which  arcslerete  the  molature  sorption  end 
shot  ten  che  duration  of  the  experiment. 

It  wee  observed  thee  :ertaln  conditioning  schemes  induce  damage  in  the  com¬ 
posite.  In  order  to  develop  time  saving,  yet  non-damaging,  laboratory  txpsrimants 
it  la  necessary  to  predict  the  stresses  which  arise  due  to  molature  and  temperature 
la  conposlts  laminates- 

In  t'.ls  paper  the  moieture  and  temperature  affects  are  Incorporated  into  e 
viscoelastic  constitutive  relation,  which  Is  based  upon  creep  data  et  various 
levels  of  r emperatursa  end  moisture  contents-  It  Is  ooted  that  both  environment el 
factors  act  as  swelling  agents,  which  introduce  stresses  In  the  presence  of  geo¬ 
metric  constraints.  On  the  other  hand,  both  moisture  end  temperature  enhance  che 
relaxation  of  stresses,  thus  compensating  for  their  initial  effects.  While  this 
'’competition"  is  accounted  for  In  the  viscoelastic  representation  it  cannot  be  con¬ 
sidered  within  the  context  of  on  clastic  reaponea.  Consequently ,  test  articles 
which  were  conditioned  In  different  manners  will  contain  disparate  residual 
stresses  and  perform  differently  under  loads. 

In  this  paper  we  consider  a  symmetric,  balanced,  cross-ply  laminate.  Results 
art  presented  for  two  conditioning  schemes  and  comparison  Is  provided  with  an 
elastic  Analysis-  Additions]  result*  and  further  details  are  given  in  Ref.  (1). 


pg£U?i inary  considerations 

.t.'f  .  rarity  ve  prevent  individually  the  major  factors  present  In  molature  coo- 

d  1  r  inning 

Moisture  t>  1  f  f  is  ion 

The  diffusion  tor  t’.ln  •  oncost  tt-  laminate*  1*  found  to  follow  Fick/s  second 
Itw  ,2).  f  3 J  *  whlih  In  one  dimenaion  read* 


*•  _  .♦ ?  a 

■7  *  u  »,7 


(1) 


The  solution  :o  (1).  und»  f  «:«*ndy  ambit  nt  conditions  Is  given  In  tvo  alternate 
forms  \U) 


«•' t .1 )  *  »,  -  («*  -  V 


"  aa  “  <®-\  “  °»o }  C-f  l 


*  er.c 


exp(-Unt») 


(2) 


(J) 


where  the  d  ico**n>  ionl  t  rifte  t*  is  given  by  t*  •  — ,  ,  and  in  (3) 

h* 

u  .  7n\UI 


!t  ran  be  shirwn  th.»t  tor 
form  la  more  idvanrakeoo 
scan:  ambient  m-lsture,  is 


c*  •  0 . (2)  converge*  rapidly  while 
in  Ij.  (2'  and  (3)  a  1*  molature, 
•'  uniform  initial  molature,  t  is  th* 


for  t*  >0.29 
m*  is  che  con- 
•oat lei 


otJuut*  dvi  u.b  i»u-  thickness,  t  1*  time  and  D  lit*  mo  incur*  diffueivity. 


1C  ha«  been  found  {s|  that  the  saturation  level  m*  depends  on  the  relative 
v  .Tidily  RH.  In  gr aphl ta/epory  eyattiua  *  linaar  relationship  a*  *  l(RH)  la  employ 
,j .  In  addition,  tha  diffusion  co«f fie  lent  D  la  temperature  dapandanl  [6).  Thla 
.pondaoce  1*  given  by  D(T)  •  Aj exp(-Bj /T) . 

Swelling  Due  to  Moisture 

The  Irani verse  •welling  of  noat  gt uphlte/epory  ay at ami  due  to  molatura  abaorp 
,.*n  la  akacchad  In  Mg.  I  (7,  8].  Tha  phenomenon  la  Chirac  tar  1  cad  by  a  moisture 
npenalon  coafficlent  Sf. 


F 1  .  1.  Transverse  Swelling  Vs.  Molatura  Content 


The  transverse  strain  due  to  txjistutf  swelling  la  expraaaed  aa  follows 


t-j-  •  dr  (o  -  tn.  )  for  m  j_  m^ 
*7  •  0  for  »  <  « > 


<*> 


In  (4)  »i  In  n  threshold  value  below  which  molatura  la  not  accompanied  by  any 
•.ransverae  availing. 

Meaauf«*n<  .it-  l  i  '  Indicate  that  the  longltud  iral  swelling  coefficient  la 
eligibly  small  and  ««  shall  take  Bl  * 

I emperor  ore  ultfusion  and  Thermal  rxpanalon 

Oats  :*n  teapcr.it urc  diffusion  ;nd.<  «t«*  tfi.it,  for  typical  laminates,  thle  pro- 
na  precr<Ms  ouch  raster  than  all  other  t  toe -dependent  proceaaes.  Thera  fora  we 
.‘mII  disregard  temperature  diffusion  and  consider  specially  unlfona,  thermally 
'  *  .  1 1  Ibratcd  states. 

In  ur  n'  i  r  »-c  t  lonai  l.rinjc  the  ihrrc.il  strains  ere  coapletaly  charac  terltad  by 
.*  lo.iglnCudlnal  and  transverse  a'fK  h  irnti  of  thermal  expansion  end  o^* 

Geometry  and  N. ■  t . it  lo n 

We  ah.il!  ion«iid.*r  syunetiK.  -J* '**».»*'  lay-ops  and  let  the  principal  material 
v.  i  co  me  id*  with  the  x-y  coordinate* .  The  x  axis  la  choaen  parallel  to  tha 
ibeta  In  th**  v  plip*.  Whenever  necessary  Je  shall  employ  subscripts  L  and  T  to 
not#  longitudinal  and  transver-**  directions.  See  Fig.  2. 


‘I 


ViicogUitlc  Behavior 

The  vltcotUitlc  response  of  graphite /epoxy  laminates  lo  the  traomna 
direction  la  of  paramount  significance  during  envlronawi ntal  conditioning.  Thm 
t 1—  dependant  traoavaraa  cosgil  lance  Sy  can  be  expraaeed  by  a  "power  lav" 
aquation  (9) 

%  h))  “  (J) 

•here  Dq  la  the  Initial  compliance,  t  la  time.  a(T,H)  la  tha  shift  factor  function 
which  Jepanda  on  humidity  arul  temperature,  and  D j  e».d  q  ara  material  cooatanta. 

The  dependanca  of  Dq  on  T  and  H  can  be  expressed  by  (9|  Oq(T»H)  •  aTl  ♦  VC  ♦  cH  4*  d. 

However,  In  order  to  simplify  computet Ions  we  ahall  approximate  Dq(T.H)  by  ite 
average  value  over  the  range  of  T  and  H  In  each  conditioning  scheme.  Thie  approxi¬ 
mation  entails  error#  of  abour  51. 

The  data  oo  tha  ehlfc-factor  function  e<T,H)  can  ba  expressed  by 

a (T , H)  -  altT)a?(H)  where 

log  s,(T)  •  aj  ♦  s„(l/T)  a5(l/T)?  and  <b) 

•j  )  -  a6mbi 


Tha  compliance*  Sj  and  Sl?  do  not  exhibit  time,  temperature  and  moleture  dspandencm. 
Elastic  Compliance!  anu  Kudu 1 1 

U«  ahall  take  all  elastic  compliances  and  taodull  as  tha  Initial  viscoelatic 
values.  Moduli  are  Inverses  of  compliances,  ss  follows 


-  -S,?/£i  ,  where  A  -  S^Dq  -  Sj|. 


<» 


Dq/6  ,  Ct  •  s^/a  .  C i j 


gO^blned  StrMg  -  Strata  jUUtioni 
tittle 

Tha  coat)  load  elastic  itrui-»C(ala  relatione  -for  the  pi  la*  considered  Kanin, 
«nd  far  a  stats  of  plana  atraaa,  ara 

9i?  *  Cl(k  ♦  Cj2«y  -  ^CL°l  *  Ci2«x)AT  “  C12®T^b  “  *|  )  *"  ^L®l  • 

“  C-jc  t  *  Cuty  -  (Cja-  ♦  C|2®»)AT  -  Cfl^fe  -  *1 )  •  Cu®*- 

<•> 

o,°  -  CW«K  ♦  Cfty  -  (C*xoT  ♦  C^a^HT  -  Ct*t(»  -  m*)  -  Cj2iL» 

«y*°  •  Ci2t*  *  cL*y  “  ^cL®l  *  cl2®t^T  _  C12®T(a  -  *i )  - 

where  superscripts  0  and  90  danota  ply  orientations .  It  should  ba  borna  La  mind 
that  a lnc«  ■  -  a(z,c).  all  strata#*  dapaod  on  a  and  c.  However  for  tha  sake  of 
r>ocaCLonal  brevity  thl*  dapendanca  la  suppressed. 

Vlacoalast lc 


In  analogy  with  (8)  tha  vlacoalast lc  atress-etralm  relations  for  combimed 
effects  (10),  expressed  In  terns  of  compliances,  ara 

*LJ«°  *  S|!-,0  *  *1°  t 

S1?o,°  ♦  ♦  D \j  U  -  dr  -  »j0 

J°y 


90 


D  0° 

S,jS ,™  ♦  Sto/"  - 


to  ,90 

,9  ^i_-t  ♦  s,,^90  -  «,* 


♦  D,  ^ 


dt 


where,  in  (9) 


l  -  C(*.t)  -  f  ,a(s)itj  fc,)>  •  C  -  C  <*•’)  '  ^ 

end 

*1°  -  •*  -  d^m  '  aLAT  »  R\90  “  -  m^  )  -  a^T  , 

*2°  -  iy  -  8T(»  -  »|  )  -  o^T  .  Rj’0  -  cy  -  6L»  -  OjAT  . 


«) 


elastic  stkess  analysis 


Consider  a  aynetric  laminate  coapojfd  of  Na  0*  piles  and  NgQ  90*  piles-  Tha 
.eainate  undergoes  n  total  temperature  change  AT,  which  consists  of  cooling  fro* 
ure  to  rooa  temperature  and  then  heating  to  the  desired  level  of  conditioning 
temperature  at  which  moisture  conditioning  la  enhanced.  Tha  laminate  also  absorb* 
noMture  dur  lng  storage  prior  to  conditioning  and  dua  to  conditioning  at  increased 
- jmldlty ■  All  this  time  the  laminate  la  free  of  external  load*. 


The  elastic  stresses  are  obtained  by  superposing  the  solution  for  a  geo¬ 
metrically  constrained  case,  In  which  tha  moisture  sorption  and  temperature  varia¬ 
tion  are  considered  without  any  deformation,  and  tha  solution  for  strain  fields 
*(c)  and  ty(c)  which  ara  uniform  throughout  tha  thickness.  The  superposition  Is 
effected  in  such  a  manner  that  the  net  resultant  forces  Px  and  Fy,  In  tha  s  and  y 
^fractions  respectively,  vanish.  The  requirements  EFX  -  0,  Ify  •  0,  raprasant  tha 
absence  of  externally  applied  load*,  but  ara  insufficient  to  account  for  detailed 
edge  effects,  and  provide  the  equations  to  determine  the  strains  t^Ct)  and  c-(t). 
Mnce  tha  moisture  profile  at  each  time  la  known  from  tha  solution  to  tha  diffusion 
equation,  denote  the  total  contribution  to  avelltQg  due  to  moisture  in,  say ,  tha 
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1th  ply  by 


•Toif)  -  \  /  <■  -  »1 )dt  (10) 

Thlcknaas  of 
1th  ply 


la  (10) ,  ih.  superscript  ♦,  which  My  bs  0*  or  90*.  Indicates  the  orientation  af 
ths  lc"  ply.  Conssqusndy,  ths  cunulstlvs  trsnsvsrso  ownlllas  dun  to  no  Is  turn  for 
all  0*  piles  Is 


.-(hOU))  -  £  s-(M?) 

1-1 


end  slnllsrly  for  ths  90*  piles 


H»0 

-  £  «TCM1°> 

i-l 

Tha  total  longltud las  l  availing  can  be  obtained  in  an  analogous  mannar.  How¬ 
ever,  since  we  cake  0^-0,  thast  analogous  terns  sell  1  vanish  in  our  csss. 

Da not*  h°  and  h90  as  the  total  chlcknaaa  of  all  of  eha  0°  and  90°  plies  re¬ 
spectively,  then  enploynenr  of  (8)  togaChar  with  tha  equilibria  raqulrsnaata 

trm  -  o,  rry  -  o  yield 

(h°CL  ♦  +  (h°  ♦  h,0)C)}cy  -  Rj 

(h°  ♦  h90)C|?c,  ♦  (h0^  ♦  h90Cl>cy  *  *j 

where 

Rj  -  tLcL(M°)  ♦  Cust(M°)  ♦  Ct.t(M9°)  ♦  CwsL«*°) 

<•  h0(Cl»l  Ku>aT)4T  <  h90(CTaT  *  Cw»,  )4T  (12) 

R.  -  Cjj^CM0)  ♦  Ctst(M°)  ♦  C,?aTOI,0>  * 

*■  h0(CTc.T  »  tl2oL>4T  «•  h,0(CLoL  ♦  C1;ST  )4T 


Th*  *  intii  :jneuub  solution  of  (II)  yields  the  valuta  of  tha  t  ine -dependent  atraina 
cx  and  v.  The  acres*  prof  lie  due  to  temperature  and  noiatura  conditioning  is  ob¬ 
tained  by  substituting  the  strains  fro®  (11)  into  (8). 

The  t ins -dependant  elastic  stresses  are  thus  determined  for  all  tints  until 
tsolscui  <  saturation  is  reached  for  each  particular  conditioning  achat* .  At  this 
etage.  the  fenpe rat ure  Is  lowered  down  to  roo®  temperature,  causing  guddsn  in¬ 
crements  In  ail  stresses.  These  increments  are  obtained  by  solving  a  ala^liflad 
version  of  (ID.  (12)  and  (8)  in  which  and  8  contain  only  aultabla  AT  tanaa. 
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VISCOELASTIC  ST* ESS  ANALYSIS 

Consider  Eqs.  (9).  Thee*  equations  contain  four  unknown  itrUHi  whoso  evalua- 
n  requires  the  solution  of  simultaneous  Integral  aquations. 

To  avoid  this  cumbersome  task  wa  Invert  the  compliances  ao4  rowTlte  (9)  to  as* 
ts  streeeaa  In  terms  of  .  Cy,  a  and  T. 

as  a  first  step  toward  this  inversion  consider  tho  Laplace  transform  of  (9). 


Uenot 


log  ?(p)  -  J  e“pt  f  (t )dt  and  p?(p)  • 


1(9) 


we  obtain  f roa  (9) 

,0  -  _  0  -0  »  _  90  *  .  90  r  90 

b,  cr  ♦  S..J  -  R,  ,  5,-0  ♦  S.o_  •  *, 

Lx  12  y  l  12  s  L  y  2 


c  .  90  A  ,  .90  5  90 

V,  ♦  -  *1 


12  x 


♦  - 1 


<n) 


Inverting  (13)  we  obtain 


3  0  .  c  =  0  a  c  5  0 
x  CLR1  ci 2*2 

a  90  - 

X 

+  Cjjlj90 

.  0  -  -  0  0 

%  -C12*l  >CrR2 

_  90 

o  • 

y 

elA« 

*  V290 

-  ST 

_  fl2 

where.  In  (14)  - 

' 1-1  *  T 

•  C12  ‘ 

6 

a*) 


,nd  4  *  SLST  '  S12  ' 

In  the  prtiint  casa,  with  and  512  being  conataat,  tha  only  Cara  which  ra¬ 


in  vtfw  nf  (5) 


mnarqurnt  ly 

i  .  u  .  D  ^ 

T  0  1  q 

P 


u» 


«r  .hall  capioy  thr  .ipproalaatc  Ion  111],  wharrby  r(t)  caa  ba  axpraaaad  fraaa 
.irvnrdlni  to  F ( t )  -  F ( p )  p  -  l/2t  Tharafora  (15)  ylalda 


.'jU  I  ’  t>0<l  »  *t  ’) 


<H> 


Siailarly,  InaCaad  o'  •.  or  shall  wrlta,  upon  lavaraloa 

r  2 


(17) 


r  trirtfona  tnvaralon  of  (9)  will  give 
bib 


Altogether,  the  L 


dl?(t )  »*»2 

(O  -  l  (t  )  1  — ♦  C.,|((t)  -  l<t)] 


-)■ 


a>>  -  /  kit 

o  ' 

0  (l  (  «?<»>  <"£<*>) 
«“('»-  J  |C12l(<t>  -  t<t)l  — It—  ♦  CTU(t)  -  t(r)l  |di 

0  1  ' 

90  90 

and  similar  expressions  for  (t)  sod  a ^  (t). 


(IS) 


btpllclcly  In  (18).  va  hava  C^t)  -  D0(l+Atl')/4  (t) .  CjU)  -  »LMCt), 
and  C12(t)  -  -  S12/d(t)  uhara  4(0  -  SjDgU'WkO)  -  s[2 

Just  as  In  the  alaadc  analysis  wa  now  hava  tarva  Rj(t),  R^t),  alac 
tfOU)  and  which  contain  known  valuaa  of  molature  m(*,t)  and  temperature 

AT(t)  and  yet  unknown  •trains  cx(t)  and  ty(t)  which  art  spatially  uniform.  Those 
unknown  atraina  aye  praaent  in  order  to  raaaaura  null  resultant  forces  £Fx(t)  • 

0,  IF  (t)  -  0.  Note  that  the  above  requirement  of  null  resultants  which  applies 
for  all  tinea,  provides  the  conditions  for  the  determination  of  ex(t)  and  «y(t)  at 
all  tinea. 


In  addition,  we  should  note  that  the  viscoelastic  analysis  considers  tbs 
Initial  Stresses  sod  atraina,  cK(0)  and  c_(0),  and  accounts  for  their  own  relaxa¬ 
tion  as  time  progresses  during  cha  conditioning  stage. 


Finally,  In  order  to  clarify  natters,  we  observe  that  during  conditioning  the 
moisture  •  va* lea  but  rhe  temperature  T  remains  constant,  therefore  AT  •  0  in  Ij®, 
etc.  during  that  stage. 


Altogether  we  hava  o  (*.t) 


■j[,LCt<<>  *  •**«>]". 


cut*  *cood^ 


•u.o)  -  »0  [ » yo '„«»  ♦  Cu«)  iy(0> 


<U  (t ) 


d»  (T) 


+[»LCL(t)  f  6TC12(()]  [ 

*  /  y<«>  - 1(,)>  ~dT“ ♦  c(2u<t)  - C(T)I  -Jr 

O  ' 

-  t (t ) J  -  oTcl2U(t)  -  t(t)l  ~ 


dr 

(19) 

,,0(..O. 

For 

dataila  sea  Baf •  !  1 ) . 


In  fl9)  m(a,0)  la  the  moisture  level  at  toe  beginning  of  the  condition log 
stags.  This  moisture  level,  et  t  -  0* .  la  taken  to  be  the  ambient  conditioning 
K>let<>re  level  at  the  outarmoat  surface  of  cha  laminate  (i  ■  ±h),  while  the  re¬ 
mainder  of  the  laminate  remataa  at  the  initial  moisture  level  due  to  storage. 

f.q .  (19).  and  Its  counterparts  for  oK^0.  oy°.  and  oy*®*  express  cha  tlma- 
dependent  viscoelastic  atraaaaa  until  saturation  time  t  •  Cj.  The  value  of  tf  of 
course  depends  on  the  condltlooing  temperature  T.  At  t  •  tf  the  teopareture  la 
Lowered  to  the  room  temperature  TR,  causing  sudden  etreaa  Increments. 

These  Inc  remanta  ate  computed  elastically  because  of  their  rapid  tfsvslopneat  t 
which  rules  out  v  lacoelaat ic lty . 

a 

Thus,  at  time  -  t^  wa  gat 


SJA 


-  v 

(M) 

♦  Cu<0)|cy(c*)  -  ^(c"»  ♦CL(0>t«k(t^)  -  S(t7» 

tMlofou*  B»pr*s»luo»  cm  b«  d.rlvM  for  o,*®(«,t().  Oy'hi.tf),  *nd 
rsv  detail*  Kef.  (1]. 


Tbe  tiw  dependant  atralna  a, (t)  and  c,(t)  art  da t era load  fro*  the  requirement 
/*  fH  T  . 


i- 


(j,t)dz  •  0  and 


i 


oy(*,t)d*  -  0  at  all  tines  t. 


THa  thus  da taro load 


S  )  and  ty (t )  ara  than  reinserted  Into  (19),  (20),  and  tha  analogous  expressions, 
tc  laid  tha  actual  strata  d iatrlbut loot . 


In  or c*r  to  per  for*  tha  computations  lndlcatad  In  (19),  v#  diacratlma  tha 
thickness  of  tha  lanlcuite  into  portion#  Az  and  tha  tin#  in  tarsal  Into  portion#  it. 
n  thickness  discretisation  la  required  in  order  to  d#tonilo#  tha  moisture  prof 11a 
Ach  elm  and  account  for  ita  of  face  on  tha  “reduced  t  ia»",  (  •  ((s,t). 

..hernora,  chta  diacrat  l  ration  la  required  when  par  forming  tha  integration  of 
«r razee*  acrosi  tha  thickness  to  obtain  ZT%  •  0  and  tty  *  0.  Satisfactory  results 
ware  obtained  by  dividing  aach  ply  into  flea  aqual  lDcraacnta  As  •  0.00023  cn. 


t-«  currant  value  c  la  required  to  evaluate  the  tlae  integral#  Ilk#,  say, 
dc  (t) 


/  CLIUt)  IClJI  ~r~*r 

0 

for  cbla  purpose  we  omploy  a  schema  similar  to  (12).  Denoting  t^  *  0,  t  •  ^  ^ 


)  * 


/dt  It)  r  k>i  dt  It) 

V«(°  - <<T)1  ~l~  *'  -L  J  C'-It(t)  * <(,)1  -fr*  4t 

0  k-1  t 

rl  y!  |CLlCl<J»l'  *  C(tkU)l  *  CLl^U)*l)  '‘^k+l*  '  ‘^k^1 

«-l  '  * 


(U) 


Expression  (21)  contains  all  previously  known  atralna  cx(t.),  t^tj) . 

,i[.)  and  tha  current,  yet  unknown,  strain  €*(1.*,.).  In  thla  Mannar  tna  unknown 
*x^tj4j)  sod,  similarly,  ty(tj^)  «re  Isolated  and  solvad  for  in  a  tins- 
•arcMag  sequence. 

•at  lifer tory  accuracy  was  obtained  by  taking  six  aqual  increment a  per  dacad* 
•*  ag  tha  logerlcholi:  cine  scale. 


hatqial  propertied  and  computations 

Tha  material  properties  which  entered  the  actual  conputatlooa  are  auaarliaf 
-•  Table  l  below.  Du e  to  che  Incompleteness  of  the  experimental  data  It  was 
”t;<»»»*ry  to  employ  data  for  two  graphite )*poxy  aystana.  Data  narked  with  a  star 
*)  refers  to  T3Q0/S208  I'roa  Ref.  (7),  otherwise  the  data  rafara  to  AS/3301-4  fro* 

*.V  (9). 

The  actual  computations  ala  at  comparing  tha  stresses  for  two  conditioning 
*=v: consent  a.  In  both  cases,  tha  Initial  stresses  are  Identical,  sod  are  ettrlbu- 
•d  to  tha  following  causes: 


Parameter 

Symbol 

Ma|t|l  *«••!« 

Holatura  dlffualvity 

Ai 

0.016713  c«,/»«c.* 

0  -  A^xpl-a  /T)  (T  In  *lt) 

»i 

*818. 0‘8* 

Coafficiaoc  of  nolatura  expan  a  len 
par  IX  weight  gain 

Longitudinal 

®L 

0.* 

Tranavaraa 

»T 

0.6099* 

"l 

0.5*50* 

Coafficiaoc  of  thermal  expansion 

Loogitud Inal 

Tranavaraa 

»T 

a 

-.5*  *  10**  c»/c»/*8» 

250*  *  L0*‘  cmJct/'V 

initial  tranavaraa  compliance 

•»2®?®  - 10 11  »T?riT 

0Q  -  aTU  ♦  bT  ♦  CH  ♦  d 

b 

.1*1*7  x  10-’  Tj— 

c 

.29075  «  10*’ 

d 

4. *82  ii  io  sr«_1 

Transverse  compliance,  creep 

°i 

.1*51  x  10*® 

Power  law  exponent 

q 

0.18 

Longitudinal  compliance 

SL 

.7653  x  lO'®**.'1 

"Cross  effect"  compliance 

S12 

-.2679  x  10‘*  IlP*'1 

Temperature  shift  factor  function 

*3 

-105.3 

tq.  <61 

«A 

6.185  x  104  *8 

*5 

-9.055  x  10*  *K2 

Moisture  shift  factor  function 

*6 

.06336 

•j  «  »fc  «  ( Xm) 

bl 

-8.9*2 

Moisture  to  R.H.  conversion  factor 

■  •  r.  «  <m.h.) 

E 

.01505* 

Table  1.  Material  Properties 


1.  Cool  -down  from  a  cure  temperature,  T  •  4)0*K  (3>0*F)  to  room  tempers- 

turf,  T  -  297*K  (7 S'F)  CUr* 

room 

2.  Uniform  moisture  saturation  during  storage  due  to  iOZ  relative  humidity. 

In  a.tuallty,  the  fMect  of  the  r « mpe rot ur e  rise  fro*  Trpt>|B  ■  75*F  to  the  condi¬ 
tioning  temperature  In  Incorporated  Into  the  initial  conditlona. 

Thf  two  conditioning  environment#  ore:  * 

1.  3  )9*  k  (liO'U  at  V8t  R.H. 

:  ,  3V,*K  U80*F>  a  r  981  R.H. 

"i r »*  magnitudes  I  material  parameters,  which  ere  dependent  on  the  conditioning 
environment,  are  calculated  for  each  environment  and  eunarlted  In  Table  2. 

Result!  are  ah»vn  f.r  [O^/VO^I^  lamlnatte.  Further  results,  for  other  lay-ups 
are  given  in  1 1  )  . 

*3» 


Coa4SUo«tBt 

Ttnoroturo 

Paraaeter 

Symbol 

IJ»*K 

3JJ*K 

Umtu 

Molatura  diffualvlty 

D 

.1000  *  10** 

•K 

•  17ft  a  i<f* 

cm’/Mt. 

Initial  transverse 

Do 

.10*3  *  10'4 

.1102  »  10-‘ 

Temperature  shift 

*1 

0.0137 

0.000*3 

factor  function 

elastic  moduli 

S. 

CT 

C12 

131.1  *  10* 

9. *5  x  106 
3.303  x  10* 

131.8  x  106 
9.17  x  106 
3.203  x  106 

k7l 

k7> 

17i 

Table  2.  Conditioning  Dependent  Material  Properties 
RESULTS  AND  DISCUSSION 

Figs.  )  and  4  exhibit  the  Jtteaaea  ox(t,c)  and  0y(*,O  due  to  aoleture  sorp¬ 
tion  at  93t  R.H.  These  figures  provide  cooper  Isons  between  conditioning  at  339*K 
tJj^hed  lines)  and  at  35S*K  (solid  lines)  sod  deeon strata  the  "competition"  be¬ 
tween  tha  affects  of  moisture  and  temperature  as  stress  Inducing  agents  on  ooe  hand 
jnd  stress-relaxing  parameter*  on  the  other  hand. 


10  34 

0 

10  34 

10  34 

0 

0<mPa) 

10  34 

•10.34 

0 

Ot  mP«l 

10.34 

IKSi ) 

msn 

(KS>1 

Fig.  1.  The  Viscoelastic  Stresses  ox  Vs.  The  Spatial  Coordinate  x  et  Various 
.cirs  t  for  |0}/90-.1B  laci’Mtr  (r  -  0  la  Symmatry  Plana)  Due  to  Moisture  Condltion- 
ng  at  98T  R.H.  at'm'K  (Solid  Lines)  snd  339-I  (Dashed  tinea),  o  In  aft  (ksl  la 
rent heses) . 


S 19 


Fig.  4.  Tha  ViacotliacU  StrtMM  o  Ve.  Tlim 
Spatial  Coordinate  t  at  Various  T Imi  t  for 
{02/902).  Laminate  (*  -  0  is  Symmetry  Plana)  Due  to 
Moiatura  Conditioning  at  98X  *.K.  at  353*K  (Solid 
Lines)  mad  339*K  (Dashed  Llama).  o  in  a?a  (bet  la 
Parenthesea) . 

0  M 

To  i4 


To  appreciate  the  total  affect  of  conditioning,  consider  the  atraaaaa  oy  at 
each  c  -nriit  toning  temperature  {Fig.  4).  During  the  early  atagaa  of  conditioning, 
t  -  1  minute,  for  which  there  hae  bean  but  minimal  added  moisture,  the  stress  pro¬ 
files  reflect  the  difference  betveeo  the  two  conditioning  temperatures,  Bowmvar, 

1a  conditioning  progresses,  say  up  to  1000  minutes,  tha  dependence  of  moisture 
sorption  and  tha  etresa  relaxation  on  the  conditioning  taaparatura  becomes  promi¬ 
nent.  At  chla  atage,  there  la  noticeably  more  moiatura  sorption  at  355*K  that  at 
339*K  am  la  evident  by  the  larger  compressive  atraaaaa  in  tha  trsnsvarts  direction 
of  tha  outer  ply.  Yet,  at  the  outer  surface  0/  Che  laminate,  where  tha  bouodarlsa 
wire  exposed  to  the  ambient  conditioning  moisture  level  throughout  tha  entire  1000 
■Inutea  of  the  ;ond 1 t loning  atage,  tha  conpreeetv*  atraaaaa  ara  about  17  psreent 
.ower  for  the  conditioning  environment  of  After  10,000  minutes  we  note  that 

tha  slopes  of  the  transverse  stresses  In  the  outer  ply  ere  of  opposite  signs.  This 
contrast  is  attributed  to  the  enhanced  relaxation  ac  333*X.  Finally,  upon  reaching 
oulatu;**  saturation  (’  22  days  at  ISb’K,  and  43  days  at  339*K),  the  strata  profilm 
In  b>ch  the  longitudinal  and  transverse  directions  la  nearly  uniform  across  tha 
thickness.  yet,  the  fact  that  tha  screes  magnitudes  are  similar  !•  purely 


fig  V  Tv*»  "treasti  ’Y  V*.  Th*  Spatial  Coordlosta  *  at  Various  Haas  t 
i>u  f  to  C  >nJ  1 1  •  on  in*  at  98*.  PH.  at  )5Vk.  Vi^coalaatic  Cssulta  (Solid  Linas) 
Vs.  LLoaar  Clast  tr  Pradlitlons  Hathtd  Lines).  <3  la  afa  (ksl  In  paraatfcaaat) . 


co lac ldentel . 


Th*  final  cool-down  to  took  taskperaturt  of  297*g  superimposes  additional 
atraaaaa  on  the  laminate.  However,  slnct  conditioning  at  a  lover  temperature  cot' 
responds  co  a  smaller  rtraaa  increment,  cha  atraaaaa  dua  to  conditioning  at  339*11 
ete  lower  than  those  atraaaaa  resulting  froa  conditioning  at  335*11.  Note  alao  that 
tha  location*  of  tanalla  and  coaprasalva  atraaaaa  ara  ravaraad. 

For  comparative  purposes,  the  alaatlc  and  viscoelastic  atraaaaa  ara  presented 
in  Pig.  3  for  conditioning  at  355*K. 

Tha  figure  ahova  that  by  discarding  time-dependent  raaponaa  and  overlooking 
street  relaxation  tha  alaatlc  analysis  pradlcta  atraaaaa  whtch  differ  from  tha 
vlacoalaatlc  raaulta  la  algn  and  ovaraat inataa  than  by  thraa  to  six  fold. 
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REMDbAl.  thermal  stresses  in  an  asymmetrical  cross-ply 

GRAPHITE/EXPOXY  LAMINATE 

Brian  Douglas  Harper* 
f.  Weitsman*’ 

Texas  ASM  Jniversity 
io 1 1 ege  Sta t : . n ,  'exas 


’his  oaper  :rese"t'.  an  exploratory  investiga¬ 
tion  of  pre  ’'esidual  t’esses  m  AS-3S0?  graphite* 
epox,  lami'jtes  due  to  cool-down  from  tneir  cure 
temcera turn  .  Emphasi  .  ■  s  pla.ed  on  the  slgnif1- 
cance  of  t : me -dependent  material  behavior  and  the 
potential  utilization  of  'his  phenomenon  to  reduce 
residual  stresses  ;iy  a  judicious  choice  of  the 
cool  down  process,  ’re  analysis  considers  the 
time -dependent  behavior  or  the  material  and  all 
calculations  employ  recent  data  on  the  thermo- 
viscoelastic  response  of  the  AS- 3502  graphite/ 
epoxy  s/stem.  'h*  viscoelastic  analysis  is  veri¬ 
fied  through  curvature  measurements  of  unsymmetric 
-ross-ply  plates  •'aoncated  from  the  AS  -  3502 
iraohite  epoxy  material. 


.  ■ Production 

’re  curinq  ev  composite  laminates  at  temp¬ 
eratures  above  tneir  service  temperature  induces 
-es’Jual  stresses  uaon  cooi  down  oue  to  the  dif¬ 
ference  in  thermal  expansion  coefficients  in  the 
'ongitudmal  ana  transverse  directions  of  the 
'am’rae.  It  "as  long  been  recognized  that  the 
presence  of  residual  stresses  greatly  affects  the 
strengtn  of  compon'e  laminates.  For  this  reason 
there  nas  been  an  increasing  interest  in  develop- 
-ent  ..’f  improved  methods  to  determine  these  stress- 
-  a r J  to  Tnir-zp  t~eir  magnitudes. 

;,nano  and  Hihr.l  utilized  an  approximate 
e’-  d  'or  iBterm.pir.;  -r.e  residual  .tresses  in 
•  r  w-met'"  .  a  i  csi-clz  Iminate.  ’hey  uOnsid- 

e  , -  m  •'  t'e  'a-rnats  at 

•e •  r  -  .  wm'.p  above  '"at  'emperature  ~0 
■>/  as',  me!  i  cere  tress  >'.ite. 

■*  ’as  v»r  ,,-,q  noted  fat  epoxies  exmpit 
pr-jnoonced  v  ..coelastic  effects,  consequently  there 
xr  its  pom  need  mo  'ustif'cati'i*  *or  an  analysis 
which  con , i ders  the  t  me  dependence  ot  the  material 
:.jr'"q  tne  entire  cool  down  orocess. 

)ouglass  and  Wei'smanu  useo  a  viscoelastic 
analysis  to  predict  the  residual  stresses  in  a 
cyimetric  . 'tss -pi  v  laminate  due  to  both  temper¬ 
ature  and  moisture  effects;  however  no  attempt  was 
made  to  verify  these  'esults  experimentally. 

s  caper  presents  a  method  for  evaluating 
tne  residua!  stresses  in  unsymmetrical  cross-ply 
laminates  which  employs  linear  viscoelasticity 
•-roudhtout  the  tool  down  stage  'Jnsyroietric 
'aminates  are  chosen  because  "heir  deformation 
involves  anti -clastic  curvatures  which  are  amenable 
*c  measurement ,  unlike  the  case  or  symmetric  lamin¬ 
ates  -everai  unsynmetr’c  c'oss-Dl  /  AS  -  3  502  graD- 
ni  te  ■  epoxy  laminates  were  fabricate.;  and  the  cur¬ 
vatures  measured  to  .  e ' -  * ,  -he  viscoelastic  analy- 
si  s 


•graduate  utuflent ,  Mecbari :a  •  Engineering 
**<’ro'es  sor  ,  Mecnar'Lu  4  Ma  t  e  r  1  u  )  ,  /enter, 
Snjj  i  r.por' eg  eparfment 


Due  to  the  significant  creep  of  epoxies,  esp¬ 
ecially  in  the  elevated-temperature  range,  there 
exists  an  opti.nal  cool  down  path  that  will  minimize 
the  residual  stresses.  This  optimal  path  is  deter¬ 
mined  for  the  AS- 3502  graphite/epoxy  system  using  a 
method  suggested  by  Weitsman’. 


Elastic  Stress  Analysis 

The  thermoelastic  stress-strain  relations  for 
an  grthotropic  lamina  under  plane  stress  conditions 
ar»* 


CijUJ 


i,  j  =  1,  2.  6  (1) 


In  (1)  and  the  sequel,  r  is  stress,  e  is 
strain,  Clj  are  stiffnesses,  a  are  coefficients  of 
of  thermal  expansion,  T  is  the  temperature  and  cT 
the  temperature  difference. 

Consider  the  cross-ply  0n/90n  unsymmetrical 
laminate  shown  in  figure  1. 


Fig.  1  Geometry  of  the  Unsyitmetric  Cross-ply 
Laminated  Plate. 


Since  in  the  principal  directions  aw 2  vanishes 
and  there  exist  no  shear  stress  to  normal  strain 
coupling  terms,  we  have 

Y  *  «  a  r  a  H 

xy  xy  lxy  u 

Referring  to  the  individual  laminae  in  Figure 
1  we  can  express  the  stress-strain  relations  as 

fol lows  - 


■0  appear  in  the  proceedings  of  the  22nd  SSD  &  Materials 
Conference,  -"Itanta,  Georgia,  April  6-8,  1981. 


TOT 


lA  *  .:‘r'  " 


90 


X 


(0) 


op 


!  n  c  • ,  i  s  upc  r r  c  "  ’  f  s  '•  ikM  c  a  t  c  y  »  / 

one^tat’or.  and  subscripts  l  and  T  denote  tbe  prop¬ 
erties  in  directions  parallel  and  transverse  to  the 
fibers  respectively. 

Following  the  assumptions  of  classical  plate 
theory*  *e  can  express  the  strains  In  terms  of  the 
roidsurface  displacements  u  ,  v  ,  and  w^  and  get, 


(3) 

« V* 


denoting  mi  dp  lane  .trains  and  curvatures  in 
the  usual  manner 


equations  ( 3 }  become 


X  1  V  ‘  :kx-  -y  *  y "  *  ^  (5) 

"he  midnisne  strains  and  curvatures  in  the 
1/90  unsyrrinetr ica  1  laminate  due  to  a  uniform  temp¬ 
erature  change  :7  can  be  determined  bv  requiring 
that  the  hp*  resultant  forces  fN'  and  bending 
•"omer-ts  V'  wn-ch  act  on  toe  plate  must  vanish 
"'is  -esults  *n  the  fol lowing  expressions 


N 


-  h 


>6) 


1  *  V  “  V  It  (CL*V  ?C12>  p 

~2  '  '  Q  l T ,  hK  *  -hK 

*“  J  *  y 

5  5['  VCt*?V  '  '  (Cl-Ct)p]at 

«ne»-e,  in  ?' 

:  1  1  [<V  CT>"  -Sj  -  f  iCfCT'2 

0  •  (Cl-C]2)Jl*  tCT+C12)aT 
3  *  (3[_-Cj2)  Ti|_ _ ( Cj *C-| ^ I  Jj 


(7) 


(3) 


Stresses  may  now  be  evaluated  by  incorporating 
the  results  from  equations  (7)  into  equations  (2). 


Viscoelastic  Stress  Analysis 


It  has  beer  noted  that  graphite/epoxy  compos¬ 
ites  exhibit  a  considerable  amount  of  time  depen¬ 
dent  mechanical  response,  especially  at  elevated 
temperatures.  This  time  dependence  may  be  approx¬ 
imated  by  a  viscoelastic  constitutive  relationship 
and  employed  to  predict  strains  and  curvatures  due 
to  cool  down  from  the  cure  temperature. 

Although  the  behaviour  of  polymeric  resins 
imiedlately  after  cure  exhibits  dependence  on  such 
complicated  factors  as  aging  time,  quenching  rates 
and  perhaps  on  a  variety  of  non-linear  effects  we 
shall  assume  herein  a  linear  and  thermorheological - 
ly  simple  (TRS)  viscoelastic  behaviour. 

Data  on  graphi te/epoxy6  indicate  that  the 
assumption  of  TRS  behavior  Involves  only  small 
errors  . 

Accordingly,  the  transient  thermal  response  can 
be  related  by  means  of  a  single  temperature-depen¬ 
dent  function  A.  T1  which  is  called  the  "shift- 
factor’  . 

The  time-dependent  portion  of  the  viscoelastic 
behavior  involves  reduced- t imes ,  which  are  denoted 
in  tne  sequel  by  t(t!  and  ;(r),  where 


tit) 


S(t)  1 _ ds _  (9) 

J  5  [T(s)J 

n  1 


The  viscoelastic  counterparts  of  the  elastic 
expressions  (2)  are: 


Jt  V  T  > 

X 

d-r 


*  c12[;(t! .  T'i] 


-  iLcLU(t:  -  -.(•)]  *  iTc12[c(t)  -  •.(-!] 

d^T(x)  I  d- 


;omb  i  n-j  equations 


9 ) ,  f S i  and  ( 6 )  we  get 


(10a) 


dc  ,  s  j 


/'.t)  =  I  '  |  -«(*)] 

0  '  J,  (x) 

*  C,[--.(t)  -  ■; ( t ) ] 


-  i(r)]  «•  iTCTr-:it) 


;  Ub  , 


,  -  .  .  1 0  .  , 

** '  so »  \  <*n<J  -  it,1  car*  oe  *  • 

7  ”  ( t )  ana  .-'’(t),  respectively,  by  interchanging 
t^e  roles  of'  and  c  in  (10' - 

i  V 


where,  in  (12),  D[C(t)  -£(t)],  P[r(t)  -f(t)],  and 
Q  CC( t )  - *" ( t ) ]  are  obtained  from  D,  P  and  Q  in  (8) 
with  moduli  C.  replaced  by  -[£(t)-£( t)] . 


In  the  AS - 3502  graphite/epoxy  system  used  in 
this  work,  the  only  compliance  which  exhibits  sign¬ 
ificant  time  dependent  behaviour  is  the  transverse 
compliance  ST°.  The  time  dependent  transverse  com¬ 
pliance  ST,  dan  be  expressed  by  a  "power  law"  equa¬ 
tion  of  tCe  form 


(13) 


wnere  t  is  time,  and  0.,  q  and  t,  are  material  con 
stants.  furthermore ,  the  shift-factor  is  given  by6 


In  the  case  of  thermorheological 1 y  simple  oe- 
navior  we  can  obtain  viscoelastic  solutions  from 
the  results  to  analogous  elasticity  problems  oy 
means  of  tne  so-called  correspondence  principle. 

The  principle  remains  valid  if  at  any  qive  time  the 
temperature  is  spatially  jrnfom.  In  our  case  this 
condition  is  satisfied  to  a  i„fftcient  degree  of 
accuracy  Jue  to  the  mqh  thermal  conductivity  of 
qrauhi  re,  ecnxy  laminates. 

this  oaper  we  ■'nail  employ  the  correspon- 
lence  principle  in  conjunction  witti  an  aDProximate 
tecnnique  for  the  inversion  of  the  Laplace- 
transformed  elasticity  solution.  Applied  together 
those  principle  and  technique  combine  to  yield 
the  quasi -elastic"  method.  '  Accordingly,  the 
viscoelastic  response  function  to  a  unit  input  is 
aooroximated  :>v  an  elasticity  solution  in  which  all 
elastic  constants  are  replaced  oy  tne  corresponding 
reduced-time  lepenaent  properties.  Once  this  unit 
response  function  is  obtained,  the  response  to  any 
general  input  may  :e  obtained  by  means  of  the  con¬ 
volution  integral 

'.msec, ient  ’  , .  we  obtain  the  *o!  lowing  expres- 
s '  c r  s  ‘o-  •ho  time  ieo“ooen‘  midpl.inn  strains  and 
■|.r  ,a  *  .  rn.  •  r  .  ’  Sw  Oi 


At(T)  ■  exp(-T/A  *  B) 


(14) 


where  T  is  the  temperature  in  degrees  Kelvin  and  A 
and  3  are  material  constants. 

Since  equations  (10)  -  (12)  contain  stiff¬ 
nesses  rather  than  compliances  it  is  necessary  to 
invert  the  compliance  matrix  to  obtain  the  time 
dependent  stiffnesses.  Performing  this  inversion 
we  get 

cT(t)  =  HIT-  ci2(t)  *  5TET 

(15) 


(t) 


$T(t) 

TTtT, 


where 


u(t) 


SLS. 


.(t)  -  s 


12 


In  order  to  perform  the  computations  Indicated 
in  (12)  we  discretize  the  time  domain  into  portions 
it  between  ti  *  0  and  the  current  value  t.  Consider 
for  example  the  convolution  integral, 

/*  AL-.U)  06) 

-'o 


Cenoting  t1  *  1  and  t  *  t().,  (16)  reduces  to 


A 


(t)  -.;(:)]  d’  ’ 


£  / 

k-1  -C 


t  km  1 


f[c!t) 


-  C(T)]  d:  =  l 


*  (17) 

i  b(V!> 

K*1  ( 


c(tk  +  1)]  *  FU(tj+1)  -  C(tk)]j 


OTft^,;  -  AT(tK) 


:  •  c.fctt;  -•»•)) 

.■  -'I  j  j 

’  -  -! t  -•;(-)  J 


Optimal  ~ime-Temperature  Path 

Reverting  to  compliances  in  place  of  stiff¬ 
nesses  in  (7)  and  employing  the  quasi -elastic 
approximate  method  we  obtain  the  following 
expression 


:(.)]  “Jjil  J- 


in  (18) 


6[J!>Li  _  n T ) S T ( t )  -  2Sp(xL  -  nT)] 

SUt)  *  14S.TTTt)  ♦  S‘  -  16S,'  jl9) 


a:  ,o  :.H  ■ '  ■■  T(  • )  -  r, 


where  T,  :s  the  stress  fro*  reference  temperature 

The  residual  stresses  in  an  elastic  material 
depend  strictly  on  the  temperature  difference  AT, 
and  not  the  time- temperature  path  employed  to 
achieve  that  temperature  difference.  A  viscoelast¬ 
ic  material,  however,  has  tune  dependent  material 
properties  which  ,ire  strpnqly  affected  by  temper¬ 
ature.  Since  the  relaxation  of  residual  stresses 
are  enhanced  at  elevated  temperatures ,  there  exists 
an  optimal  coo!  down  path  T(t!  which  minimizes  the 
residual  stress. 

With  prescribed  'nitial  ^elevated' temperature 
’inai  temperature  Tp  and  cooling  time  t^  a  sol¬ 
ution  for  the  optimal  path  ’it!  can  be  obtained  for 
a  wide  class  of  ‘unctions  K ( t }  in  t 1 3 ) .  We  take 
*j  to  be  cure  temperature,  at  wnich  the  laminate 
is  assumed  stress  free,  and  Tp  as  the  room  temper¬ 
ature. 

the  optimal  path  *(t'  which  minimizes  the  cur¬ 
vature  k  can  be  derived  bv  a  method  developed  by 
Weitsman^. 

Accordingly,  it  can  be  shown  that  this  optimal 
path  contains  discontinuities  at  t  «  0  and  t  •  tf , 
"he  first  temperature  drop,  from  *.  to  7„,  is  deter- 
minted  from  the  ‘-jnscondental  eduation 


•  r  as  exponential  temperature  shift  factor 
•’ unc •  >cn  as  in  rogation  ,141,  equation  (TO)  re¬ 
duces  to  T„  -  -  5. 

iii  the  time  interval  between  the  initial  and 
f'nal  discontinuities  ’t  was  shuwn  that  the  optimal 
path  T(r!  is  joverned  by  tne  equation 


:  (t,)  -lit!]  AT[T(t)] 
-•;(t"'l  A-i  'li  t '  ja"t 


in  ..'O'  ana  ,.'11,  primes  'ridicate  derivatives  with 
rpspect  ‘o  ‘he  lrqumenf 

'he  iDtimal  time-temperature  path  is  constiuct- 
"d  by  iterat'nn  through  • he- emp I ovment  of  (21). 
Jividioq  the  ‘imp  period  into  n  equal  SUb- 
■mterva'is  t  the  tera'ior  on  (21'  may  be  carried 
Out  bv  selecting  a  juess  value  T(tf)  denoted  by  I, . 
at  ’  '  equation  qives 


and  the  discretization  of  tf  into  increments  At 

yields 


T(tf  -  At)  .  T  t  — -L-L -  IM 

ATtTf )A-(Tf)  F  W 

With  both  Tf  and  T(tr  -  At)  known  we  extra¬ 
polate  (21 )  back  to  t  -  if  -  2 At  to  determine 
Ti'-c  -  2dt),  Continued  exprapolation  of  (21) 
back  to  t  •  o  yields  a  value  of  T(o)  that  will 
generally  be  different  from  T„  determined  bv  (20). 
If  t;o;  ■  T„a  new  guess  value  of  Tf  is  selected 
less  than  the  original  Tf  anq  vlce  versa. 

This  iteration  is  continued  until  a  value  of 
%  Is  found  which  gives  T ( 0 )  c  T0  to  some  desired 
atcuracy. 

The  value  of  T,  thus  determined  will  generally 
be  different  from  the  room  temperature,  which  there 
by  determines  the  discontinuity  In  the  optimal  time 
temperature  path  at  time  t  ■  t,. 


Summary  of  Material  Parameters* 


The  properties  of  the  AS -35- J  graph! te/epoxy 
system  used  in  this  study  are  presented  in  Table  1 


Table  1.  Material  Properties  . 


Parameter  Symbol 

Magni tude 

COMPLIANCE 

longitudinal 
"Cross  Effect" 

sL 

.552 

Sp 

-  .160 

Transverse 

sj* 

6.12 

COEFFICIENT  OF 

THERMAL  EXPANSION 

Longi tudlnal 

-  ,27 

Transverse 

aT 

21  .8 

’EMPERATURE  SHIFT 
r ACTOR  FUNCTION 

Equation  (14) 

A 

6.0 

B 

49.7 

TIME  DEPENDENT 

TRANSV.  COMPLIANCE 

Equation  (13) 

T„ 

.222 

ror  0  <  log  t  <  5.8 

ID? 

.618 

Id1 

.00667 

For  v  <  log  t  <  8.4 

ID, 

.506 

Id’ 

.0217 

log  t  <  8.4 

IP, 

.305 

lq’ 

.048 

In  Table  1,  all  compliances  are  in  10” 
in  JK.  TS  is  inin-6  ps  1  _1  min. "A  and 


A  is  in  JK,  D|  is  in  10 
in  minutes 


*  To  check  tne  fabrication  process  of  our  speci¬ 
mens  we  determined  experimentally  the  values  of  E,  , 
as  well  as  .’  u  our  moduli  were  slightly 
higher  than  those  recorded  In  literature.  We  also 
found  -i  ult  »  233.5  ksi  »  5,89  ksi  as  com¬ 

pared  to  218.4  ksi  and  6.31  ksi  ,  respectively. 


Experimental  Determination  of  Curvatures 

Unsymmetric  0°2/90^2  6"  square  plates  were  made 
from  AS-3502  graphi te/epoxy  material.  These  spe¬ 
cific  dimensions  were  selected  in  order  to  obtain 
thermal  strains  and  consequent  out-of-plane  deflec¬ 
tions  which  were  sufficiently  large  for  ease  of 
measurement,  and  yet  small  enough  relative  to  the 
plate  thickness  to  maintain  the  validity  of  the 
xinematic  assumptions  of  classical  plate  theory. 

Twenty  such  panels  were  manufactured  and 
mtasured.  These  panels  were  divided  among  five 
groups,  each  group  undergoing  its  own  cool  down 
history  from  cure  temperature  to  room  temperature. 
These  time-temperature  histories  are  shown  in 
Figures  2  and  3. (Note  that  in  Fig.  3  the  time  is 
given  in  hours,  while  in  Fig.  2,  which  shows  the 
fast  cooling-path  ' A" ,  the  time  is  in  minutes.) 

Upon  cool  down  tne  curvatures  of  the  plates 
were  determined  by  securing  the  plate-specimens 
to  three  fixed  supports  which  defined  a  reference 
plane.  The  deflections  o£  the  plates  relative  to 
tnis  reference  plane  were  'hen  measured  using  dial 
qauges . 

1“  view  of  the  last  two  of  equations  ;4)  we 

ha  .  e 


Figure  3.  Cool-Down  Paths  "B"  thru  "E" 

The  locations  of  the  points  where  w  was  measur¬ 
ed  on  the  surface  of  the  plate  are  shown  in  Fig.  4, 
where  all  distances  are  in  inches. 


w  =  i  -  l  k/  (24) 


The  measured  values  of  w  thus  enable  use  to 
calculate  K  and  K  .  The  curvatures  were  cal¬ 
culated  bet/*  with  Knd  without  the  assumption  that 
k 

'  y 


Time  (mmules l 

Figure  3.  Cool  Down  Path  "A" 


y 


Figure  4.  location  of  Points  where 
Deflections  were  Measured 


Computations 

The  computations  performed  In  this  study  were 
aimed  at  two  purposes: 

1.  To  compare  experimentally  determined  cur¬ 
vatures  to  those  predicted  by  (12)  employ¬ 
ing  the  paths  A-E  shown  in  Figs.  2  and  3. 


1 


2. 


To  determine  the  optimal  time-temperature 
paths  that  will  minimize  the  curvatures 
for  cooling  times  of  50,  100,  and  200 
minutes . 

All  calculations  employed  the  material  prop¬ 
erties  for  the  AS-3502  graphi te/epoxy  system. 

Results 

In  order  to  assess  reproduc i bi 1 i ty  and  data 
scatter  several  plates  were  cooled  down  along  the 
same  time-temperature  path.  We  therefore  had  9 
samples  cooled  along  path  "B",  2  samples  along 
path  "C"  and  3  samples  in  each  of  the  remaining 
cool-down  paths  "A",  "0",  and  "E",  for  a  total  of 
20  samples. 

Unfortunately,  the  relatively  thick  24  ply, 

6"  square  plate  presented  an  unforseen  problem. 

It  appears  that  the  thermal  excursion  from  cure  to 
room  temperature  gives  rise  to  stresses  that  are 
close  to  the  transverse  strength  of  the  laminas, 
thus  making  them  highly  susceptible  to  f-acture. 

In  fact  all  of  the  samples  except  those  cooled 
according  to  path  C  contained  cracks. 

Table  2  contains  the  measured  curvatures  as 
well  as  a  list  of  the  number  of  cracks  detected  in 
each  sample.  The  table  contains  also  the  theo¬ 
retical  curvature  computed  according  to  equation 
(12).  In  all  cases  the  cracks  occurred  along 
directions  paralled  to  the  fibers.  The  differ¬ 
ence  in  the  magnitude  of  K  and  k  can  be  attri- 
ed  to  the  uneven  number  ofxcracks-yi n  the  0°  and 
90°  plies.  Due  to  this  uneveness  it  was  found 
worthwhile  to  employ  the  average  of  K  and  K 
(K,  in  Table  2)  when  making  comparisons  to  tfle 
thioretlcal  curvature. 

It  may  be  noted  from  Table  2  that  the  effect 
of  the  cracks  is  to  diminish  the  magnitudes  of 
the  observed  curvature  to  values  lower  than  those 
predicted  by  equation  (12).  With  increasing 
nunber  of  cracks,  the  discrepancy  between  observed 
and  theoretical  curvature  is  found  to  increase. 

The  samples  containing  few  or  no  cracks  (i.e., 
samples  4,  13,  13  and  14)  had  average  curvatures 
which  were  very  close  to  those  predicted  by  (12). 


In  'able  2  there  seems  to  be  no  obvious  re- 
’ at  ion  between  the  cool  down  time  and  the  number 
jf  cracks  obtained.  This  is  most  apparent  by  ob¬ 
serving  that  the  samples  In  the  12  hour  coo)  down 
path  0  contain,  on  the  average,  more  cracks  than 
the  samples  In  the  4  hour  cool  down  path  B. 

It  is  quite  possible  that  the  cracks  were 
caused  by  the  fabrication  procedure.  When  a 
sample  is  layed  up,  a  cork  dam  is  placed  around 
it.  to  keep  the  resin  from  flowing  when  it  is  in  a 
nearly  liquid  state  uimediately  prior  to  cure. 

It  wa;  noticed  that  the  samples  in  cool  down  pro¬ 
files  A,  B,  0  and  E  cracked  when  this  cork  was 
removed  from  the  edge  of  the  samples.  Apparently 
the  samples  were  so  close  to  failure  that  even 
the  small  amount  of  pressure  applied  in  removing 
the  cork  sufficed  to  propagate  the  cracks.  To 
circumvent  this  factor  we  placed  a  teflon  sheet 
between  the  cork  and  the  sample  when  employing 
cool  down  profile  C,  which  allowed  these  plates  to 
be  effectively  separated  from  the  cork  dam  without 
the  formation  of  any  cracks. 

To  get  a  better  understanding  of  the  effect  of 
the  cracks  on  the  observed  curvature,  cracks  were 
purposely  induced  into  one  of  the  uncracked 
samples  (sample  14)  and  the  subsequent  change  in 
curvature  was  recorded.  Table  3  contains  the 
number  of  cracks  in  sample  14  with  the  resulting 
average  curvature.  Also  included  is  the  fraction 
of  the  uncracked  curvature  measured  orginally. 

The  number  of  cracks  in  sample  14  versus  the 
fraction  of  the  uncrqcked  curvature  resulting  from 
tnose  cracks  is  shown  in  Fig.  5. 

In  order  to  salvage  as  much  Information  as 
possible  from  the  flawed  specimens  we  employed  the 
results  depicted  in  Fig.  5  to  retrace  the  un¬ 
cracked  curvature  of  the  specimens  in  cool-down 
profiles  A,  B,  D  and  E.  These  retraced  uncracked 
curvatures  are  presented  in  Table  4,  where  the 
theoretical  values  predicted  by  equation  (12)  are 
listed  for  purpose  of  comparison. 


Table  2.  Experimental  and  Theoretical  Curvatures. 


Sample  Cool  Path  No.  Cracks 


K  K. 

Hfr  A 


On*'1)  (In’1)  (In*1) 


K  (Theory) 
(In*1) 


1 

A 

17 

.0218 

-.0158 

.0188 

.0249 

2 

A 

13 

.0186 

-.0178 

.0182 

.0252 

3 

A 

13 

.0190 

-.0218 

.0204 

.0253 

4 

B 

5 

.0211 

-.0271 

.0241 

.0250 

5 

B 

10 

.0194 

-.0236 

.0215 

.0241 

6 

S 

10 

.0209 

-.0213 

.0211 

.0243 

7 

B 

8 

.0200 

-.0178 

.0189 

.0241 

8 

8 

7 

.0247 

-.0164 

.0206 

.0234 

9 

B 

16 

.0221 

-.0138 

.0179 

.0243 

10 

8 

3 

.0278 

-.0207 

.0243 

.0247 

11 

B 

12 

.0204 

-.0218 

.0211 

.0240 

12 

B 

12 

.0151 

-.0247 

.0199 

.0246 

13 

C 

0 

.0274 

-.0226 

.0250 

.0252 

14 

C 

0 

.0272 

-.0226 

.0249 

.0256 

15 

0 

13 

.0175 

-.0208 

.0191 

.0247 

16 

0 

16 

.0183 

-.0215 

.01 99 

.0247 

17 

0 

30 

.0147 

-.0173 

.0160 

.0241 

18 

E 

8 

.0233 

-.0202 

.0217 

.0243 

19 

E 

7 

.0237 

-.0210 

.0223 

.0252 

20 

r 

9 

.0220 

-.0223 

.0221 

.0249 

\ 


Table 

3  . 

The 

Lfloct  ot  Cracks  n  the  Curvature  ot 

Specimen  14 . 

Number  ot 

Cracks 

(In  S  Fraction 

of  Uncrackf  K 

O 

.0249  1 

3 

.02)9 

.960 

5 

.0240 

.964 

7 

.0227 

.912 

s) 

•  J22  5 

.904 

10 

.020) 

.811 

12 

.0206 

.827 

14 

.0194 

.779 

16 

.0176 

.707 

18 

.0175 

.  "01 

•3  -> 

.0163 

.  bb3 

27 

.  01 3fc 

.627 

Table  4.  Retraced-Experimenta 1  Curvatures  and 
Theoretical  Values 


>pec imen 

Predicted  K  { In  * ) 

Theoretical  K  (In  ' 

— 

.0252 

.02.4 

2 

.  ii24>! 

.  12  52 

) 

.02-- 

.025  3 

u 

.  020° 

.02  50 

3 

•  J223 

.02-1 

r> 

.0243 

.02-3 

7 

.0248 

.0241 

3 

.0257 

.0234 

9 

.0242 

.0243 

10 

.0261 

.0247 

11 

.0228 

.0240 

1 2 

.0255 

.0246 

15 

.0239 

.0247 

16 

.0270 

.0247 

17 

.0255 

.0241 

18 

.0240 

.024  3 

19 

.0242 

.0252 

20 

.92-9 

.0249 

Numtui  >1  Cr»c k i 

f'aure  5.  Fraction  of  Uncracked  Curvature  Versus 
Number  of  Cracks  for  Sample  14. 

To  exclude  the  effects  of  variation  of  pla*e 
thickness  which  occur  from  sample  to  sample  we 
consider  the  dimensionless  quantity  hK  * 


*  Equation  (12)  shows  that  hKx  remains  constant  for 
each  specified  cooling  history. 


The  average  of  the  quantity  hKfl  for  each  cool  down 
path  is  listed  in  Table  5  together  with  the  corre¬ 
sponding  theoretical  value  of  hK. 

Very  good  agreement  exists  between  this  average 
observed  hK^  and  the  theoretically  predicted  value. 
Discrepancies  range  between  4%  for  cool  down  path 
D  and  less  than  IX  for  path  B.  Note  the  general 
agreement  with  viscoelastic  predictions,  as  compar¬ 
ed  with  up  to  12X  departure  from  the  linear-elastic 
resul t . 

The  optimal  time-temperature  path  was  calcu¬ 
lated  for  cooling  times  Lf  *  50,  100  and  200  min. 
using  equations  (20)  and  (21)  and  the  iterative 
scheme  outlined  in  (23). 

The  three  optimal  time-temperature  paths 
calculated  are  presented  in  Fig.  6  along  with  the 
re'ulting  time -dependent  curvature  predicted  by 
(12).  Discontinuities  occurred  at  t  *  0  and 
t  ■  tf,  with  a  nearly  linear  path  during 
tf  >  t  '  0. 

Figure  7  presents  the  optimal  curvature  versus 
the  logarithm  of  the  cool-down  time  tf  (tf  in 
minutes).  The  elastic  curvature  calculated  using 


Curvature  fMK#)  at 


1 


Table  5.  Comparison  between  Averaged  hK^  and  Theory. 


P  ich 

Measured  ind  Ad  lusted  *iK  (in/ln) 

Theoretical  hKx  fin/ In) 

Viscoelastic  E Las  tic 

A 

.00167  i 

.00004 

.00169 

.00186 

R 

.00162  * 

.noon 

.00163 

.00186 

.■HUM)  t 

.00002 

.00162 

.00196 

P 

.00167  t 

.00011 

.00160 

.00186 

E 

.00156  i 

. 00004 

.00159 

.00196 

(7),  is  shown  In  dashed  line.  As  expected  from 
viscoelasticity,  longer  cool-down  times  yield 
greater  stress  relaxation  with  a  subsequent  reduc¬ 
tion  of  the  optimal  curvature. 

Concluding  Remarks 

This  paper  presented  an  analysis  and  experi¬ 
mental  results  of  the  effects  of  viscoelastic 
response  of  composite  laminates  during  the  thermal 
cool-down  stage.  Although  we  must  view  the  present 
data  as  preliminary,  it  provides  encouraging 
evidence  that  the  above-mentioned  effects  are 
detectable  by  direct  measurements.  Somewnat 
inadvertently,  the  fracture  and  failure  which 
persistently  occurred  in  the  experiments.  Indicate 
that  the  residual  thermal  stresses  are  of  severe 
magnitudes  and  must  not  be  ignored  in  laminate 
designs. 

Further  studies,  based  upon  this  paper  are 
currently  In  progress. 
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A  RAPIDLY  CONVERGENT  SCHEME  TO  COMPUTE  MOISTURE  PROFILES 
IN  COMPOSITE  MATERIALS  UNDER  FLUCTUATING  AMBIENT  CONDITIONS 


by 


Y.  Weitsman* 


Abstract 


Tnis  paper  presents  a  highly  efficient  numerical  scheme  to  compute 
the  moisture  distribution  in  composite  materials  and  adhesive  joints 
under  time  varying  ambient  relative  humidities  and  temperatures .  The 
moisture  diffusion  is  assumed  to  follow  Fick's  laws.  It  is  shown  that  by 
appropriate  switching  among  the  various  forms  of  the  analytic  solutions, 
all  involving  infinite  series,  it  is  possible  to  attain  extremely  high 
accuracy  by  means  of  a  meagre  number  of  terms. 

An  example  is  provided  to  illustrate  the  method. 


♦Professor,  Mechanics  and  Material  Center 
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BASIC  CONSIDERATIONS 


Consider  a  moisture  sorption  process  that  is  described  by  the 
classical  diffusion  laws.  In  the  one  dimensional  case  we  have 

=  D  i!m  ( i ) 

?t  3x? 

to  which  we  must  attach  initial  and  boundary  conditions. 

In  (1)  m  =  m(x,t)  is  moisture  content,  x  is  the  spatial  coordinate, 
t  is  time  and  D  is  the  coefficient  of  moisture  diffusion. 

It  has  been  observed  ^  ^  that  the  equilibrium  moisture  content 
depends  on  the  ambient  relative  humidity,  and  we  shall  also  assume  that 
the  boundary  conditions  are  determined  by  the  same  quantity.  Further-, 
more,  the  moisture  diffusivity  was  found  to  be  most  sensitive  to  tempera¬ 
ture  ^  Several  empirical  relationships  were  proposed,  and  we 


shall  employ 

K  =  Cra 

(2) 

D(T)  =  Dp  exp(A/TR  -  A/T) 

(3) 

In  (2)  and  (3)  M  is  the  equilibrium  moisture  content,  r  the 
ambient,  relative  humidity,  T  the  temperature,  TR  the  reference  temperature 
and  A,  C,  a  the  material  constants. 

In  accordance  with  previous  analyses  ^  ^  we  can  uncouple  the 
process  of  heat  diffusion  from  all  other  time-dependent  material  processes, 
e.g.  moisture-diffusion  or  stress-relaxation.  This  simplification  is 
justified  because  for  all  practical  temperature  fluctuations  and  geomet''ical 
dimensions  the  time  required  to  reach  thermal  equilibrium  is  several  orders 
of  magnitude  shorter  than  the  time-scales  for  moisture  diffusion  or  for 
relaxation  response.  Consequently,  we  consider  spatially  uniform  temperature 
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profiles,  namely  T  =  T(t)  as  prescribed  by  the  fluctuations  in  ambient 
temperature,  when  analyzing  transient  moisture  diffusion. 


SYMMETRIC  EXPOSURE 

Consider  an  infinite  plate  of  thickness  21.  Let  -  L  <  x  <  L  and 
assume  an  initial  uniform  moisture  distribution  m  .  When  the  plate  is 
exposed  to  an  elevated  ambient  relative  humidity  the  boundary  moisture  is 
given  by  p,  namely  m(x  =  +  L,  t)  =  p.  Due  to  the  symmetry  of  the  present 
problem  it  suffices  to  analyze  only  the  region  0  <  x  <  L. 

For  constant  p  the  moisture  content  m(x,t)  is  given  by  well  known 
expressions  ^  Since  we  aim  at  extending  those  expressions  to  the 

case  of  fluctuating  p ( t )  and  temperature  T(t)  we  choose  to  represent  them 
in  the  following  form 


m(x,t)  -  I 0 ( x , t )  =  p  I(x,t) 


(4) 


The  functions  IQ(x,t)  and  I(x,t)  take  two  alternate  forms 


I0(x,t)  = 


In  (5) 


C(x,t) 

or  ,  I  ( x , t )  = 

1  -  E(x,t) 


1  -  C(x,t) 

or  (5) 

E(x,t) 


E 

n=l 


C(x,t)  =  -2}  (-l)n+1  cos(pnx/L)  exp(-pn  t*) 


and 
E  ( x ,  t ) 


C-) 

E 

n=l 


(-1) 


n+1 


erfc  I  — - -x— -  |  +  erfc 

2/P 


(6) 


2n  -  1  +  x/L 

2/F 


(7) 


with  -  (2n  -  1 ) -/ 2 

* 

and  t  -  Dt/L 


(8) 
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The  complementary  error  function  erfc(z)  decays  rapidly  with  z. 

r9i  .1 

Us  asymptotic  value  is  given  by  1  J  erfc  z'c(/iT z)  exp(-z  ) 

_ 1 6 

consequently,  for  computational  precision  of  0(10"  )  -  as  obtains  in 

douUe  precision"  routines  in  digital  computers  -  we  can  set 
erfc  z  0  for  <:  >  5.877.  In  the  sequal  we  shall  designate  this  number 


bv  • . 

The  rapid  decay  of  erfc  z  implies  that  series  (7)  converges 
rapidly  for  short  times.  On  the  other  hand  it  is  obvious  that  series 
(6)  converges  rapidly  for  long  times.  To  achieve  computational 
efficiency  we  should  therefore  switch  among  the  two  forms  of  equations 


(5) 


1  6 


Straightforward  arithmetics  yields  that  accuracy  of  0(10  )  is 


maintained  by  the  following  set  of  rules 

2  /  v  2 


for 


for 


for 


i  - 

Tr 


21 


<  t*  <  (  -?J  |  use  i  terms  in  series  (/) 


t  <  use  five  terms  in  series  16) 

9- 


(9a 

(9b 


- 2  <  t  <  — - ^ - 7  use  i  terms  in  series  (6)  (9c 

(21  ♦  1)  (2i  -  1) 


In  (9a)  and  (9c)  i  1,  2,  3,  4.  Also  Q  -  14.93  and  X  -  5.877. 

★ 

For  t  ■>  Q  the  moisture  distribution  is  uniform  to  within 

o(io'1" )• 

It  follows  from  (9a)  that  we  never  need  more  than  the  four  follow¬ 
ing  terms  in  series  (7): 


erfc  |  i-"— */k  )  f  erfc  (  -  -—-A  j  .  erfc  [  AA 

2/F  }  \  2/r*  \  2/A 


erfc 


It  can  be  noted  that  the  form  of  expressions  (9)  remains  valid  for 
any  desired  accuracy  e,  except  that  X  and  Q  depend  on  Obviously  for 


) 

) 


3  +  x  /l 
2vP 
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a  smaller  accuracy  we  require  even  fewer  terms  in  (6)  and  (7). 

Consider  now  the  case  of  fluctuating  temperatures,  T  =  T(t).  In 

view  of  (3)  the  diffusivity  D  is  now  time  dependent  and  the  non  dimensional 

★ 

time  t  in  (8)  becomes  a  complicated  function  of  real-time  t.  However, 

*  ? 

if  we  consider  t  =  DRt/L  at  the  reference  temperature  T  =  TR  then  in 

r  l  ol  * 

analogy  with  thermoviscoelasticity  we  can  replace  t  with  the  reduced 

★ 

non-dimensional  time  f  whenever  T  =  T(t)  as  follows 


L 


*  f 
J  0 


exp[A/TR  -  A/T(s)]ds 


(10) 


The  moisture  distribution  under  fluctuating  temperatures  is  given 

•k  ★ 

by  (4)  with  t  replaced  by  f,  .  In  view  of  the  si ngle-val uedness  of 
★  ★ 

(t)  it  is  always  possible  to  convert  the  results  back  to  real 

time  t. 

Consider  next  the  case  of  fluctuating  ambient  relative  humidity 
r  =  r(t) .  By  equation  (?)  this  implies  i.  .  (t ) .  Wnen  i.  -  ;.( t)  and  T 
equation  (4)  yield,,  upon  employment  of  the  superposi t inn  integral 


-t. 


m(x,t)  -  :n0 1 0 ( x ’ r  ) 


f  !!>,•  (t)  -  b  (:)]  (  )  d : 

0 

Equation  (11)  must  of  course  be  evaluated  numerically. 


(ID 


T(t) 


NON-SYKMETRIC  EXPOSURE 

Consider  now  an  infinite  plate  of  thickness  L  whose  faces  x:0  and  x=L 
are  exposed  to  different  relative  humidities  which  fluctuate  independently 
of  each  other.  We  still  assume  that  all  temperature  fluctuations  are 
spatially  uniform  within  the  entire  plate. 
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The  solution  to  differing,  but  constant  boundary  conditions 
m(o,t)  =  u°  and  m(L,t)  =  u*"  with  null  initial  moisture  m(x,o)  =  o  can 
be  expressed  as  follows 


m(x,t)  =  u°  H0(x,t*)  +  uL  H. ( x , t* ) 


(12) 


wi  th 


In  (13) 


S0(x,t  ) 


Hn(x,t  )  =  lor 


U0(x,t  ) 


HL(x,t  )  = (or 


SL(x,t  ) 


UL(x,t  ) 


(13) 


So(x,t  )  =  1 


SL(x,t 


|  y^n"sin~TL  exP(~n?1t?t  ) 
n=l 

oo 

=  r +  |^^xP(-n,'V) 


*  X/J  *  * 

U„(x,t  )  -  erfc  (■-—  )  +  Vn( x ,t  )-W„(x,t  ) 

?/r* 


★  1  w  /  1  *  * 

U.(x,t  )  =  erfc  (---—-)  +  V,  ( x , t  )  -  W.  (x,t  ) 
L  2.T*  L  L 


(14) 


(15) 


The  functions  V  and  W  in  (15)  represent  the  following  infinite 
series 


Vo(x,t*) 


erfc 


?n+x/l_\ 

/ 


W  o  (  x  ,  t  ) 
VL(x,t  ) 
WL(x,t*) 


E  -  (“) 


E 

n r  1 

E 

ri  -  0 


erfc  /2nil+x/L\ 

V  2  fir) 


(16) 
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Expressions  (14)  are  available  in  the  literature 


,  while  (15) 


[8] 

is  obtained  by  means  of  a  straightforward  Laplace  transform  and  inversion 
method. 

Maximal  efficiency  in  evaluating  H0  and  is  again  obtained  by 

switching  between  their  alternate  forms  given  in  (13)  and  detailed  in 

(14)  -  (16),  because  (14)  is  efficient  for  long  times  and  (15)  is  advantageous 

_  i  c 

for  short  times,  for  instance,  for  an  accuracy  of  0(10  )  we  never  need 

more  than  four  terms  in  each  of  U„,  U.  Sn  and  S.  as  listed  in  Table  1  below. 

L  >  L 

For  a  lesser  accuracy  the  number  of  terms  is  of  course  smaller. 


Range 

1 _ _ _ 

Largest  Number 

of  Terms  in  Each  Series 

Total  Number 
Terms  ' 
in  (1?) 

■9 

Wo 

WL 

So 

SL 

*  , 

O  t  <(t 

\ 

0 

•- 

1 

0 

0 

0 

0 

0 

2 

1  2  * 
(])  t  • 

-I  ? 

rJ  \ 

i 

0 

0 

0 

4 

1 

(?V>  t  ■ 

(;■) 

l 

1 

1 

1 

0 

6 

'  •  ★ 

s )  t  ' 

t 

2 

2 

1 

1 

0 

0 

8 

1 1 
r\ 

0 

0 

0 

0 

i 

1 

2i 

(i 

Ml 

i 

(i -1,3,3, 

1,0) 

_ 

Table  1:  Number  of  Terms  Required  in  Various 
Truncated  Series  to  Attain  Accuracy  of  0(10~lf>)  in 
Moisture  Profile  ( X-5 ,877  ,  R  =  16/r?  log  e) 

For  fluctuating  boundary  conditions  n°(t)  and  p^(t),  and  with  varying 

★ 

temperature  T(t)  we  employ  the  reduced  time  i,  given  in  (in)  and  a  superposition 
integral  analogous  to  (11)  to  get 


/I 

j  H  o  [  * ,  r(t)  - 

U  { 

+  nL[x,  "(t)-a i )]  ; l-  j  cii 

THE  NUMERICAL  SCHEME 

To  compute  the  moisture  m(x,t)  we  divide  the  time-span  of  interest  tf 
into  n,  not  necessarily  equal,  sub-interval  s.  These  intervals  ---  t^  -  t^ 

(i  =  1,  ...n)  with  t  =  0  and  t  =  t.  should  be  selected  in  a  manner 

x  9  ' 5  7  o  n  t 

that  both  the  ambient  moistures  |j°(t)  and  jiL  ( t )  as  well  as  the  temperature 

T(t)  are  represented  to  within  a  satisfactory  approximation  by  the  "stair- 

★ 

case"  functions 

t“(t!  =  T(t).T] 


for 


t  ■  ^  <  t  '  t  .  (i  =  1,  2,  ...n) 


Mote  that  in  the  symmetric  case  p  =  m  . 

n  O 

Denote  rj .  -  exp(A/T^  -  A/T.) 
then  (10)  yields 


-R-  E  <*. 


1,2,  ...n) 


k- 1 


Thereby 


.  .  -  r,. 

1 1  1 


V  E  ltk  -  lk-l>  3k 


k  jfl 


( j  -  0 ,  1 ,  ...  i - 1 ,  t  1  n 


(18) 


(19) 


(20) 


*  Obviously,  only  one  ambient  moisture  i ( t)  is  involved  in 
the  symmetric  case. 
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The  integrals  (11)  and  (17)  are  now  represented  respectively  by  the  sums 


mlX’t-j )  =  Vo(x’V  +  2  (“j 

j  =  l 


"j-l)  I(x’Sj) 


and 


m(x,t  • ) 


U].,)  Ho  ( x ) 


(21) 


(22) 


Expressions  (21)  and  (22)  remain  valid  for  any  interrediate  time 

t  ■  where  t.  ^  <  t..  <  t .  ,  provided  we  substitute  the  value  of  t.  in  place 

of  t.  in  (19)  and  (20)  as  well  as  in  (21)  and  (22). 

Computational  efficiency  is  achieved  by  switching  between  the  two 

alternate  forms  given  in  (5)  and  (13)  which  is  accomplished  by  testing  the 

ranges  of  f ^  and  r,.j  according  to  rules  (9a)  -  (9c)  or  in  Table  1,  respect- 
■*  *  * 

ively.  Obviously  p.  and  t..  must  replace  t  in  equations  (9)  and  in  Table  1. 

A  NUMERICAL  EXAMPLE!  FOR  THE  SYMMETRIC  CASE 

To  illustrate  the  method  we  consider  the  case  of  a  sixteen  ply 

5203/T300  graphite/epoxy  laminate  with  L  =  0.04".  For  this  material 

_  s  ; 

Dd  ;  1.5019  x  10  in  /'min  and  A  =  6340. 

K 

The  composite  laminate  was  considered  to  be  exposed  to  fluctuating 
ambient  relative  humidity,  which  is  reflected  as  a  fluctuating  boundary 
moisture  ;i,  and  to  fluctuating  temperatures. 

Two  cases  were  considered.  In  the  first  case  both  the  ambient  RH 
and  temperature  fluctuated  in  phase  while  in  the  second  situation  the 
fluctuation  was  out  of  phase.  Specifically,  in  both  cases  the  ambient 
moisture  content  fluctuated  between  1%  weight-gain  and  1/2:,  weight-gain 
every  5000  minutes.  In  case  1  the  temperature  varied  from  350°K  to  297°K 


every  5000  minutes  while  in  case  ?  the  temperature  fluctuated  between 
297 °K  and  350''K  with  the  same  frequency  of  5000  minutes  (but  out  of  phase 
with  the  moisture).  Case  1  is  shown  by  solid  lines  and  case  2  is  marked 
by  dashed  1 i nos  in  fig.  1  . 

The  results,  exhibited  in  Fig.  1.,  show  the  variation  of  moisture 
level  with  time  at  a  station  located  at  x  ~  0.035". 

Note  that  sharp  slopes  in  m(x,t)  vs.  t  occur  during  the  high-temperature 
time  intervals.  Consequently  the  in-phase  case  approaches  the  saturation 
level  of  m(x,t)  -  IT.  Conversely,  for  the  case  that  peak  levels  of 
T(t)  and  p(t)  are  out  of  phase  the  moisture  level  at  x  -  0.035"  approaches 
0.5*.  The  details  are  shown  by  the  heavy  lines  in  Fig.  1. 


Figure  Title: 

Moisture  Levels  at  x  5  0.035"  vs.  Time  in  a  0.08"  Thick  5208/T300 
Graph i te/Epoxy  Laminate  That  is  Fxposed  Symmetrically  to  Two  Cases  of  fluctuating 
Ambient  Relative-Humidity  and  Temperature. 

Case  1:  R.W.  In-Phase  With  Temperature 
Case  ?:  R.H.  Out-of-Phase  With  Temperature 
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Abstract 

Three-dimensional  constitutive  equations  are  derived  for  viscoelastic  compo¬ 
site  materials  with  time-dependent  damage.  Starting  with  the  assumption  that  the 
material  is  linearly  viscoelastic  when  the  damage  is  constant  and  that  the  damage 
consists  of  flaws  as  characterized  by  internal  displacement  discontinuities,  we 
investigate  the  effect  of  time-dependent  growth  and  healing  of  these  flaws  on  the 
global  stress-strain  equations  for  materials  under  transient  stresses  and  tempera¬ 
ture.  Prediction  of  flaw  growth  itself  is  then  discussed.  Some  concluding  remarks 
pertain  to  the  characterization  of  nonlinear  viscoelastic  materials  including  the 
effect  of  damage. 

Introduction 

The  the rmome chan i e a .1  response  of  many  materials  is  significantly  effected  by 
the  development  of  micro-flaws.  With  viscoelastic  materials  the  current  state  de¬ 
pends  on  the.  history  of  this  damage,  and  therefore  the  overall  material  response 
may  be  much  more  complex  than  for  elastic  materials  or  for  viscoelastic  materials 
without  damage.  Another  complication  results  when  the  previously  developed  flaw 
surfaces  rejoin  and  thereby  produce  an  increase  in  material  stiffness;  the  inter¬ 
face  contact  depends  on  the  history  of  flaw  development  and  global  thermomechanical 
inputs . 

Particles,  fibers,  and  lamina  interfaces  usually  serve  as  both  micro-flaw 
sources  (partly  as  a  result  of  local  stress  concentrations)  and  as  arrest  points 
after  local  stable  or  unstable  growth.  Therefore,  composite  materials  having  a 
brittle  or  weak  matrix  may  undergo  a  considerable  amount  of  global  softening  prior 
to  overall  fracture.  Considering  the  complex  microstructure  of  these  materials, 
one  may  expect  that  any  realistic  global  constitutive  equation  which  accounts  for 
damage  will  be  extremely  involved.  However,  as  shown  in  this  paper,  it  is  possible 
to  rigorously  develop  explicit,  realistic  equations  by  accounting  for  certain  sim- 
plif.'ing  features  of  many  composites,  such  as  a  matrix  that  is  relatively  soft  com¬ 
pared  to  the  reinforcing  material.  The  approach  followed  here  is  much  more  general 
ti  an  used  previously  for  a  particulate  composite  under  uniaxial  stress  (Schapery, 
ll)74a);  hut  lor  uniaxial  stress  and  without  rejoining  of  flaw  surfaces  the  present 
equations  take  the  same  form  as  the  earlier  ones. 

Inasmuch  as  the  purpose  of  this  paper  is  to  help  provide  a  basis  for  discussion 
at  the  Workshop  hv  describing  the  writer's  approach,  we  shall  not  attempt  here  to 
review  relevant  work  of  other  investigators.  Instead,  it  is  hoped  that  the  Workshop 
will  serve  this  purpose. 


♦Prepared  for  presentation  at  the  National  Science  Foundation  Damage  Workshop, 
Cincinnati,  May  4-7,  1980 


Linear  Viscoelastic  Constitutive  Equations  with  Damage 


Consider  an  isotropic  or  anisotropic  composite  material  element.  We  assume 
that  it  is  globally  homogeneous  with  or  without  damage;  i.e.  the  scale  of  stress/ 
strain  nonuniformit Les  due  to  any  physical  source  (cracks,  voids,  particles, 
fibers,  etc.)  is  assumed  small  compared  to  the  size  of  the  material  element.  We 
further  assume  that  the  material  is  linearly  viscoelastic  except  for  damage. 

Using  single  index  notation  for  the  global  or  average  stresses,  o^,  and 
strains,  c j ,  (Sokolnikof f ,  1956),  general  linear  relations  between  these  variables 
and  a  uniform  temperature  change,  AT,  may  be  written  in  the  form  of  nereditary 
integrals, 


ei  =  K/j  f  +  i,  j  =  1,...  6  (1) 

where  the  braces  are  abreviated  notation  for  a  hereditary  integral;  viz., 

{  Scr  |  H  f  S  (t  ,x)  |2-  dx  (2) 

o_ 

and  a  repeated  index  is  to  be  summed  over  its  range.  We  suppose  that  =  op  = 

AT  0  for  t<0;  the  lower  limit  in  Kq .  (2)  is  O'  rather  than  0  in  order  to  allow 
lor  step-function  inputs  at  t  =  0.  T.ach  quantity  Sjj  =  Sjj(t,r)  is  a  "creep  com¬ 
pliance,"  which  is  equal  to  the  strain  tj  at  time  t  due  to  a  unit  value  of  stress 
o  ^  applied  at  time  x.  Similarly,  the  thermal  coefficient  =  a^(t,r)  is  also  a 
creep  compliance  since  it  is  the  strain  c-  at  t  due  to  a  unit  value  of  AT  applied 
at  t.  If  all  of  these  compliances  are  functions  of  the  time  difference  t  -  t, 
instead  of  t  and  r  separately,  the  material  is  said  to  be  "nonaging;"  otherwise 
"aging"  exists.  This  aging  may  he  the  result  of  one  or  more  processes,  including 
chemical  reactions.  However,  we  assume  the  compliances  do  not  depend  on  stress, 
and  therefore  damage  ls  not  yet  included. 

The  damage  to  he  considered  here  consists  of  surfaces  of  internal  displace¬ 
ment  discontinuities  as  defined  by  a  set  of  discrete  parameters  um(m  =  7,  8,  ..., 

K) .  These  quantities  represent  the  three  components  of  relative  displacement 
(opening  and  sliding)  at  as  many  points  on  adjacent  internal  surfaces  as  needed  to 
accurately  represent  the  relative  movement  between  all  material  points  which  were 
together  before  the  damage  occurred,  all  surfaces  which  develop  during  the  time 
period  of  interest  are  to  be  included.  The  resultant  force  acting  between  a  pair 
of  these  points  prior  to  complete  separation  is  denoted  by  three  components  fm. 

This  procedure  oi  modeling  the  traction  and  relative  displacement  distribu¬ 
tions  along  current  and  future  surfaces  of  displacement  discontinuity  by  dis¬ 
placements  and  forces  at  discrete  material  points  is,  of  course,  what  one  would 
follow  in  a  finite  element  representation  of  the  entire  material  clement.  However, 
we  do  not  require  that  the  continuum  itself  be  represented  in  this  fashion. 

Rather,  it  is  assumed  only  that  the*  continuum  is  a  linear  viscoelastic  (aging  or 
nonaging)  material. 

There  is  assumed  to  be  a  thin  layer  along  the  surfaces  of  flaw  prolongation 
in  which  ail  large  strains  and  nonlinear  material  behavior  exists.  This  is  the 
same  assumption  used  by  Schaperv  (1975)  in  deriving  equations  for  predicting  speed 
o i  individual  cracks  in  linear  viscoelastic  media.  The  mechanical  and  failure 
behavior  of  this  thin  layer  is  defined  by  the  N-6  functionals, 

f  =  F (u  ,t);  m,n  =  7,  ...N  (3) 

m  n 

expressing  the  forces  between  adj  >  ent  material  points  at  the  continuum  boundary 
on  each  side  of  t  lie  thin  layer  of  nonlinear,  failing  material  as  functions  of 


displacement  history  and  possibly  time  (to  account  for  chemical  aging,  diffusion  of 
liquid  at  crack  faces,  etc.)*  In  predicting  global  mechanical  behavior  we  shall 
neglect  flaw-edge  details  and  simply  specify  that  either  a  force  is  zero  (e.g:.  a  free 
crack  surface) or  else  the  conjugate  relative  displacement  is  zero  (e.g.  a  point  in 
the  continuum  ahead  of  a  crack  tip).  Other,  more  general  cases  could  be  consider¬ 
ed,  including  friction  between  crack  faces,  pressure  of  fluid  in  cracks,  etc.;  but 
for  now  this  simplification  will  be  used.  Later  we  shall  briefly  consider  the 
effect  of  subsequent  contact  between  the  faces  and  healing.  It  is  proposed  to  use 
more  detailed  behavior  at  each  flaw-edge  as  given  by  Eq .  (3),  in  order  to  predict 
their  initiation  and  speed. 

It  should  be  emphasized  that  crack  growth  is  contained  in  the  following 
analysis  as  only  a  special  case,  and  not  the  only  case.  Local  failure  at  isolated 
individual  or  groups  of  points  is  taken  into  account.  In  principle,  therefore, 
phenomena  such  as  failure  of  individual  fibers  in  a  fibrous  composite,  groups  of 
extended  polymer  chains,  and  rubber  matrix  material  between  near-contact  points  of 
hard  particles  in  a  highly-f illed  rubber  can  be  treated.  For  example,  in  the 
latter  case  some  of  the  forces  fmwould  represent  forces  acting  between  particles  at 
the  points  of  near-contact,  and  the  rubber  between  any  pair  of  these  points  would 
be  characterized  by  one  of  the  functionals  in  F.q .  (3). 

Let  us  now  write  out  different  forms  of  the  linear  viscoelastic  constitutive 
equations  for  the  continuum  by  starting  with  the  generalized  form  of  Eq .  (1)  in 

which  global  strains  f:  and  internal  forces  tn  are  expressed  as  linear  functionals 
of  global  stresses  at,  relative  displacements  um,  and  temperature  change: 

f.  =  js.,6.}  +  is.  u  i  +  |  a.  AT  |  (4a) 

f  =  |  S  .6  .  |  +  |  S  u  |  )•  j  a  AT  }  (4b) 

m  1  mj  J  >  1  inn  n'  <  m  ’ 

Wien  all  un  arc  zero  the  material  is  by  definition  undamaged  and  Eq .  (1 )  is  re¬ 
covered  from  Eq .  (4a);  the  compliances  Sjj  and  ct  j  in  Eq .  (4a)  are  the  same  as 
those  in  Eq .  (1).  Another  useful  form  of  these  relations  is 


o,  =  |  C..  £.{+{(:.  n  \  -  {p.A'i'i 

i  1  ij  j  '  'inn'  *  l  ' 

f  =  ]e  .e .  Me  U  (  -  ]  B  AT  } 

m  '  mj  j  ’  '  ran  n  '  1  m  ’ 


where  Cjj  defines  the  functionals  which  arc  the  inverses  of  those  associated  with 
S  ,  ;  viz.. 


h  ’  lstj  tVt.ll 

reduce  to  relations  of  the  tvpe  in  En .  (24).  The  remaining  kernels  l'. 
obey  similar  relations;  e.g.  11 

-  |  S .  u  {  =  I  S .  .  |  C .  u  }  1 

1  inn'  '  l  j  '  jn  n '  ' 


The  functions  Cjj  define  mechanical  behavior  in  the  undamaged  state  (un  =  0).  They 
are  called  relaxation  moduli  as  they  are  the  stresses  due  to  unit  strains  applied 
at  t  =  t .  tn  general,  we  shall  refer  to  all  kernels  or  material  functions  in 
F.q.  (3)  as  relaxation  functions. 

Next,  we  shall  order  the  forces  fm  such  that  m  increases  with  increasing  time 
of  local  taiiure,  where  the  time  of  local  failure  is  defined  to  be  the  time  at 


which  fm  first  vanishes  when  u^  t  0.  Inasmuch  as  both  normal  and  shearing  forces 
are  included  in  the  sec  fm  and  failure  may  occur  at  more  than  one  point  at  a  given 
time,  the  number  of  discrete  failure  times  will  be  less  than  the  number  of  nodal 
points. 

With  this  in  mind,  suppose  that  tK(K  =1,  2,  ...)  is  the  Kth  distinct  failure 
time,  where  t-j  =  fg  =  ...f^  =  0,  and  all  material  points  corresponding  to  k  +  1 , 
k  +  2,  ...  N  have  not  yet  failed.  As  discussed  previously,  we  shall  predict  global 
response  by  neglecting  the  effect  of  a  partial  local  failure  (fm,  u^  +  0) . 

Write  out  Eq .  (4)  for  this  state  at  t  =  t^: 


where 


iVjl  + 1  sip"pt  +  IV*! 

j  S  .a . I  +  |  S  u  f  +  ja  AT { 
»  qj  j  '  <  qp  p>  1  q  1 

IS  .a  i  +  I  S  u  l  +  Ju  AT  [ 

r  I  r j  j *  1  rp  p 1  [  r  * 


'  1 

0  = 


p,  q=7,  8,  ...  k 

r  =  k  +  1 ,  k+2,  ...N 


(8) 

(9) 

(10) 


(11) 


Equation  (9)  provides  a  set  of  k  -  6  equations  from  which  the  k  -  6  values  of  Up 
may  be  found.  Substitution  of  these  displacements  into  Eq .  (8)  then  yields  the 
desired  global  constitutive  equations  for  the  damaged  material. 

In  order  to  obtain  some  explicit  results  we  shall  introduce  a  certain 
simplification  which  is  applicable  to  many  composite  and  monolithic  materials; 
this  point  will  be  argued  following  the  analysis.  Referring  to  Eq .  (5),  we  assume 
that  all  relaxation  functions  are  proportional  to  a  single  relaxation  modulus, 

E  -  E ( t , t ) (  for  uniaxial  loading  (say). 


nb 


=  C 


E_ 
ah  E,, 


B 

a 


(12) 


where  a,  b  ~  1,  2,  ...  N,  and  F.p  :-=  ER(tp,  Tr);  the  modulus  Er  is  a  constant  (the 
modulus  at  reference  times  t^,  r^)  which  is  introduced  so  that  the  dimensions  of 
tiie  constants  C°b  and  8°  will  lie  the  same  as  the  original  relaxation  functions. 


Suh^Cid  te  Eq .  (12) 

into  F.q . 

(5)  and  obtain  a 

set  of  equations  for  an 

"equivalent" 

elastic  material , 
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«  C°  .e‘?  + 

r  n  C 

C  .  U 

-  6°At? 

(13a) 
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After  replacing  all  strains,  displacements,  and  temperature  change  in 
Kqs .  (A)  and  (8-10)  by  the  corresponding  quantities  in  Eq .  (14)  only  constant  com¬ 
pliances  remain  (because  F.qs.  (4)  and  (8)— (If);  characterize  the  same  material 
element  as  I'.q.  (13)).  Obviously.  Eqs.  (R)-(10)  become 
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(15a) 

(15b) 

(15c) 


where 


[s°.]  =  fc!.l  1  ,-S!  =  S°.C!  ,  a!  =  S° . 6 ° 

L  ijJ  L  1JJ  ln  l3  Jn  1  iJ  J 

S°.  =  C°.S!  ,  S°  =  C°.S°  +  C°  ,  a°  =  C°.a°  - 

mj  mi  ij  mn  mi  in  mn  m  mi  l  m 

and  recall  that 

i,.j  -■=  1,  ...  6;  m, n  =  7,  ....  N;  p,q  =  7,  ....  k;  r=k  +  l,  N. 

The  relative  displacements  follow  from  Eq.  (15b)  and  the  definition  ["t°  1  =  1 


rT°  l=rso 
L  pqJ  L  pq. 


T°  S°  .o  .  -  T°  a°ATe 

pq  ij  j  pq  q 


and  substitution  into  Eq .  (15a)  yields  the  desired  constitutive  equations  at  the 


t  tmc  t  -  t^.: 


■  (sh +  asii)h +  f; +  4“i)“r‘  (i8> 

where  AS .  .  -  -  S°  T°  S°  .  ,  Act.  r  -  S!  T°  o°  (19) 

ij  ip  pq  q.)  l  ip  pq  q 

As  there  normally  will  he  a  large  number  of  points  of  local  failure,  it  is 
desirable  to  express  Fq .  (18)  in  terms  of  distribution  functions  for  AS-ji  .  Thus, 
let 

n,.(S,  tc)  dSdt  .  and  n.(S,  t.)  dSdt. 
ilt  1  it  1 

b  the  number  of  material  points  which  contribute  to  AS^.  and  Aa^  (respectively)  an 
amount  between  S  and  S  +  dS  when  the  local  failure  time  Is  between  tf  and  tf  +  dtf. 
Hence,  at  Lhe  current  time  t, 


,  /7\ 


.  .  (S ,  t  f ) SdSdt . 
ij  f  f 


-<*>  0 


An  .  = 


) SdSdt, 


(20b) 


Equal  Inti  (18)  can  now  be  written  ns 


e 
c  . 
1 


(21) 


=  (s!j  +/  Fij(tf)dtf)°j  +  (ftl+/  Fi(tf)dtf)^ 


where  new  distribution  functions  have  been  introduced, 

J  n  _  (  S ,  t  f ) SdS 


W  s 


K.(tf)  i-  f  n.(S,  tf)SdS 


6  jC- 

Kecall  that  (  .  and  AT  have  been  defined  in  Kq.  (14);  explicitly, 

i  Cl  3c- 

=  rl  T)  ardT 

R  o 


ATC 


■u 


a  at 

E(t,  t)  dr 


(22a) 


(22b) 


(23a) 

(23b) 


where  Hr  is  the  modulus  at  arbitrar i Iv  selected  reference  times  tR  and  tr.  It  is 
v.-rv  interesting  to  observe  th.it  the  constitutive  equations  (4)  and  (5)  for  a 
vis  oelastic  material  with  damage  have,  as  a  direct  result  of  Kq .  (12),  reduced  to 
those  for  a  linear  elastic  material  with  time-dependent  properties;  but  the  strain 
and  temperature  histories  are  in  general  different  from  the  actual  ones. 


1 1  t_ec  r  _of _ Flaw  (.'losing  and  Healing.  liquation  (21)  is  based  on  the  assumption  that 

the  tones  fm  are  gero  up  to  the  current  tin.-  after  local  failure  at  the  associated 
material  point  occurs.  However,  in  many  cases  this  will  not  be  true,  as  with  load 
rev rsn 1  and/or  the  application  of  a  large  external  pressure.  For  the  latter 
ns.',  inierfjc  ini  normal  and  shear  stresses  may  exist.  If  the  forces  f^  where 
this  occurs  are  known,  we  can  readily  correct  Kq .  (18)  tor  them.  Referring  to 
lq.  (18b),  we  see  that  they  will  be  accounted  for  if  the  associated  a"  in  Eq .  (19) 
are  replaced  by  -  (f^/AT*7).  On  the  other  hand,  the  flaw  closure  upon  unloading 
will  produce  contact  forces  which  depend  on  loading  and  temperature  history.  Also, 
the  (lavs  may  again  he  capable  of  supporting  tensile  forces  if  rebonding  occurs. 

In  this  latter  case  the  time  at  which  contact  is  established  may  be  appreciably 
if  f iv ted  i>v  the  i n t e rmo 1 ec u l a r  forces  of  attraction  if  they  are  sufficiently  large; 
indeed,  the  closing  rate  is  governed  by  an  equation  that  is  analogous  to  that  for 
e rack  growth  (Schaperv,  19/0 . 

In  order  to  deal  with  this  problem  involving  interfacial  contact  and  possible 
healing,  let  us  start  with  the  linear  constitutive  equations  in  the  form. 
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(24b) 


The  material  functions  are  different  from  those  in„Eq.  {it);  of  course,  they  can  be 
expressed  in  terms  of  the  latter  ones.  Note  that  and  ap  define  global  mechani¬ 

cal  response  for  the  maximum  amount  of  damage  (all  fm  =  0) ,  whereas  Sjj  and 
define  response  without  damage  (all  um  =  0) .  Let  us  again  introduce  the  simplifi¬ 
cation  given  in  Eq .  (12).  However,  rather  than  using  modified  strains  and 
temperature,  we  shall  apply  an  equivalent  modification  to  the  stresses  and  forces. 
First,  define  the  mechanical  creep  compliance  D  =  D(t,  t)  for  uniaxial  loading  by 

c  =  ( D( Ke> }  (25) 

Inasmuch  a»  1)  and  E  are  independent  of  strain  history,  we  can  use  in  Eq .  (25)  the 
unit-step  strain  history  f  =  H(t  -  t „ ) ,  where  H  (t  -  td)  vanishes  when  t  <  f0  and 
is  uni  tv  when  t  '  ic.  Equation  (25)  reduces  to  Lite  integral  equation, 


H ( t  -  f o)  = 


D(t,  f)  F. (t  ,  t  0 )  dx 


The  lower  limit  is  i  0  ,  rather  than  t0;  the  relaxation  modulus  E(x,  t0)  is  discon- 
t  inuous  at  t  =  t0,  and  therefore  the  contribution  from  the  singularity  in  3£/3t 
trust  be  included.  For  a  nonaging  material  L  =  E(t  -  t) ,  and  Eq .  (26)  may  be  easily 
Lap!.-..- t  tans  formed  , 

1  =  s  D  E  (27) 


whet  e  t  if  overbar  denotes  a  Laplace  transform  and  s  is  the  transform  parameter. 
Next ,  dot i ne 


^  l, 

I  "k  I  Jj  •>  "K  l  n  J 


UK  ^  FK1 


L.p.i.it  ior.  •  (5),  (12),  and  (28)  imply  l.q.  (7'.)  reduces  to  the  equations  for  an  equiva¬ 
lent.  clast  ic  material, 


r.  =  S'!  ,  c  .  +  S"  f  t-  ft®  AT 
i  i  i  j  inn  i 


u  S°  .O  +  S°  f  '  +  a°AT 
in  m  )  )  rnn  n  m 


wh.et  c  Si-  and  a®,  (with  a,  h  «  1,  .  .  .  ,  N)  are  constants, 

.  1  i»  <\  v 


T 1 1  i • .  choice  oi  variables  will  simplify  the  subsequent  contact  analysis. 
Howe.,  r,  it  is  first  of  interest  to  consider  an  implication  of  Eq .  (29a)  when  re¬ 
written  ia  terms  of  i  *[ ,  ATr ,  Oj  ,  and  f^.  For  a  state  of  constant  damage  (as 

d<.  lined  !■••  the  particular  set  of  f  which,  vanish)  Eqs.  (29a)  and  (21)  must  be  in 

agreement ,  regardless  of  whether  or  not  the  mechanical  variables  are  constant  in 

t ime.  lie;, 

1  * 

Si  ’  si,  w  Fn(S)d'f  •  a\  -“i  +j  Fi(tf)dtf 


a  . 
i 


(31) 


whore  t  is  the  time  at  which  the  constant  damage  state  is  first  reached. 

Now,  suppose  that  a  constant  damage  state  is  maintained  for  a  period  of  time 
which  is  long  enough  tor  the  values  of  f®  in  Eq .  (28)  to  essentially  vanish;  the 
stress  history  need  not  be  constant  during  this  time. 

The  quantity  {Df}  can  be  interpreted  as  the  mechanical  strain  due  to  a 
uniaxial  stress  history  o  -  f  (t)  for  an  undamaged  specimen;  therefore,  the 
behavior  of  { D f n }  when  fn  =  0  is  the  same  as  the  time-dependence  of  strain  in  a 
so-called  recovery  period.  Consequently,  it  is  necessary  to  assume  that  the 
material  element  under  consideration  is  such  that  the  strain  due  to  stress 
eventually  approaches  zero  after  stress  is  removed;  a  crosslinked  polymer  under 
count  ant  or  increasing  temperature  exhibits  such  behavior.  If  this  condition  is 
not  met,  then  the  long-time  value  of  f  (assuming  it  exists)  should  be  subtracted 
from  tf  and  the  difference  used  in  a  modified  version  of  the  subsequent  analysis. 

let  us  further  assume  that  um  =  0  at  all  points  where  contact  occurs;  i.e., 
material  points  on  adjacent  surfaces  are  assumed  to  rejoin  with  the  same  points  as 
before  local  failure.  In  this  case,  the  problem  is  completely  analogous  to  the 
problem  of  damage  growth,  in  which  the  roles  of  f®  and  if  are  interchanged  (cf. 
Eqs.  (lb)  and  (29).  Of  course  the  order  in  which  local  failure  occurs  is  not 
necessarily  the  same  as  that  tor  local  contact;  it  is  not  necessary  to  reorder  fn 
and  11^  to  reflect  this  behavior.  Equation  (29b)  can  be  solved  for  all  f^  where 


out. ut  occurs  since  th< 


are  zero.  If  desired,  one  could  then  solve  the  inte¬ 


gral  equation  in  Eq .  (281  for  the  contact  force  history,  fn  =  fn(t) ,  using  the 

'ii'uslv  derived  f'n 


u'i.ih  l  v  ut‘i  ivru  1  II  * 

In  view  of  the  ah 
iiuv the  same  form  as 


ove  considerations,  the  r 
t hose  in  fqs .  (2 1 ) - (2 3) : 


ting  constitutive  equations  will 


a +- 
l 


f  W' U,)‘T 

A  • 


(32) 


w  l  ( 


d  i  s t  r  ibut  i  ':a  1  unc  t  i  ons 
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n  .  .  ( S  ,  t 
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r.(t  ) 

1  c 


fn  t  ,s 

J  1  c 


1  SdS 


(33b) 


n  (S,  t  USdt  and  n  (S,  t  )dSdt 
II  c  c  i  c  i 


(39) 


t.  the  number  ol  material  point  s  wh  i<  !i  contribute  to  S  j  j  and  Acij  (respectively) 

:n  .r.  'writ  between  S  and  S  +  dS  when  tile  local  contact  time  is  between  t  and  t,.  + 

c  l 


Can  a  i  t  ut  ive  Kqqi.it  ion  ,  with  Cojnhjanqd  Damage  Growth  and  Flaw  C.ontac  t/Heal  ing  .  Both 
set  ■  at  equations  for  damage  growth,  Eq .  (21),  and  subsequent  rejoining  of  flaw 
sm  t.i.as,  Eq .  (321,  are  contained  in  the  following  set: 


In  the  development  of  the  theory,  tho  distribution  functions  in  Eq  ,  Ob)  were 
assumed  to  vanish  at  certain  times,  depending  on  whether  or  not  there  was  flaw 
sin  lace  contact  or  the  damage  was  constant.  However,  Eq .  05)  may  have  more 
general  nppl  ieahi 1  it  v.  11  internal  failures  and  contact  oo  nr  s imultaneously  at 
dill, Tent  points,  and  there  is  not  significant  intern,.  Lion  between  these  processes, 
tele  can  show  that  Eq  .  (V>)  will  be  appli,  ible. 

We  have  not  vet  discussed  the  mathemat ;.m1  proper; ics  of  tho  various  relaxa¬ 
tion  and  compliance  functions  wh  u  !i  appear  thru  ;  ghoul  *  tie  analysis.  Although  they 
will  not  be  listed  here,  it  is  of  interest  to  ohse: .-,  •'■at  since  the  material  is 
linearly  viscoelastic  for  a  fixed  amount  ot  damage,  ■  .  '  ibriura  and  nonequ  i  i  ihr  i  um 

t  lit1  rtm 'dynamics  provide  ,  onsiderible  expl  ic  it  inform',:  ;  ■  •  •  V.  •  fang, 

Si  haperv  ,  ldf,4)  .  for  ex  imp.  1  e  ,  S|-  is  symmetric  and  pes  i  t  ive  definite.  Morvver,  the 
a:  i  j.  (  15)  is  mi  -h  t  ha  t  material  syrm-trv  :  ements  may  be  easily  intro¬ 

duced.  Indeed,  restrictions  on  SO  and  i  j  cat:,  over  lircctly  from  linear 
clast  i,  ity  theory .  A 1  so ,  through,  the  integrals  ter  damage  and  contact  effects  in 
Kq .  (hit,  a  material's  symmetry  can  change •  For  example,  without  damage  the 
nutter  i  il  may  he  isot  route  (as  character  ir.ed  by  S”j  and  i  .  It  flaws  develop  into 
a  regular  pattern  (e. g. ,  an  array  of  parallel  cracks),  so  ,ts  to  produce  an  ortho¬ 
tropic  material,  say,  the  distribution  functions  in  Fq.  (in)  will  possess  the  same, 
t rans format  ion  properties  is  linear  elastic  moduli  of  orthotropic  materials. 

The  effect  of  temperature  on  mechanical  properties  is  included  in  the  present 
theorv.  For  example,  with  constant  or  transient  temperatures  the  aging  term  of  the 
various  material  functions  (u.g.  lift,  O),  permits  us  to  introduce  thermo rhee 1 ng i- 
c.i  1  1  .•  simple  behavior,  as  well  as  other  more  involved  behavior  (e.g.  Schapery , 

I  •»  '  .'•>'>  ;  this  is  ccomp  I  i  shed  through  an  appropriate  choice  ot  D(t,  T  )  . 

t  he  t  i  n.i  1  matter  in  this  section,  let  iu,  -ashler  the  basis  for  F.c  .  (Id). 

A;. -iime  tho  compos  i  t e  material  consists  of  u  .  c  a-  v  i  sooe  last  ic  phase.  Further  — 
pore,  siinpose  that  tliis  phase  i  -  isot  rope  o.  i  r.  is  a  t  irr.evise  constant  Poisson's 
ratio  uni  thermal  expinsion  coefficient,  wh : ,  ii  is  a  good  assumption  for  many 
material  (Sch.iperv.  lO/.'.h  )  .  Tin-  other  phases  are  o..s- timed  to  be  either  relatively 
vt,r v  r  ig  id  ot  s"tt  tholes  and  <  racks)  .  ihen ,  Fq.  (12)  can  h-  established  directly 
from  d  i  men.  i  on  a  i  analysis.  Ii  the  compos  j  t.  e  coil  line  continuous  and  straight, 
at  1 1  t  t  ii  .a  .,  it  is  not  necessary  to  i n *.  rod' ;c e  *  ne  severe  assumption  that  their 
.,  x  i  il  modulus  is  infinite.  Father,  the  tie..  in  Fq .  (It)  cm  he  redefined  Oj 
end  i  j.  nod  i  t  ied  a.  i  ord  ingl v )  so  that  it  does  not  include  the  axial  force  in  the 

t  i  I’ers  ;  (t  ii  •  t  :  t  ,f  broken  t  i :  r is  -.•.e.  i  1  v  p  ,■.!:■  t,\l  by  means  :  this  :  ■  ;  •  t  Lned 


s  t  r  e  s  ,  . 
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Stochastic  Models  of  Microcracking  and  Failure 

The  constitutive  relations,  Eq.  (35),  are  expressed  in  terras  of  distribution 
functions  for  the  time-dependent  number  of  local  failure  and  interfacial  contact 
sites.  These  functions  in  turn  depend  on  the  stress  and  temperature  histories,  and 
in  this  section  we  shall  briefly  review  the  writer's  approach  to  the  prediction  of 
this  dependence  as  well  as  the  prediction  of  global  failure. 


Analysis  of  Microcracking.  First,  it  is  noted  that  in  principle  the  failure  and 
contact  processes  can  be  predicted  from  Eqs.  (13)  for  the  continuum  (or  the  more 
general  set,  Eq.  (5)),  together  with  local  constitutive  equations  for  the  material 
at  the  failure  sites.  However,  considerable  simplification  is  needed  for  real 
materials  considering  their  highly  complex  microstructure.  We  treated  this  prob¬ 
lem  for  local  failure  in  past  work  (Schapery,  1974  a,  c)  by  applying  viscoelastic 
crack  growth  theory  to  isolated  preexisting  flaws;  based  on  experimental  data, 
crack  speed  was  assumed  to  obey  a  power  law  in  the  local  stress  intensity  factor 
for  the  opening  mode  of  growth.  Typically,  these  initial  flaws  are  predicted  to 
become  unstable  with  very  little  growth.  Therefore,  their  influence  on  global 
response  was  neglected  until  the  time  of  instability  (local  failure).  The  global 
softening  effect  was  taken  into  account  after  an  assumed  period  of  rapid  growth  and 
arrest.  The  analysis,  including  consideration  of  the  statistical  distribution  of 
initial  flaw  sizes,  microstructure  geometry,  etc.,  results  in  the  prediction  that 
global  softening  functions  (such  as  F. .  and  F  )  depend  on  only  the  current  Lebesgue 
norm  of  stress.  With  later  consideration  of  mixed-mode  growth  and  growth- 
retardation  effects  due  to  large  strains,  a  generalized  damage  parameter  was  pro¬ 
posed  for  isotropic  media  under  global  proportional  loading  (Schapery,  1978), 

t 

/q 

W(t’)— •  dt'  (37) 

e  " 
o 

where  q  and  p  are  positive  (and  typically  >>1).  W(t)  is  a  positive  function  of  time 
through  dependence  on  local  fracture  energy,  temperature-dependent  material  param¬ 
eters,  chemical  aging,  etc.  The  stress  a  and  strain  z  are  global  variables.  The 
use  of  total  strain  was  based  on  solid  propellant  data,  although  other  Integrand 
variables  (e.g.  strain  due  to  damage)  may  be  more  appropriate  for  different  materi¬ 
als;  but  use  of  a  function  of  L^.  itself  as  a  factor  in  the  integrand  has  no  essen¬ 
tial  effect  on  Eq.  (21)  becausefcq.  (37)  can  then  be  solved  explicitly  for  L^,  and 
the  resulting  effect  is  identical  to  that  without  the  factor. 

Considerable  simplification  results  if  the  distribution  functions  in  Eq.  (36) 
depend  on  time  through  a  single  parameter,  such  as  L, ;  of  course,  for  general  load¬ 
ing  conditions  one  would  have  to  use  at  least  suitably  defined  stress  and  strain 
invariants  instead  of  a  and  e.  In  this  event,  and  assuming  contact  effects  can  be 
characterized  using  an  analogous  parameter  L  (but  with  values  of  p  and  q  different 
from  those  in  Eq .  (37)) Eq.  (36)  becomes 

su  (t-  T)  5  Tu  |>(t)>  VT)] 

a*  (t,  1)  s  [lc<c),  Lf(x3  (38) 

where  T. .  and  T.  are  to  be  interpreted  as  material  functions  of  L^.  and  L  .  These 
function^  may  be  found  using  experimental  results  from  mechanical  tests;  Beckwith 
(1974)  demonstrates  this  in  which  a  one-dimensional,  isothermal  version  of  Eq.  (35) 


with  damage  only  (Schapery  1974a)  is  employed  using  results  from  multiple-step  creep 
and  recovery  tests. 

A  similar  one-dimensional  case  of  Eq.  (35).  together  with  a  damage  parameter 
similar  to  Eq.  (37),  was  recently  used  by  Schapery  (1979)  to  characterize  and  pre¬ 
dict  the  behavior  of  solid  propellant  with  non-decreasing  strain  input  under 
isothermal  and  nonisothermal  conditions.  Specifically,  the  stress  is 

o  -  °0/Td  (39) 

where  0  is  the  linear  viscoelastic  stress  a  =  (E(e  -  aAT) }  (without  damage),  W  is 

a  function  of  temperature,  and  is  a  funct?on  of  the  damage  parameter  Eq.  (37); 

the  strain  factor  was  not  restricted  to  a  power  law,  but  this  form  with  p^q^lO  fits 

the  data  quite  well.  Note  that  by  solving  Eq.  (39)  for  strain,  the  result  has  the 

same  form  as  Eq.  (35).  A  softening  function  T,  consistent  with  experimental  data  is 

a 

CL,/  q 

Td  =  e  (40) 

where  C  is  a  positive  constant.  Equations  (39)  and  (40)  may  be  used  to  obtain 
stress  as  an  explicit  function  of  strain. 


o 

where 

L  =f  W(t’)(ao/f)qdf  (42) 

Jo  ° 

in  which  a  >0,  and  f  =  f(e)  replaces  eP^q.  Equation  (41)  predicts  some  very  inter¬ 
esting  types  of  behavior  which  have  been  reported  for  solid  propellant .  For  example; 
with  q  >>  1,  L  >  0,  and  CL  >>  1,  the  stress  is  practically  independent  of  strain 
history.  Also°  given  the  strain  history  in  the  form  e  =  e.  g(t)  and  p  =  q,  we  find 
o/eA  is  independent  of  e  ,  but  the  material  is  nonlinear;  this  type  of  behavior  has 
been  reported  by  Farris  L1971)  for  solid  propellant  at  small  strains. 


=  o 


1  +  CL 

o 


-1/q 


(41) 


Stochastic  Model  for  Global  Fracture.  If  one  approximates  the  global  failure  of  a 
composite  or  monolithic  material  as  being  due  to  the  growth  of  one  dominant  flaw, 
then  a  relatively  simple  probabilistic  fracture  theory  results  by  using  the  same 
type  of  power  law  flaw  growth  model  used  to  develop  Eq.  (37).  The  principal  result 
is  expressed  by  the  equation  (Schapery,  1974c) 


Pf(0<t<tT)  =  PL  (43) 


where  P^  is  the  probability  of  failure,  P^  is  the  master  cumulative  distribution 

function  for  creep-rupture  tests  or  constant  amplitude  fatigue  tests  and  is 

the  final  time  of  interest.  Now,  with  the  condition  that  one  must  set  L  *=  L 

(=  largest  value  of  L  up  to  the  current  time)  whenever  the  following  equation3* 

predicts  L  <  L  ,  we  use  in  Eq .  (43)  the  expression 
r  max 


t 


L(t)  =  log 


(44) 


with  Lt  =  L(tT>;  also,  is  a  function  of  material  properties  and  m  and  k  are 
constants . 

This  result  can  be  viewed  as  an  extension  of  the  analysis  of  Halpin  et  al. 
(1973)  to  time-dependent  stresses  and  temperatures.  The  coefficient  c  is  not 
necessarily  constant;  for  example,  it  may  depend  on  strain,  damage,  tempera¬ 
ture,  and  complex  frequency  effects.  Special  cases  of  Eq.  (43)  have  met  with 
some  success  with  solid  propellant  and  fibrous  composites;  but  more  study  is 
needed  before  its  range  of  validity  can  be  established. 

Except  for  the  first  term  in  Eq.  (44),  which  is  often  negligible,  the  para¬ 
meters  for  global  fracture,  L,  and  for  local  damage,  L^,  are  essentially  the 
same.  This  is  a  direct  result  of  having  used  for  both  predictions  a  crack  speed 
equation  of  the  form  da/dt~K^,  where  is  the  opening-mode  stress  intensity 
factor.  Global  fracture  is  assumed  to  result  from  unstable  growth  of  a  dominant 
crack  in  the  opening  mode  of  deformation.  Therefore,  one  can  imagine  many 
situations  in  which  this  model  would  not  apply.  Nevertheless,  it  does  provide 
a  convenient  reference  case  against  which  more  involved  behavior  can  be  compared. 


Concluding  Remarks 

So  far  in  this  paper  we  have  considered  material  behavior  which  is  linearly 
viscoelastic  except  for  damage.  The  problem  of  developing  explicit  nonlinear 
viscoelastic  constitutive  equations  (with  constant  or  varying  damage)  which  are 
both  realistic  and  useful  is  of  course  much  more  difficult.  However,  a  study 
of  actual  behavior  does  reveal  certain  simplicity  which,  if  introduced  in  a  damage 
theory,  results  in  equations  which  are  not  much  more  complicated  than  Eq.  (21)  or 
(35).  Consider,  for  example,  the  behavior  of  carbon-black  filled  rubber.  As  a 
summary  of  the  findings  of  several  investigators  on  large  strain,  uniaxial  stress- 
strain  behavior  of  rubber,  we  may  write  (Mullins,  1969): 

C  =  F(amax)  g(0)  (*5> 

where  e  and  a  are  "engineering"  strain  and  stress,  respectively,  F  reflects  the 
effect  of  damage  in  that  it  is  a  function  of  the  maximum  value  of  stress,  o  , 
(considering  the  entire  history  of  loading)  and  g(o)  is  a  nonlinear  functior»aof 
stress;  viscoelastic  and  healing  effects  are  not  included  in  this  expression.  The 
coefficient  F  is  also  a  function  of  the  volume  fraction  of  particles,  but  their 
effect  on  g(a)  is  very  small. 

A  generalization  of  Eq .  (45)  that  contains  Eq.  (39)  for  viscoelastic 
behavior  as  a  special  case  may  be  obtained  by  replacing  o  in  Eq.  (39)  by  g(o). 

On  the  basis  of  this  finding,  one  is  lead  to  generalize  Eq.  (35)  by  proposing  the 
following  three-dimensional  constitutive  equation  for  anisotropic  or  isotropic 
materials , 


e  .  . 
-U 


(46) 


where  £..  and  a.,  are  suitably  defined  strain  and  stress  tensors  for  finite 
strain, x^nd  W  =1w(o..,  T,  damage)  is  the  Gibbs  free  energy  function  for  an 
elastic  material  wiit\  damage.  Similar,  slightly  more  involved  equations  have 
been  developed  in  part  from  nonequilibrium  thermodynamics  and  applied  to  solid 
propellant  (Schapery,  1973).  However,  their  range  of  validity  is  essentially 
unknown  as  experimental  data  are  still  very  limited.  It  would  certainly  be 
helpful  to  have  a  good  physical  model  which  both  predicts  the  type  of  simplicity 
exemplified  by  Eq.  (45)  and  is  valid  for  multiaxial  stress  states.  In  this 
regard,  it  is  of  interest  to  observe  that  we  can  derive  Eq .  (45)  for  nonlinear 
elastic  materials  by  assuming  the  flaws  do  not  interact  and  the  maximum  stress 
affects  their  growth  but  not  the  specific  process  involved  in  local  failure. 
Assumption  cf  a  distribution  of  local  strengths  or  initial  flaw  sizes  would  re¬ 
sult  in  the  dependence  on  maximum^tress ,  as  would  the  assumption  that  the 
coefficient  F  is  a  function  of  L,  ^  with  q  -*  «,  where  p  and  W(t')  are  finite 
(cf.  Eq.  (37)). 

Finally,  we  note  that  the  form  of  the  parameters  L,  in  Eq.  (37)  and  L  in 
Eq.  (44),  and  their  dependence  on  stress,  are  predicted  from  linear  viscoelastic 
fracture  mechanics  theory.  Dependence  of  the  integrand  on  global  strain  (and/or 
parameters  defining  the  global  damage  state)  would  be  predicted  from  the  linear 
theory  if  the  local  fracture  energy  for  the  material  at  a  crack  tip  depends  on 
these  quantities.  For  example,  with  solid  propellant  there  is  a  broad  distri¬ 
bution  of  particle  sizes,  and  the  material  at  one  crack  tip  may  contain  many 
much  smaller  flaws;  thus,  the  damage  state  of  the  composite  material  could  be 
expected  to  effect  the  fracture  energy  for  any  single  crack.  Consistent  with 
this  observation  is  the  fact  that  crack  speed  in  solid  propellant  becomes  inde¬ 
pendent  of  strain  when  the  strain  exceeds  the  value  at  which  new  vacoules  form 
(Schapery,  1979);  the  speed,  however,  continues  to  increase  rapidly  with  the 
stress  intensity  factor.  The  strain  may  affect  not  only  the  local  fracture 
properties,  but  also  the  nonlinear  form  of  the  energy  available  for  driving 
cracks.  This  latter  case  for  large  applied  strains  is  illustrated  by  Andrews  (1968) 
with  globally  elastic  materials  and  by  Brockway  and  Schapery  (1978)  with  visco¬ 
elastic  materials. 
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ABSTRACT 

A  method  using  a  longitudinally  split,  laminated  beam  specimen  is  developed 
to  obtain  delamination  fracture  toughness  as  a  function  of  the  rate  of  crack  pro¬ 
pagation.  First,  the  relation  between  energy  release  rate,  applied  displacement . 
and  various  laminate  parameters  is  derived  using  a  large  displacement,  small 
strain  theory.  Experiments  employing  glass/epoxy  composites  with  axially 
oriented  fibers  and  of  three  different  thicknesses  ai.d  a  wide  range  of  loading 
rates  are  then  described.  Although  the  beam  deflections  and  rotations  are  very 
large,  good  agreement  between  measured  and  predicted  beam  compliance  is 
demonstrated.  The  energy  release  rale  G  and  crack  speed  a  are  shown  to  obey  the 
power  law  G  •v  a";  essentially  (he  same  result  is  obtained  for  all  three  laminate 
thicknesses. 


INTRODUCTION 

Predictions  of  the  useful  life  for  composite  material  components  require  a 
comprehensive  understanding  of  the  material’s  response  to  complex  load 
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histories  in  anticipated  service  environments  with  the  presence  of  defects  such 
as  microcracks,  voids,  and  deiaminations.  Such  defects  can  occur  during 
manufacturing  or  may  develop  in  service,  causing  structural  degradation  or 
failures  at  stresses  well  below  the  strength  levels  expected  for  defect  free 
material.  Linear  elastic  fracture  mechanics  (LEFM)  has  been  developed  to  deal 
with  crack-like  defects  by  relating  defect  geometry  and  design  stress  to  a 
material  response,  normally  called  the  fracture  toughness.  The  fracture 
toughness  of  a  material  is  usually  characterized  by  critical  energy  release  rate 
Gc  or  the  critical  stress  intensity  factor  Kc. 

In  an  ideal,  monolithic,  isotropic  material  the  fracture  toughness  is  in¬ 
dependent  of  the  orientation  of  the  crack  plane  for  a  given  mode  of  deforma¬ 
tion  at  the  crack  tip,  such  as  the  opening  mode  (Mode  1).  However,  in  an 
anisotropic  composite  material,  the  material  response  may  vary  considerably 
depending  on  the  plane  of  fracture  and  the  energy  dissipative  processes  involv¬ 
ed.  Work  to  date  for  Mode  I  has  indicated  Klc  values  of  15-30  MPasA m  for 
center  notched  tensile  specimens  of  fiber-reinforced  plastic  laminates  [1 J;  these 
compare  favorably  with  aluminum  alloys  often  used  in  the  aircraft  industry 
which  have  K,c  values  of  23-44  MPa\/" m  12]. 

Growth  of  interlaminar  flaws  (delamination)  is  an  important  part  of  the 
failure  process  in  many  laminates  [3,4).  Compressive  fatigue  appears  to  be  an 
especially  severe  type  of  loading  in  producing  delaminations  [5,6];  out-of¬ 
plane  stresses  developed  through  compressive  loading  and  local  buckling 
are  thought  to  be  the  primary  cause  of  such  delamination  type  fractures.  These 
observations  indicate  that  the  delamination  fracture  toughness  may  be  the 
critical  toughness  parameter  for  fatigue  stressing  where  in-plane  stresses  are 
compressive. 

While  considerable  effort  has  been  expended  to  define  fracture  toughness 
for  tensile  loading  of  laminates  with  flaws  normal  to  lamina  planes,  very  little 
has  been  done  to  better  define  and  understand  the  delamination  fracture 
behavior.  Apparently,  the  few  studies  that  have  been  conducted  to 
characterize  delamination  have  utilized  surface  notched  specimens  (cf.  Figure 
la)  that  give  a  mixed  tensile  and  shear  stress  state  at  the  delamination  crack  tip 
[3,  7,  8];  precluded,  therefore,  is  determination  of  Ktt.  or  G,c  for  the  opening 
mode  of  delamination  in  which  only  tensile  stresses  across  the  crack  plane  exist 
near  the  crack  tip.  In  order  to  fully  characterize  delamination  fracture 
toughness  we  believe  it  is  necessary  to  study  the  effect  of  various  proportions 
of  tensile  and  shearing  stresses  at  the  crack  tip,  including  pure  tension. 

The  objective  of  the  investigation  described  herein  has  been  to  develop  an 
experimental  approach  with  the  associated  analysis  to  obtain  the  fracture 
toughness  for  the  opening  mode  of  delamination.  An  axially-split  beam 
geometry.  Figure  lb,  was  chosen  to  give  essentially  pure  opening  mode  frac¬ 
ture.  To  support  the  study  of  thin  laminates,  nonlinear  beam  theory  was  used 
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p 


Figure  I.  Specimens  used  in 
delamination  fractur ■■  toughness 
studies:  (a)  Wang/Mandell 
Specimen;  (b)  Split  laminate 
specimen  used  in  our  tests. 


in  the  analysis.  Experimental  measurements 
were  made  on  a  unidirectional  glass/epoxy 
composite  (Scotchply)  with  axially  oriented 
fibers  to  verify  the  analysis  and  allow  deter¬ 
mination  of  the  delammation  fracture 
toughness  of  this  material.  The  analytical 
approach  for  elastic  behavior  is  presented  in 
the  next  section,  followed  by  a  generaliza¬ 
tion  for  limited  viscoelastic  behavior  and  a 
description  of  the  experimental  program. 
The  experimental  results  are  then  used  with 
the  theory  to  characterize  analytically  the 
fracture  behavior  and  check  for  internal  con¬ 
sistency  of  the  results. 

The  thin,  split  beam  geometry  gives  rise 
to  stable  crack  growth,  and  therefore  is  par¬ 
ticularly  suited  to  determine  the  relation  bet¬ 
ween  slow  crack  speed  and  energy  release 
rate.  A  portion  of  our  study  therefore  has 


been  devoted  to  characterizing  the  fracture  toughness  of  a  viscoelastic 


laminate. 


ANALYTICAL  PROCEDURES  FOR  ELASTIC  BEHAVIOR 

Consider  the  beam  specimen  in  Figure  lb,  in  which  the  delamination  crack 
tip  is  at  the  point  x  =  a.  The  energy  release  rate  associated  with  a  virtual  crack 
growth  of  an  amount  da  is,  by  definition,  the  mechanical  energy  that  becomes 
available  at  the  crack  tip  per  unit  area  of  new  surface.  In  terms  of  the  total 
strain  energy  in  an  elastic  beam,  W,  this  release  rate  is  [91, 


/  3  W 

B  3  a 


(1) 


where  the  derivative  is  evaluated  for  constant  beam  tip  displacement  A.  Also, 
B  is  specimen  width  normal  to  the  page.  The  value  of  C  at  which  the  crack  ac¬ 
tually  starts  to  propogate  is  the  critical  energy  release  rate  or  fracture 
toughness,  Gc. 

Evaluation  of  Equation  (1)  for  our  specimen  geometry  is  based  on  an  ap¬ 
proximate  analysis  in  which  strain  energy  is  determined  for  the  cantilevered 
beam  in  Figure  2;  the  relevant  parameters  are  tip  displacement,  A,  beam 
length,  L,  area  moment  of  inertia,  /,  and  the  axial  modulus  of  elasticity,  E. 
The  analysis  for  a  nonlinear  beam  in  terms  of  these  parameters  is  presented 
next.  Primary  results  of  the  theory  are  given  in  dimensionless  form  in  a  table 
for  general  use  in  reducing  data. 
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Strain  Energy  in  a  Nonlinear  Beam 

The  strain  energy  due  to  bending  of  a  beam  of  length 
L  is,  according  to  elementary  beam  theory, 


W,  =  Vl 


L 

f 


M 1 

~eT 


ds 


(2) 


Figure  2.  Nomenclature 
used  in  nonlinear 
analysis  of  cantilevered 
beam. 


where  M  is  the  local  bending  moment. 

As  a  result  of  the  low  flexural  rigidity  of  thin  com¬ 
posite  specimens,  large  deflections  and  rotations  may 
be  present  during  the  test  and  need  to  be  accounted 
for  in  the  analysis.  Such  a  correction  to  linear  beam  theory  has  been  made  by 
Bisshopp  and  Drucker  [10],  in  which  they  allowed  for  arbitrarily  large  rota¬ 
tions  and  bending  deflections.  A  linear  stress-strain  equation  and  small  strains 
were  assumed. 

The  beam  nomenclature  and  geometric  relationships  are  shown  in  Figure  2, 
and  the  principal  results  are  [10]: 


fsin$o  —  sin4>)^  d§ 


o 


(3a) 


and 


A_ 

T 


1  \E1  ]  '  C  sin|  d$ 

sf2  [  PL 2  J  J  (sin$0  —  sin^)'* 

O 


(3b) 


where  is  the  angle  of  the  tangent  at  the  loaded  end  (cf.  Figure  2),  $  is  the 
angle  at  intermediate  points,  P  is  the  load,  and  the  remaining  terms  were  de¬ 
fined  previously.  After  determining  the  appropriate  change  of  variables  and 
transformations,  Bisshopp  and  Drucker  were  able  to  rearrange  Equation  (3) 
into  the  form 
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M  »  fw  -  m  e.) 


and 


A_=  !  _2  [  E(k)  -  E{k,  9,)  T 

L  l  F(k)-  F(k.  6 ,)  J 


(4a) 


(4b) 


where  F(Ar)  and  £(&)  are  the  complete  elliptic  integrals  of  the  first  and  second 
kind,  respectively;  also,  ©i)  and  E(k,  0,)  are  the  corresponding  incomplete 
elliptic  integrals.  The  elliptic  parameters  k  and  6,  are  related  to  <j>0: 


k  a  _L_(l+  sin^o)'"5 ,  0.  ®  sin‘‘  (\f2  k)  '  (5) 

v  2 

Notice  that  Equation  (5)  implies  the  results  in  Equation  (4)  are  functions  of  on¬ 
ly  one  parameter,  say  $0;  hence,  they  are  implicitly  related  with  the  cor¬ 
respondence  presented  in  the  second  and  third  columns  in  Table  1  and  in 
Figure  3. 

Table  I.  Nonlinear  Beam  Variables. 


♦o 

PL* 

t 

WL 

GBl* 

GB 

(df?9> 

FT 

L 

7n 

281 

2P 

Linear  ( 
Theory  j 

4 

L 

3,  4  •} 

IT 

1(X)2 
2V  l  ' 

f<f) 

10 

.3530 

•  1160 

.0203 

.0613 

.1737 

20 

.7306 

.2302 

.0817 

.2499 

.3421 

30 

1.1626 

.3406 

.1854 

.5814 

.5000 

40 

1.6923 

.4455 

.3338 

1.0880 

.6429 

50 

2.  3922 

.5437 

.5320 

1.6330 

.7662 

60 

3.4054 

.6340 

.7901 

2.9502 

.8663 

70 

5.0812 

.7167 

1.1331 

4.7758 

.9399 

80 

8.6787 

.7948 

1.6486 

8.5558 

.9858 

88 

14.9058 

.8472 

2.2411 

14.9293 

1.0002 

88 

20.7282 

.8711 

2.2600 

20.9704 

1.0117 

90 

- 

1.0000 

- 
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Next,  the  strain  energy  stored  in  both  cantilevers  comprising  the  split  beam 
is  calculated.  Combining  the  general  definition  for  strain  energy  W,  in  a  beam, 
Equation  (2),  with  the  specific  results  from  nonlinear  beam  analysis  given  in 
[10],  one  may  express  the  total  strain  energy  W  =  2W,  in  terms  of  elliptic  in¬ 
tegrals  as  follows: 


WL 

2EI 


=  21  F[k)  ~  F(k.9,)\' 


E(k )  -  £(k,e,) 


<*'  -  .)] 


(6) 


(6) 


L  m  —  F(k.e ,) 

This  nondimensional  result  is  presented  in 
Figure  4  and  the  fourth  column  in  Table  1. 
For  comparison,  we  record  also  the  linear 
solutions. 


PL1 

El 


=  3 


A 
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WL 
2  El 


2  L 


V) 


G  s— 


i  aw 

TUT 


wl 

2  El 


=  JI&/L) 


Energy  Release  Rate  G  for  a  Nonlinear  Beam 

The  energy  release  rate  as  defined  in  Equation 
(1)  can  now  be  evaluated  using  Equation  (6). 
However,  in  keeping  with  beam  notation,  the 
length  L  is  Used  instead  of  the  symbol  ,,a”  to 
denote  crack  length;  hence 


Figure  3.  Nondimensional  load- 
deflection  curve  for  cantilevered 
beam.  Symbols  indicate  ex¬ 
perimental  data  from  split 
laminate. 


(8) 


Prediction  of  G  may  be  easily  accomplished 
by  first  writing, 


(9) 


Figure  4.  Nondimensional  strain 
energy  for  cantilevered  beam. 


where  /  =  is  the  solid  curve  in  Figure  -1.  Thus, 


1  dW 

JUT 

Next,  define  S, 


El  y  .  A  +  _A_  df{fS/L) 
L‘~  (L  L  3(A/L) 


(10) 
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S-d  {log  f)/ d  (log  (A /£)) 
Therefore,  from  Equations  (8)  —  (II), 


G  = 


El  WL 


BL'  El 


0  +  S)  = 


W 

~BL 


(1  +  S) 


(11) 


(12) 


A  rearrangement  of  Equation  (12)  yields  a  nondimensional  form  of  the  energy 
release  rate, 


GBL1 
IE l 


WL 

1EI 


(1  +  S) 


(13) 


Since  both  WL/EI  and  S  in  Equation  (13)  are  functions  of  A /L,  GBL'/EI  is 
an  implicit  function  of  A /L;  it  has  been  evaluated  and  is  presented  in  Table  1 
and  in  Figure  5.  Linear  theory  yields  S  =  2,  and  therefore  from  Equations  (7) 
and  (13). 


GBL J  _  9  A 
2  El  2  '  L  ’ 


(14) 


which  is  shown  in  Figure  5.  The  last  column  in 
Table  1  is  the  ratio  of  the  fifth  to  the  second 
column;  this  ratio,  GB/2P,  enables  the  calcula¬ 
tion  of  G  without  using  the  flexural  rigidity,  El. 

Suppose  that  during  quasi-static  crack 
growth,  one  determines  A  and  L  at  any  given 
time  in  a  delamination  fracture  test.  One  may 
then  use  Figure  5  or  Table  1  to  determine  Gc 
directly  because  Gc  =  G  for  a  slowly  growing 
crack.  That  slow,  controlled  crack  growth  can 
indeed  be  achieved  with  this  test  is  discussed 
next. 

Figure  S.  Nondimensional  energy 
Stability  of  the  Crack  Growth  re!ease  ra"  f°r  ^milevered  beam. 


It  is  clear  that  increasing  the  load  or  grip  displacement  for  a  fixed  crack 
length  will  increase  the  strain  energy  stored  in  the  specimen.  From  Figure  4  it  is 
also  apparent  that  crack  extension  (increasing  L)  which  occurs  at  constant  grip 
displacement  results  in  a  decrease  in  the  strain  energy  stored  in  the  specimen, 
as  shown  schematically  in  Figure  6.  As  the  specimen  with  fixed  crack  size  L,  is 
deformed  from  A,  to  AJf  the  strain  energy  obviously  increases.  The  strain 
energy  release  rate  also  increases  as  the  test  goes  from  a  to  b  (cf.  Table  1). 
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Suppose  at  b  the  critical  energy  release  rate  is  ex¬ 
ceeded  slightly  and  crack  extension  for  constant  A 
occurs,  changing  the  crack  length  from  L,  to  L,. 

At  point  c  the  energy  release  rate  is  again  insuffi¬ 
cient  to  give  crack  extension,  so  that  further 
loading  to  point  d  is  necessary,  which  again  cor¬ 
responds  to  G  being  slightly  larger  than  the  critical 
energy  release  rate,  Gc.  Crack  extension  occurs 
from  L j  to  L,  as  one  moves  along  d-e.  We  may 
conclude  from  these  considerations  that  crack 
growth  is  stable  in  a  controlled  displacement  test  if 
Gc  does  not  change  with  L  (or,  at  least,  if  it  does 
not  decrease  rapidly  with  L). 

In  an  actual  test  the  loading  and  crack  exten¬ 
sion  occur  more  or  less  continuously,  instead  of  in 
steps.  Indeed,  with  negligible  kinetic  energy 
(quasi-static  loading  and  slow  crack  growth)  the 
energy  release  rate  is  essentially  equal  to  Gc  at  all  times.  The  path  a-c-e  for  G  = 
Gc  represented  by  the  dotted  line  (not  necessarily  straight)  in  Figure  6  would  be 
followed  in  an  actual  experiment. 

In  summary,  the  relationships  between  A/L  and  P,  W  and  G  have  been 
determined  in  terms  of  elliptic  integrals  using  nonlinear  beam  analysis.  These 
results  are  summarized  in  Figures  3,  4  and  5  and  Table  1;  they  will  be  used  to 
analyze  the  experimental  results  following  the  next  section,  which  is  concerned 
with  extension  of  these  results  to  include  rate  effects  at  the  crack  tip. 

EFFECT  OF  LOCAL  VISCOELASTIC  BEHAVIOR 

The  above  theoretical  prediction  of  energy  release  rate,  G,  can  be  used  for 
viscoelastic  materials  if  the  axial  stress-strain  relation  is  essentially  elastic  and 
viscoelastic  effects  arc  limited  to  a  small  zone  around  the  crack  tip.  In  many 
cases,  these  conditions  will  be  met  for  thin  laminates  consisting  of  a  polymeric 
matrix  and  a  large  volume  fraction  of  continuous  glass,  graphite,  or  boron 
fibers  oriented  in  the  axial  direction.  If  the  matrix  is  very  soft,  large  inelastic 
shear  deformations  and  possibly  significant  microcracking  may  occur  away 
from  the  crack  plane.  Assuming  these  complications  do  not  exist,  the  material 
in  the  neighborhood  of  the  crack  tip  must  absorb  essentially  all  of  the  available 
energy  during  slow  crack  growth,  and  we  may  write 

G  =  2F,  (15) 

The  quantity  is  the  fracture  energy;  it  is  the  energy  required  to  produce  a 
unit  of  new  surface  at  the  crack  speed  L.  This  energy  may  include  energy  of 


Figure  6.  Schematic  of  strain 
energy  versus  crack  length  for 
constant  grip  displacement. 
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polymer  chain  failure,  energy  absorbed  in  local  (near-crack  tip)  viscoelastic 
deformation  processes,  fiber  fracture,  and  fiber-matrix  debonding.  In  analogy 
with  the  discussion  in  [12]  for  monolithic  materials,  we  conclude  if  the  speed 
does  not  change  appreciably  in  the  time  required  for  the  crack  to  propagate  a 
distance  equal  to  the  scale  of  the  local  deformation  and  failure  processes, 
will  depend  on  the  instantaneous  crack  speed  but  not  on  the  history  of  loading 
and  not  on  the  rate  of  change  of  L.  If  the  speed  dependence  of  is  due 
primarily  to  viscoelastic  behavior  of  the  intact  matrix  near  the  tip,  we  further 
expect  it  to  obey  a  power  law,  <t>b^  L",  where  n  is  constant  [12).  Thus,  in  view 
of  Equation  (15),  the  energy  release  rate  will  obey  the  same  power  law,  C  ~ 
Ln. 

Whether  or  not  the  result  for  the  split  beam  (viz.,  G  =  G(L))  can  be  used 
with  other  geometries,  such  as  growth  of  a  through-crack  in  a  tensile  sample 
with  fibers  at  90°  to  the  load,  depends  at  least  on  the  scale  of  the  zone  of  frac¬ 
ture  at  the  crack  tip  and  the  extent  of  viscoelastic  behavior.  It  is  expected  that 
it  will  apply  as  long  as  viscoelastic  and  failure  processes  are  highly  localized  to 
the  crack  tip.  With  large  scale  viscoelastic  effects,  the  value  of  T6  may  not 
change,  but  one  cannot  use  Equation  (15)  as  it  relates  f6  to  an  elastic  energy 
release  rate. 

The  effect  of  crack  speed  is  illustrated  schematically  in  Figure  6,  where  dif¬ 
ferent  paths  are  followed  for  different  crack  speeds.  For  the  usual  case  in  which 
rffy'dLX),  the  path  will  move  upward,  as  shown  in  the  figure. 

Finally,  it  is  to  be  noted  that  if  is  constant.  Equation  (15)  may  be  replac¬ 
ed  by  the  equivalent  statement  G  =  which  is  the  critical  energy  release  rate. 
We  prefer  to  not  use  the  symbol  Gc  instead  of  2rfc  for  materials  in  which  rfc 
depends  on  crack  speed  because  the  term  “critical  energy  release  rate”  normal¬ 
ly  refers  to  a  quantity  that  defines  the  boundary  between  no-growth  and 
growth. 

EXPERIMENTAL  PROCEDURES  AND  RESULTS 

The  experimental  program  to  be  described  in  this  section  had  two  objec¬ 
tives:  (1)  to  experimentally  verify  the  nonlinear  beam  analysis;  and  (2)  to  ob¬ 
tain  the  relation  between  energy  release  rate  and  crack  growth  rate.  First, 
specimen  preparation  and  testing  will  be  described  and  typical  raw  data 
presented.  Reduction  of  the  data  to  energy  release  rate  values  using  the 
previously  described  nonlinear  beam  analysis  is  then  accomplished,  and  the 
results  are  presented  in  tabular  and  graphical  form. 

Specimen  Preparation 

Scotchply,  which  is  an  E-glass  reinforced  type  1003  epoxy,  was  selected  for 
this  study  because  it  is  translucent,  which  makes  the  visual  measurement  of  the 
moving  crack  front  during  testing  quite  easy.  Eight,  twelve,  and  sixteen  ply 
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30.5  cm  x  30.5  cm  panels  were  laid  up  using  prepreg  tape  containing  con¬ 
tinuous  fiberglass  filaments.  Teflon  was  inserted  along  one  edge  of  each  panel 
to  provide  the  initial  midplane  delamination  crack.  All  lamina  had  the  same 
fiber  orientation  in  the  laminates.  After  curing  in  accordance  with  the  sup¬ 
pliers  specifications,  the  30.5  cm  square  laminates  were  cut  parallel  to  the  fiber 
direction  into  2.54  cm  wide  unidirectional  test  strips.  The  cured  lamina 
thickness  was  about  0.021  cm,  giving  laminate  thicknesses  of  0.170,  0.259  and 
0.338  cm,  respectively,  for  the  8,  12  and  16  ply  specimens. 


Test  Procedure 

The  specimens  were  tested  in  an  ambient  en¬ 
vironment  (approximately  75  °F  and  50%  RH)  us¬ 
ing  an  Instron  tensile  test  machine  with  a  one 
thousand  pound  load  cell;  a  twenty  pound  full 
scale  range  was  used.  Special  grips  were  designed 
to  allow  the  toad  to  be  applied  along  a  fixed  line  of 
action  while  permitting  a  virtually  free  rotation  of 
the  loaded  ends  of  the  split  beam  (cf.  Figure  7),  as 
was  assumed  in  the  analysis.  Load  was  monitored 
as  a  function  of  time  using  a  strip  chart  recorder. 
The  displacement  versus  time  was  calculated  from 
the  crosshead  speed.  The  crack  position  as  a  func¬ 
tion  of  time  was  noted  visually  as  the  crack  front 
passed  reference  marks  on  the  specimens.  Typical 
experimental  results  are  presented  in  Table  2  for  a 
twelve  .  y  specimen  tested  at  a  crosshead  speed  of 
2.5  cm/minute.  Five  each  of  the  eight,  twelve  and 


Figure  7.  Split  laminate 
specimen  showing  (a)  hinge  at¬ 
tachment  used  fur  loading  and 
lb)  deformation  in  crack 
growth  test. 


sixteen  ply  specimens  were  tested  at  a  crosshead  speed  of  2.5  cm/minute.  Two 


specimens  were  tested  with  increasing  crosshead  speed. 


Data  Reduction 


The  measured  data  included  load,  P,  grip  displacement,  2A,  and  crack 
length,  L;  they  can  be  easily  used  to  give  instantaneous  energy  release  rate,  C, 
using  the  analysis  described  in  an  earlier  section,  especially  the  results  in  the 
last  column  in  Table  1.  After  determining  A /L,  Lagrangian  interpolation  was 
used  to  calculate  both  Pl.'/EI  and  GB/2P.  The  results  of  these  interpolations, 
along  with  the  subsequent  calculation  of  the  energy  release  rate,  G,  are  sum¬ 
marized  in  Table  2.  The  crack  gro-'th  rate  was  calculated  from  the  measured 
values  of  grip  displacement  versus  crack  extension,  L,  knowing  the  rate  of  grip 
displacement.  From  each  set  of  measured  values  of  P,  2A,  and  L,  the  fifth  col¬ 
umn  in  Table  2  was  used  to  derive  individual  values  of  F.I\  they  were  averaged 
with  respect  to  length  for  each  laminate  for  later  use. 
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Table  2.  Experimentally  Measured  Values  of  Crack  Length.  L,  Load,  P,  and  Grip  Displacement. 
2A,  for  Twelve-Ply  Specimen.  Also  Tabulated  are  Various  Dimensionless  Terms  Used  in  Deter¬ 
mining  G  and  Checking  Theory. 


L 

(cm) 

2e 

(cm) 

P 

(N) 

4/L 

PL2 

FT 

(Theory) 

G8 

(Theory) 

G 

N/m 

fi.  9 

5.6 

27.0 

.32 

1.07 

.470 

991 

to.? 

7.2 

24.4 

.36 

1.26 

.530 

1010 

It.  4 

9.2 

22.1 

.40 

1.47 

.583 

1020 

12.7 

11.1 

20.6 

.44 

1.66 

.636 

1020 

14.0 

13.1 

IB. 9 

.47 

1.83 

.676 

993 

15.2 

15.1 

17.8 

.50 

2.04 

.713 

987 

16.5 

17.6 

17.0 

.53 

2.23 

.750 

1010 

17.8 

19.8 

16.3 

.56 

2.54 

.786 

1010 

19.1 

22.4 

15.8 

.59 

2.85 

.819 

1020 

20.3 

24.7 

14.8 

.61 

3.10 

.843 

978 

21.6 

26.9 

14.1 

.62 

3.22 

.853 

949 

22.9 

28.9 

12.9 

.63 

3.34 

,B6? 

878 

24.1 

32.1 

13.2 

.66 

3.81 

.899 

930 

25.4 

35.1 

13.4 

.69 

4.46 

.921 

970 

DISCUSSION 


Load-Deflection  Relationship 

Results  taken  from  8,  12,  and  16  ply  specimens  have  been  plotted  in  Figure 
3.  While  PL’ /El  can  be  predicted  by  the  analysis,  as  shown  by  the  solid  line, 
the  values  may  also  be  determined  directly  from  measured  quantities.  The  ac¬ 
curacy  of  the  analysis  is  clearly  demonstrated  by  the  fact  that  the  experimental 
data  points  for  all  three  thicknesses  fall  nicely  on  the  curve  predicted  from  the 
analysis  when  a  constant,  average  value  of  El  for  each  laminate  is  used. 
Another  way  of  demonstrating  the  accuracy  of  the  analysis  is  to  predict  the 
crack  length  L  at  various  times  in  a  test  using  the  measured  values  of  P  and  2A 
and  the  average  El.  The  predicted  values  of  L  are  compared  with  the  ex¬ 
perimental  results  in  Table  3  for  a  twelve-ply  specimen.  The  agreement  is  seen 
to  be  excellent,  due  in  large  part  to  the  insensitivity  of  L  to  the  experimentally 
determined  parameters. 

This  method  of  determining  the  instantaneous  crack  length  should  be  very 
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Table  ).  Companion  of  Measured  Values  of  Crack  Length,  L,  To  Value 
Calculated  Using  Nonlinear  Beam  Analysis  and  Measured  Values  of  Crip 
Displacement,  its,  and  Load,  P. 


Measured  Crack  Length  (an)  Calculated  Crack  Length  (cm) 


8.89 


10.16 

11.43 

12.70 

13.97 

15.24 

16.51 

17.78 

19.05 

20.32 

21.59 

22.86 

24.13 

25.40 


8.91 
10.13 
11 .48 
12  70 
13.97 
15.21 

16.54 
17.70 
19.05 
20.35 
21.56 
22.91 
24.25 

25.55 


useful  in  tests  of  opaque  materials,  such  as  graphite/epoxy  laminates. 
However,  for  high  accuracy  the  value  of  El  employed  in  the  prediction  of 
crack  length  should  be  obtained  directly  from  the  laminate  using  one  or  more 
pre-determined  lengths  and  theory  (i.e.  second  and  third  columns  of  Table  1). 
This  is  necessitated  by  the  fact  that  predicted  values  of  El  differ  somewhat 
from  the  measurements,  possibly  due  to  nonuniformity  in  the  fiber  distribu¬ 
tion.  Further,  the  photomicrographs  of  the  twelve  ply  laminate  in  Figure  8 
show  that  the  local  distribution  is  quite  heterogeneous.  In  many  cases  the 
fibers  appear  to  be  in  contact,  which  could  be  a  significant  factor  in  producing 
the  scatter  in  the  delamination  fracture  energy  discussed  next. 

Energy  Release  Rate-Crack  Speed  Relationship 

The  energy  release  rate  is  plotted  against  crack  speed  in  Figure  9.  The  data 
can  be  approximated  by  a  power  law  G  ~  L° '  Schapery  ( 1 2]  has  predicted  that 
crack  growth  rates  in  viscoelastic  material  m lj  be  described  by  an  equation  of 
the  form  L  ~  Kj  where  q  =  2(1  +  1/m)  if  the  intrinsic  fracture  energy  and 
strength  at  the  crack  tip  are  independent  of  crack  speed  and  the  creep  com¬ 
pliance,  D,  obeys  the  power  law  D  ~  tm.  For  many  glassy  polymers,  including 
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the  epoxy  resin  in  Scotchply,  the  creep  compliance 
is  of  the  form  D  =  D0  +  £>,/"  [13].  The  value  of 
D,  increases  with  increased  molecular  mobility, 
such  as  may  result  from  the  increase  in  volume  due 
to  the  high  triaxial  stresses  near  a  crack  tip.  If  D,  is 
sufficiently  large  then  m  =*  n.  Inasmuch  as  K,  ~ 
G"  (cf.  Equation  (16))  and  G  ~  L° we  obtain  q 
—  20  and  n  ss  0. 1  if  m  cs  n.  In  order  to  check  this 
value  for  n,  creep/recovery  tests  were  run  on 
unidirectional  tensile  specimens  of  Scotchply  with 
a  90°  Tiber  angle  to  determine  the  time  dependence 
of  the  matrix.  An  exponent  of  n  ot  0.05  —  0.10 
was  determined  for  these  tests,  which  is  consistent 
with  the  fracture  results  for  delamination  fracture 
toughness.  However,  the  data  are  so  limited  that 
one  should  consider  the  relationship  q  ^ r  2(1  + 
1/n)  as  very  tentative. 

Fracture  toughness  of  unidirectional  Scotchply 
has  been  measured  by  Wu  [14]  using  a  center- 
notched  specimen  with  90°  fibers.  Averaging  out  rate  effects  in  order  to  obtain 
a  single  critical  stress  intensity  factor,  he  obtained  for  Kk  a  value  of  1.9 
MPa\f  m.  These  results  may  be  compared  to  our  results  if  G  values  are  express¬ 
ed  in  terms  of  K ,  using  the  relationship 


Figure  8.  Photomicrographs  of 
Scotchply  showing  distribution 
of  glass  fibers  in  epoxy  matrix: 
(!)  J 2  plies,  500x  ;  (b)  12  plies, 
lOOx 


K,  =  (EG)" 


06) 


where  E  is  the  effective  modulus  and  can  be  shown  to  be  approximately  equal 
to  one  half  the  transverse  modulus  E> ,  for  an  orthotropic  material  [15].  Utiliz¬ 
ing  the  manufacturer’s  predicted  value  for  E„  of  9.7  x  )  O' MPa  and  our 
observed  range  of  G  values  of  525  —  1000  N/m  gives  a  predicted  range  of 
values  for  A',  of  1  6  —  2.2  MPasf  m  which  brackets  Wu’s  result  for  Krc.  These 
values  may  also  be  compared  with  values  measured  for  metals  and  various 
plastics,  as  summarized  in  Table  4. 

Several  simplifying  assumptions  were  made  implicitly  in  the  analysis  that 
appear  to  be  justified  by  the  good  agreement  between  the  measurements  and 
predictions  from  the  nonlinear  beam  analysis,  as  seen  in  Figure  3.  These 
assumptions-tnclude  plane  stress  and  “rigid  wall”  behavior  at  the  beam  end 
point,  which  in  our  specimen  is  the  crack  front.  The  actual  behavior  in  this 
crack  tip  region  may  only  be  described  by  a  more  complex  three  dimensional 
analysis  using  finite  elements;  however,  the  departure  from  these  idealizations 
appears  to  have  a  negligible  effect  on  the  overall  load-deformation  relation  (cf. 
Figure  3),  and  consequently  a  small  effect  on  strain  energy  stored,  in  that  it 
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Figure  9.  Energy  release  rate  versus  crack  speed  for  three  different  thicknesses  of  Scotchply. 


Table  4.  Typical  Fracture  Toughness  Parameters  IK/c  and  CqI 


H*ter(»l 


1C 

HP..S 


5c 

*H/m 


HtUli  (5.  101 
T1-6JM-4V 
7075- T651 
4340 
?0?4  T3 


Tjjgnjoplij t t  c_ jjittrlih  HI.  1?1 


Polywethyl  Me the cry late 

1.6  - 

1.9 

1.1 

Polystyrene 

0.98  . 

1.1 

0.35 

Polyvinyl  Chloride 

1.6  - 

2.3 

1.2 

Nylon  -  6,  6 

0.51  - 

0  83 

0.25 

Poly*t6yl»n» 


us 

171 

24 

8 

60  -  99 

18  -  49 

44 

V 

0  83  -  l.Z 


5.0 
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may  effectively  be  neglected.  Significant  breakage  of  glass  fibers  away  froM 
the  crack  plane  could  also  have  introduced  an  error  into  the  bending  rigidity 
and  energy  release  rate.  Sample  calculations  indicate  that  the  axial  stress  in  thtl 
beam  never  approached  the  failure  value;  furthermore,  almost  no  fiber  frac¬ 
ture  was  noted  in  the  tests  other  than  along  the  crack  plane.  Finally,  as  a  result 
of  transverse  strains  (normal  to  the  loading  direction),  the  beam  has  a  conn 
pound  curvature  which  would  tend  to  stiffen  the  beam  and  depend  or  beam 
length.  Again,  in  view  of  the  agreement  in  Figure  3,  this  effect  was  apparently 
negligible  for  the  specimens  used. 

CONCLUSIONS 

A  simple  approach  to  the  determination  of  delamination  fracture  toughness 
in  the  opening  mode  has  been  developed  using  a  split  beam  and  a  nonlinear 
analysis.  The  analysis  has  been  confirmed  with  experimental  measurements  on 
Scotchply  using  specimens  of  three  thicknesses,  tested  over  a  wide  range  of 
crack  growth  rates.  The  measured  range  of  energy  release  rates  of  525—1000 
N/m  for  the  range  of  crack  growth  rates  studied  is  shown  to  be  consistent  with 
predictions  from  an  idealized  viscoelastic  crack  growth  theory  and  viscoelastic 
behavior  of  the  resin  as  determined  from  creep/recovery  tests. 
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A _ _ _ 

ABSTRACT 

Pertinent  results  front  a  generalized  J  Integral  theory  for  nonlinear  visco¬ 
elastic  media  are  first  reviewed,  and  then  some  special  cases  for  power-law  ma¬ 
terials  are  given  to  illustrate  their  application  in  crack  growth  and  failure 
analysis.  This  theory  was  recently  developed  by  the  author;  it  is  a  generaliza¬ 
tion  of  the  familiar  J  integral  theory  which  is  customarily  restricted  to  non¬ 
linear  elastic  and  viscous  materials  and  to  materials  obeying  the  classical  de¬ 
formation  theory  of  plasticity.  Some  experimental  results  from  cyclic  and  tran¬ 
sient  loading  of  fibrous  and  particulate  composites  are  then  given  and,  together 
with  the  theory,  are  used  in  a  tentative  interpretation  of  the  viscoelastic 
fracture  process. 


notch  or  crack  size 
crack  speed 

time-scale  shift  factor  for  temperature  and  moisture  effects 

unspecified  constants 

creep  compliance  for  uniaxial  loading 

superscript  denoting  pseudo-elastic  variables 

reference  modulus 

indices  corresponding  to  coordinate  directions  (1,2,3) 

J  Integrals  for  elastic,  viscous,  and  viscoelastic  media,  respec¬ 
tively 

stress  Intensity  factor  for  opening  mode 
exponent  relating  crack  speed  and  Jv 
vector  normal  to  J  integral  path 

exponents  defining  time-dependence  and  nonlinearity,  respectively 
potential  energy  of  body  and  applied  loads 
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designates  location  of  mathematical  crack  tip 

exponent  relating  crack  speed  and  applied  stress 

surface  tractions 

displacements 

energy  density 

Cartesian  coordinates 

length  of  failure  zone 

path  used  to  evaluate  J  Integral 

intrinsic  fracture  energy 

strains  and  strain  rates,  respectively 

stresses 

magnitude  of  normal  stress  at  crack  tip 
constant  (yield  stress  if  N<<1) 
time  variable  of  integration 
pseudo  strain  energy  density 
pseudo  complementary  energy  density 


INTRODUCTION 

Time-dependent  deformation  behavior  of  the  matrix  constituent  in  many  com¬ 
posites  is  often  very  pronounced,  especially  in  elevated  temperature  environ¬ 
ments;  exposure  of  polymer-matrix  composites  to  high  relative  humidity  leads  to 
similar  behavior  and  aggravates  the  effect  of  temperature  if  sufficient  time  is 
provided  for  absorption  of  the  water  vapor.  This  time-dependence,  or  visco¬ 
elasticity,  is  an  important  factor  in  determining  the  rate  of  growth  of  micro- 
and  macrocracks  (Including  delamination  cracks).  Complicating  the  composite 
response  characteristics  is  nonlinearity  on  the  macroscale  (although  often  due 
in  part  to  matrix  microcracking)  and  on  the  microscales  of  the  reinforcing  fibers 
or  particles  and  the  crack  tip  structure.  Development  of  realistic  theoretical 
models  ">f  damage  growth  and  failure  therefore  requires  in  many  cases  that  one 
account  jor  both  time-dependence  and  nonlinearity  in  the  deformation  behavior. 

In  spite  of  the  complexity  of  the  problem,  it  is  believed  possible  to  develop 
at  this  time  practical  models  which  will  be  useful  in  aiding  our  understanding 
of  structural  response  of  composites  so  that  Improved  materials,  structural  de¬ 
sign  methods,  and  reliability  can  be  realized. 

In  this  paper  we  briefly  describe  and  then  apply  a  new  fracture  theory  (1! ) , 
based  on  the  so-called  J  integral  (2),  for  analyzing  micro-  and  macrocracking  in 
nonlinear  viscoelastic  composites. 

Following  (2),  the  J  integral  is  defined  for  two-dimensional  problems  by 
r 

J  -^.(Wdx^  — da)  (1) 

where,  as  shown  in  Fig.l,  I'  is  a  curve  surrounding  the  notch  tip,  and  repeated 
indices  imply  suimnation  over  their  range.  The  quantity  W»W(e  )  is  defined  such 
that  J 


aij  ’  3W/3cij 


(2) 


where  o  and  e  are  the  stress  and  strain  tensors,  respectively.  Thus,  W  is 
the  strHn  energy  density  or  an  analogous  potential  function  in  the  deformation 
plasticity  theory.  Also,  T  is  the  traction  vector,  in  which  The  dis¬ 

placement  vector  is  denoted  oy  u^.  A  very  useful  feature  of  J  19  tnsf  its  value 


is  independent  of  path  f  if  W  does  not  depend  on  the  coordinate  x.  (cf.  Fig.l)  ~ 
other  than  through  the  dependence  of  e  on  x^  the  material  may  be  anisotropic 
and  have  properties  that  vary  with  x,.  :A  recond  Important  characteristic,  espe¬ 
cially  for  experimental  characterization  of  fracture,  is  that  if  the  notch  tip  is 
advanced  an  amount  da  without  change  in  tip  structure, 

J  -  -3P/5a.  (3) 

which  is  the  rate  of  decrease  of  potential  energy  (of  the  body  and  applied  loads) 
with  respect  to  notch  length  per  unit  thickness  normal  to  the  page. 


Fig.l  Crack  or  notch  in  elastic  material,  with  tip  at  Pt;  highly  damaged  material 
is  shaded. 

Over  ten  years  ago.  Rice  (2)  developed  and  applied  the  J  integral  for  two- 
dimensional  deformation  fields  (plane  strain,  generalized  plane  stress,  and 
antiplane  strain)  in  nonlinear  elastic  media  with  notches  or  cracks.  This 
work,  which  was  theoretically  limited  to  small  strains,  lead  the  way  to  major 
developments  during  the  past  decade  in  the  engineering  characterization  and 
analysis  of  fracture  initiation  in  metals  with  both  small  scale  and  large 
scale  plastic  deformations  (e.g.,  it)  •  (Although  the  J  integral  was  seeming¬ 
ly  developed  for  elastic  materials,  the  constitutive  equation  for  the  deforma¬ 
tion  theory  of  plasticity  without  unloading  is  anlogous  to  that  for  elastic 
materials  in  the  sense  that  it  may  be  written  as  in  Eq.(2);  consequently,  ap¬ 
plication  to  fracture  initiation  in  ductile  metals  has  been  successful.) 

The  J  integral  can  be  interpreted  as  a  conservation  law  for  mechanical 
energy.  As  such,  one  would  expect  to  be  able  to  deduce  a  three-dimensional 
generalization  with  large  strains  if  suitable  measures  of  stress  and  strain 
are  used.  Indeed,  such  a  version  does  exist;  as  noted  by  Knowles  and  Stern¬ 
berg  (5),  in  their  study  of  finite  elastic  deformations  near  a  crack  tip, 
this  conservation  law  is  implicit  in  an  early  result  due  to  Eshelby  (6)  in 
the  theory  of  continuous  dislocations.  The  J  Integral  with  large  plastic 
deformations  was  used  by  McMeeking  (7)  in  a  study  of  phenomena  at  crack 
tips,  including  void  growth.  Rice  (8)  has  reviewed  applications  to  stable 
crack  growth  in  rate-independent  media. 

An  integral,  say  J',vlth  properties  analogous  to  those  for  J  may  be 
defined  for  nonlinear  viscous  media  when  the  constitutive  equation  can  be 
expressed  in  the  form 

<ri}  -  3H’/9ci3  (4) 

where  W'»w'(e  ,)  and  e  is  strain  rate.  The  same  definition  as  in  Eq.(l)  is 
used  for  J’,  except  strain  rates  and  velocities  replace  strains  and  displace¬ 
ments.  Landes  and  Begley  (9)  showed  that  crack  speed  a  in  two  different 
specimen  geometries  of  superalloy  at  high  temperature  may  be  correlated  in 
terras  of  J';  viz.,  a»a(J').  Sakata  and  Finnie  (10)  and  Ohji.  et  al.  (11) 
independently  studied  stress  intensity  factor  K.  and  J',  respectively,  as 
fracture  parameters  for  nonlinear  viscous  materials.  Many  papers  have  now 
appeared  which  provide  strong  experimental  support  for  the  hypothesis  that 


a=a(J')  for  creeping  metals  at  high  temperatures  under  constant  and  variable 
loading,  and  the  reader  is  referred  to  (12-14)  for  recent  work  and  citations  of 
other  studies. 

In  order  to  provide  a  rigorous  theoretical  argument  that  J'  is  the  appropri¬ 
ate  basic  parameter  for  correlating  fracture  data,  and  to  relate  observed  behavior 
to  basic  material  properties,  it  is  believed  necessary  to  investigate  crack-tip 
behavior  under  LocaLly  finite  stresses  in  which  a  finite  scale  for  damage  and 
material  separation  exists  (e.g.,  the  length  a  in  Fig.l.).  To  the  author's  know¬ 
ledge,  previous  theoretical  studies  involving  nonlinear  viscous  or  viscoelastic 
behavior  outside  of  the  failure  zone  have  introduced  strain  rates  which  are 
infinite  at  the  crack  tip  for  mathematically  sharp  cracks.  The  basic  shortcoming 
of  such  singular  fields  can  be  illustrated  by  analysis  of  linear  viscoelastic 
media,  as  reviewed  by  Rice  ( 8) .  For  example,  starting  with  a  solution  for  finite 
stresses,  Schapery  (15)  proved  chat  if  the  magnitude  of  the  crack-tip  stress  is 
mathematically  increased  without  limit,  the  scale  a  vanishes  and  crack  growth  be¬ 
comes  independent  of  the  linear  viscoelastic  properties;  only  the  initial  value- 
of  compliance  or  modulus  remains.  Associated  infinite  9train  rates  produce  the 
physically  unacceptable  prediction.  Numerous  publications  on  the  analysis  of 
crack  growth  in  linear  viscoelastic  media,  based  on  a  finite  crack-tip  stress, 
now  exist  (e.g.,  see  (8)  and  (1S-20) ) . 

Very  recently,  Schapery  (_2)  extended  the  J  integral  concept  to  nonlinear 
viscoelastic  media  for  both  initiation  and  growth  of  cracks.  This  work  is  based 
on  finite  crack  tip  stresses,  so  that  one  may  derive  basic  failure  behavior  re¬ 
lated  to  the  physical  characteristics  which  define  the  zone  of  intensive  damage 
at  the  crack  tip.  This  region  is  called  the  "failure  zone"  (whose  length  is  a 
in  Fig.l)  rather  than  damage  zone;  the  theory  in  (1.)  allows  for  certain  kinds  of 
damage  in  the  nonlinear  material  surrounding  the  crack  tip,  such  as  the  type  of 
stable  microcrncking  and  void  growth  that  develops  in  many  composites  if  the 
stress  level  is  not  too  high.  For  propagating  cracks,  a  layer  of  c he  highly 
damaged  material  that  was  once  in  the  failure  zone  will  be  left  on  the  crack 
faces.  If  was  assumed  in  (1)  that  this  layer  is  thin  enough  to  be  neglected 
in  predicting  stresses  and  strains  in  the  neighborhood  of  the  crack  tip;  thus, 
a  is  the  only  scale  parameter  that  reflects  the  local  failure  process  for  pro¬ 
pagating  cracks  (as  well  as  that  for  cracks  which  just  start  to  grow,  corre¬ 
sponding  to  the  fracture  Initiation  condition). 

In  the  present  paper  we  first  give  the  constitutive  equation  used  in 
developing  the  generalized  J  integral  theory.  Discussed  next  is  the  relation 
of  crack  speed  to  the  generalized  J  Ir.egral,  designated  by  for  a  visco¬ 
elastic  material.  An  explicit  relation  between  crack  speed  and  applied  stress 
is  then  stated  for  a  single  crack  in  a  power  law  material;  both  time-dependence 
and  stress-dependent  nonlinearity  are  assumed  to  obey  power  laws.  The  two 
independent  exponents  combine  to  yield  one  exponent  defining  the  crack  speed 
as  a  power  law  in  stress;  for  fatigue  loading  of  monolithic  materials  and 
certain  types  of  composites,  the  same  exponent  relates  number  of  cycles  to 
failure  and  applied  stress.  When  used  with  theory,  the  values  of  experimen¬ 
tally  determined  exponents  are  very  helpful  in  identifying  primary  physical 
mechanisms  affecting  crack  growth  and  failure.  This  point  is  illustrated  by 
using  data  on  fatigue  failure  of  graphite  fiber-reinforced  epoxy  laminates 
and  macrocrack  propagation  in  composite  solid  propellant  (which  is  rubber 
with  a  high  volume  fraction  of  hard  particles). 

NONLINEAR  CONSTITUTIVE  EQUATION 

The  relation  between  stress  and  strain  tensors  is  assumed  in  che  form 
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where  the  e*  are  termed  pseudo-9trains ,  and  are  functions  of  stresses  (as  veil 
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on  the  particular  material)  through  a 


(6) 


If  Eq. (6)  can  be  inverted  to  yiqld  o  as  a  function  of  , 
that  a  potential  function  exists  with  the  property 


then  Eq.(6)  implies 
that 


♦-•c  +  01jeUT.;  .  _  (8) 

Viscoelastic  behavior  is  defined  by  the  creep  compliance,  D.  It  may  be  observed 
thac  the  notation  D(t,x)  is  often  used  instead  of  D(t-x,t);  but  the  notation 
is  equivalent  and  allows  for  aging.  Thermorheologically  simple  behavior,  a  form 
of  temperature  dependence  commonly  obeyed  by  polymers,  is  a  special  case  of  this 
representation  for  creep  compliance,  with  constant  or  transient  temperature.  The 
modulus  is  a  free  parameter  which  makes  Eq.(5)  dimensionally  correct  as  $  has 
the  dimensions  of  modulus;  may  be  selected  as  desired  to  simplify  the  notation 
in  fractura  results.  That  Eq.(5)  is  a  good  approximation  for  many  nonlinear  vis¬ 
coelastic  materials  is  discussed  in  ( 1) ;  also,  it  is  shown  in  (1)  that  may  de¬ 
pend  on  certain  important  types  of  damage,  such  as  distributed  microcracfes. 

If  we  assume  D(t-t,t)  is  proportional  to  t-x,  Eq.(5)  reduces  the  constitutive 
equation  for  a  linear  o*r  nonlinear  viscous  material. 

COPLAN AR  CRACK  GROWTH  IN  POWER-LAW  MATERIALS  _ _ _ _ 

By  means.,  of .  an  .elastic-viscoelastic  correspondence  principle  established  in 
(l5  and  approximations  similar  to  those  used  in  (16 ,  Part  II)  for  linear  theory,  it 
is  found  that  if  a  is  small  compare d~td_other  relevant  dimensions  such  as-crack  lengthT- 
2r  -  ERD(a/3a,t)Jv  (9) 

where  I\  the  intrinsic  fracture  energy,  is  the  mechanical  work  that  the  continuum 
does  on  the  failure  zone  when  one  unit  of  new  surface  is  formed.  The  generalized 
integral,  J  ,  is  defined  as  in  Eq.(l),  but  ♦  replaces  W  and  pseudo-displacements, 
u^  (corresponding  to  pseudo  strains,  Eq.(6))  replace  u  .  In  deriving  Eq.(9),  it 
is  assumed  crack  speed  is  essentially  constant  during  the  time  a/a,  which  is  that 
required  for  crack  growth  of  an  amount  a.  The  time  t  in  the  argument  of  D  re¬ 
flects,  for  example,  aging  and  transient  temperature  and/or  moisture  content;  but 
such  effects  are  assumed  constant  in  the  time  interval  a/a.  For  simplicity, 
these  transient  effects  will  be  neglected  in  all  subsequent  discussion. 

Let  us  now  introduce  a  power  law  for  creep  compliance, 

D(t-t)  -  D1(C-r)n  (10) 

where  D  and  n  are  positive  constants.  Also,  assume  $  is  a  homogeneous  func¬ 
tion  of^degree  N+l  in  pseudo  strain. 


«(ce 


|  •(Cjj) 


where  c  and  N  are  constants.  The  material  defined  by  $  may  be  anisotropic  and 
physically  nonhomogeneous .  ! 

In  order  to  illustrate  the  significance  of  the  exponent  in  Eq.(ll),  consider 
a  uniformly  Stressed  specimen  under  a  time-dependent  uniaxial  stress  o  (or,  more 
generally,  under  raultiaxial  proportional  loading).  We  find  the  strain  in  the  di¬ 
rection  of  loading  to  be 


,erdiJ L(c't)" 


d(o/o  ) 


where  a  is  a  constant.  If  N  is  small,  o  can  be  interpreted  as  a  yield  stress. 
In  a  creep  test  a  is  constant  and  Eq .(12)  reduces  to 


W  (a/flo> 


(13) 


Let  us  introduce  temperature  T  and  moisture  M  effects  in  standard  notation  for 
poLyrrers  by  writing 


Vl  '  a™<  <M) 

where  a^-a^CT.M)  is  che  tine-scale  shift  factor.  Thus,  the  creep  3train  is, 

(15) 


c>  (t/ara)n(o/0o)1/N 


The  requirement  of  a  finite  crack  tip  stress,  whose  magnitude  at  the  point 
PT  in  Fig.X  is  denoted  by  o  ,  lead9  to  the  relation  (1) , 
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(16) 


where  is  a  dimensionless  constant. 

In  general  o  and  r  may  depend  on  crack  speed.  If,  however,  they  are  con¬ 
stant,  Eqs.(9),  (lO) ,  (14),  and  (16)  yield  simply 
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where  C  and  k  are  positive  constants,  with 
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T~and  d  (instead  of- tr )  are  constant, ~ 
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(19) 


Finally,  if  a  and  crack  opening  displacement  Av^  (st  the  left  end  of  the  failure 
zone  in  Fig.l)  are  constant  (1), 


k  •  — - 

n(l  +  S) 


(20) 


The  constant  C  is  different  for  each  of  these  cases.  It  should  be  added  that  if 
the  spacewise  distribution  of  the  3tress  in  the  failure  zone  varies  greatly  with 
speed,  Che  change  will  be  reflected  as  speed  dependence  of  C;  but  this  effect 
will  be  neglected  here.  Also,  the  assumption  of  constant  ff  and  Av^  results  in 
Eq . (18) ,  and  therefore  does  not  represent  another  Independent  case. 

For  bodies  on  which  surface  tractions  are  specified,  instead  of  displace¬ 
ments,  J  as  a  function  of  these  tractions  is  the  same  as  for  an  elastic  material 
For  an  isolated  crack  in  a  macroscopically  homogeneous,  nonlinear  power  law  body 
under  proportional  loading  (  !_) ,  (assuming  a«a), 

i+  (i/8)  I 

Jv  ‘  aC2°o(~)  •  0  >  0  (21): 

o 


where  isAa  dimensionless  constant  and  o»o(t)  is  the  remote  stress  in  a  given 
direction;  N  is  for  the  far-field,  which  may  differ  from  N  for  the  near-tip  material 
Equations  (17)  and  (21)  yield 


a  «  C 


(°_)q 
3  a-,  \y  ; 
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(22) 


where  is  another  constant  and 

q  =  k(l  +  N)/N  (23) 

When  the  geometry  of  a  body  is  that  in  Fig.  2,  which  is  a  sheet  (or  thick 
slab)  that  is  loaded  vertically  through  rigid  clamps,  the  form  of  J  is  the  same 
as  in  Eq.  (21),  except  h  replaces  a.  Assuming  a<<a<h/2,  the  coefficient  is 
essentially  constant ,  where  o  Is  equal  to  the  stress  o  in  the  uniformly  stressed 
portion  of  the  strip;  as  a  good  approximation,  at  leagTT  if  N*l,  a  may  be  taken  a9 


(24) 


the  total  vertical  load  Pv  divided  by  the  uncracVed  area  (1).  Thus 
i_k 
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Given  o(t)  or  Pv(t),  Eqs. (22)  and  (24)  may  be  integrated  to  predict  crack 
size.  If  the  stress  is  applied  for  a  long  enough  period,  a  failure  time  t  is 
found.  Suppose  k>l  in  Eq.(22);  then  a-~  when  f 
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(25) 
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where  C,.  is  a  constant  proportional  to  a  ;a  is  initial  crack  size.  The  same 
type  of  result  is  derived  from  Eq.(24)  b8t,C  ois  to  be  replaced  by  E  and 
failure  occurs  when  2a-L.  It  should  be  added  that  Eqs. (21)-(25)  are  not  restrict¬ 
ed  to  plane  stress  or  plane  strain  problems*  The  axisymmetric  problem  of  a  pen- 
ney-shaped  crack  is  also  included,  in  which  2a  and  L  are  crack  and  sample 
diameters,  respectively.  r 

When  the  applied  load  is  constant  (creep  test)  and  the  physical  environment 
is  constant  (constant  value  of  a_.) .  Eq.(25)  yields  t,-v0“q.  Similarly  for  a 
cyclic  tensile  fatigue  test,  N  Z  o"q,  where  N  is  tfce  number  of  cycles  to  fail¬ 
ure  and  o  is  stress  amplitude;  it  is  assumed  for  simplicity  here  that  the  wave 
Shape  and  amplitudedo  not  change  with  time.  Thus,  at  least  for  the  two  idealized 

problems  considered  above,  there  is  a  simple  relationship  between- the  ToiT-Tog - 

sl^p,e.<Urom« crack  speed,  creep-rupture,  or  fatigue  data  and  the  basic  material 
exponents ,  as  given  by  Eq.(23). 

^We  have  compared  the  value  of  the  exponent  q  found' from  crack  growth  and - ■ 

failure  data  in  tests  conducted  on  different  polymeric  materials.  It  is  often 
found  that  the  value  corresponding  to  k  in  Eq.(18)  agrees  with  the  data  Some 
examples  are  dlscussec  in  the  next  section.  In  contrast,  for  metals  undergoing 
•viscous  creep  (n-1) ,  Eq.(20)  appears  to  provide  the  best  agreement  (1);  typically 
for  metals  N<<1,  and  thus  there  is  little  difference  between  Eqs.(19y  and  (20). 
EXAMPLES  AND  CONCLUDING  OBSERVATIONS 


The  generalized  J  theory  may  be  employed  in  different  ways  to  characterize 
and  predict  failure  behavior  of  materials  (1).  Also,  as  noted  above.  It  provides 
a  direct  relation  between  fracture  response  and  basic  creep  characteristics-  only 
this  aspect  will  be  examined  here.  ’  7 

Crack  speed  and  overall  specimen  fracture  was  studied  in  (16,  Part  III)  in 
which  the  linear  theory  was  applied  to  a  polyurethene  rubber.  For  the  material 
used,  n-0.5,  N-1,  and  Eqs. (18)  and  (23)  predict  q-6;  this  value  of  q  is  in  ex¬ 
cellent  agreement  with  the  data,  providing  good  evidence  that  r  and  a  are  con¬ 
stant.  The  theory  in  (16)  did  not  account  for  the  local  nonlinearity®  as  defined 
by  N;  as  shown  here,  even  though  this  exponent  may  not  be  unity,  it  has  no  effect 
on  q  (unless,  of  course,  Eq.(20)  were  to  apply).  i 

Application  of  the  theory  to  composite  materials  is  not  necessarily  straight- 
fo-ward  and.  In  fact,  considerable  additional  analysis  may  be  required  to  relate 
overall  loads  and  deformations  Co  J  for  microscale  phenomena  governed  by  Eq.  (9) 
and  more  specifically  by  Eq.(17).  This  complexity  is  due  in  part  to  the  fact  that 
the  creep  compliance  D  in  Eq.(5)  for  the  neighborhood  of  cTack  tips  may  be  dif- 
ferentfor  the  macroscale  of  a  composite  material.  On  the  other  hand,  if  the 
matrix  constituent  of  a  composite  is  homogeneous  and  obeys  Eq.(5),  and  the  rein¬ 
forcement  phase  is  essentially  rigid,  the  creep  compliance  will  ba  the  same  re¬ 
gardless  of  the  scale  of  the  deformation  process.  Stress  and  environment-induced 
phase  changes  (e.g.,  glass-to-rubber  transition  at  crack  tips)  introduces 
additional  complexity.  Temporarily  setting  aside  such  concerns,  let  us  consider 
the  data  in  Figs. 3  and  4.  For  the  fatigue  loading  case.  Fig. 3,  and  using  either  a 
dominant  flaw  model,  Eq.(25),  or  a  matrix  degradation  model  with  uniformly  distri¬ 
buted  microflaws  (23)  (each  of  which  obeys  Eq.(25))  we  anticipate  that  q-2(l+l/n); 
it  is  assumed  that  any  global  nonlinearity  is  due  to  microcracking,  and  therefore  ' 
N-1.  Inasmuch  as  the  zero-time  (glassy)  compliance  is  neglected  in  this  compari¬ 
son,  and  the  value  of  q  Is  based  on  k  in  Eq.(18),  we  suggest  that  the  cracking 
process  is  controlled  by  polymer  matrix  above  the  glass  transition  temperature;  the 
—very  high,  stresses  in  the  matrix  around  microcracks. are  thus.assumed  to  lower  the  _ 


glass  transition  temperature  below  the  local  temperature  (which  may  be  elevated  due 
to  heat  generation)^)*  For  a  fiber-dominated  laminate  the  value  of  q  is  larger 
than  in  Fig. 3  (21,  Phase  III).  As  indicated  in  Fig. 4,  q-2(l+l/n)  for  the  solid 
propellant  if  the  data  are  normalized  to  constant  strain  level  through  the  func¬ 
tion  f;  data  actually  obtained  in  constant  strain  tests  obey  q»2(l+l/n)  (22) .  It 
is  believed  the  interaction  of  far-field  microcracking  with  the  macrocrack  causes 
this  strain-dependence  (1), 

Clearly,  there  are  many  factors  that  have  to  be  accounted  for  in  the  devel¬ 
opment  of  models  for  micro-  and  macro-cracking  of  composites.  But,  it  is  believed 
that  the  fracture  theory  in  this  paper  will  prove  useful  in  a  systematic  mechanics 
approach  in  view  of  the  present  encouraging  results  and  the  broader  theoretical 
basis  developed  in  (1) . 
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Fig. 2  A  cracked  strip  or  slab.  Also,  an  elementary  model  for  raicrocracking  in 
the  soft  matrix  phase  of  a  composite. 


2  >  « 


LOG  NUMBER  OF  CYCLES 

Fig. 3  The  S-N  curve  of  [+45/902 ]g  tensile  fatigue  coupons  of  graphite/ 
epoxy  (AS/3501-6);  creep  exponent  from  (21^  Phase  II).  Specimen  mois¬ 
ture  and  temperature  levels  are  in  equilibrium  with  indicated  environ¬ 
ments.  Frequency  ■  3  HZ,  stress  ratio  ■  0.1.  After  (21,  Phase  III). 
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Fig. 4  Crack  propagation  in  strip-biaxial  specimens  of  an  HTPB  propellant  (e  is 
in  percent  and  cross-head  rate  is  constant  and  given  in  Inches/aln) . 

After  (22). 


